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A.^ESCRIPTION  OF  THE  SCIENTIFIC  RESEARCH  GOALS 

The  thrust  of  this  project  is  to  come  up  with  the  design  and  specification  of  advanced 
composite  materials  using  dielectric  ceramics,  chiral  polymers  or  piezoelectric  materials  as 
inclusions  with  polymer-based  materials  as  the  binders.  The  research  is  theoretical  and 
complements  the  associated  experimental  and  developmental  programs  on  such  materials  at  the 
Research  Center  for  the  Engineering  of  Electronic  and  Acoustic  Materials  and  the  Materials 
Research  Laboratory.  Basic  research  on  the  interaction  of  electromagnetic  and  acoustic  fields  with 
composite  materials  at  wavelengths  comparable  to  inclusion  size  and  spacing  were  to  be  studied  to 
come  up  with  structure-  property  relations  that  offer  better  possibilities  for  controling  material 
behavior  than  those  permitted  by  simple  mixing  rules.  Emphasis  was  on  chiral  composites, 
piezoelectric  composites,  fractal  composites,  and  ferrite  composites. 


B.  SIGNIFICANT  RESULTS  IN  THE  PAST  YEAR 

Much  progress  has  been  achieved  during  the  three  years  of  this  program  as  evidenced  by  the 
numerous  publications  in  refereed  journals. 

Basic  research  on  the  electromagnetics  of  chiral  media  have  now  been  applied  to  chiral 
waveguides,  chiral  mirrors  and  their  applications,  Green's  functions  for  chiral  media  which  were 
then  used  in  the  solution  of  several  problems  of  practical  interest.  The  first  monograph  on  this 
subject  entitled  "  Time  Harmonic  Electromagnetic  Fields  in  Chiral  Media  "  ,  co-authored  by 
Lakhtakia,  Varadan  and  Varadan  was  published  by  Springer  -  Verlag  in  1989.  Multiple  scattering 
calculations  have  been  extended  to  composites  containing  chiral  inclusions,  as  well  as  anisotropic 
inclusions.  We  have  demonstrated  that  chirality  significantly  alters  the  absorption  characteristics  of 
an  otherwise  low  loss  composite.  This  has  also  been  verified  by  experiments  on  chiral  composites 
designed  and  prepared  at  our  laboratory  and  measured  in  the  8  -  18  GHz  frequency  range.  The 
chirality  parameter  (3  has  been  measured  for  the  first  time  for  chiral  materials.  This  parameter  has 
not  been  measured  for  optically  active  materials,  although  that  phenomenon  has  been  known  and 
studied  for  over  one  hundred  years. 

Composites  containing  anisotropic  inclusions  have  been  modeled,  treating  for  the  first  time 
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the  effects  of  inclusion  anisotropy.  A  new  coupled  dipole  model  was  formulated  and  implemented 
numerically  to  compute  the  single  scatterer  response.  These  were  then  included  in  the  multiple 
scattering  computations.  The  results  were  then  compared  with  models  which  neglect  the  anisotropy 
of  the  inclusions.  The  differences  are  tremendous  and  wrong  approximations  can  lead  to  gross 
errors. 

The  phenomenon  of  back  scattering  enhancement  is  of  great  practical  interest  because  it  may 
permit  imaging  through  composite  media.  This  particular  problem  is  engaging  the  interest  of 
several  theorists  and  experimentalists.  Using  our  multiple  scattering  formulation,  we  have  been 
able  to  obtain  good  agreement  with  enhancement  observed  experimentally  on  random  systems  in 
the  optical  region. 

The  theoretical  problems  associated  with  the  microwave  sintering  of  ceramics  fit  very  well 
into  the  scope  of  our  multiple  scattering  studies.  Green  ceramics  in  the  form  of  porous  pellets  have 
a  high  dielectric  loss  tangent,  leading  to  absorption  of  microwave  power  and  resulting  increase  of 
temperature  and  loss  tangent.  This  creates  a  condition  for  rapid  increase  of  temperature  often 
referred  to  as  temperature  runaway.  By  theoretically  analyzing  the  effective  properties  of  the 
porous  ceramic,  the  thermal  runaway  condition  can  be  predicted  by  including  thermal  effects  via  a 
diffusion  equation  and  hot  wall  boundary  conditions  correlated  with  experimental  observations. 
Research  is  ongoing  in  this  area. 

Inverse  problems  associated  with  ferrite  composites  have  been  studies.  Ferrite  composites 
have  both  dielectric  and  magnetic  properties  which  are  complex  and  frequency  dependent.  From 
measurements  of  the  propagation  constant,  it  is  not  possible  to  infer  both  properties.  An  inverse 
procedure  has  been  proposed  and  has  been  tested  against  experimental  data.  It  can  be  used  to  infer 
constituent  properties  from  the  data.  This  is  a  convenient  way  of  obtaining  constituent  properties, 
especially  for  constituents  in  powder  form. 

C.  PLANS  FOR  NEXT  YEAR'S  RESEARCH 

This  is  the  final  report  for  the  three  year  program.  The  Center  will  continue  development  of 
further  applications  of  novel  chiral  composite  as  well  as  piezoelectric  composite  coatings  for 
various  types  of  sensors  and  coatings. 
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BRIEF  RESEARCH  REPORTS 
Electromagnetic  Response  of  Chiral  Media 

Electromagnetic  waves  in  a  chiral  medium  are  circularly  birefringent.  In  other  words,  left-  and 
right-circularly  polarized  waves,  with  different  phase  velocities,  exist  in  this  class  of  media.  In 
order  to  describe  the  electromagnetic  properties  of  isotropic  chiral  media,  the  usual  constitutive 
equations,  D  =  eE  and  B  =  jxH,  are  inadequate  because  they  admit  to  a  single  phase  velocity. 
Instead,  D  =  e  [E  +  PVxE]  and  B  =  ji  [H  +  pVxH];  these  constitutive  relations  are  symmetric 
under  time-reversality-  For  a  right-handed  medium,  the  chirality  parameter  p  >  0;  while  for  a 
left-handed  medium,  (5  <  0.  Several  theoretical  aspects  of  electromagnetic  wave  propagation  and 
scattering  in  isotropic  chiral  media  have  been  explored,  and  have  been  described  in  the  previous 
Annual  Report. 

A  major  activity  during  the  present  period  was  the  publication  of  a  book  on  electromagnetic 
chirality  in  the  Lecture  Notes  in  Physics  series  published  by  Springer-Verlag.  This  is  the  first  book 
on  the  subject,  and  it  aims  to  summarize  the  recent  developments  in  this  exciting  new  area. 
Undertook  as  a  service  to  the  electromagnetics  community,  this  book  lays  down  the  foundations  of 
electromagnetic  theory  as  applicable  to  chiral  media.  [A.  Lakhtakia,  V.K.  Varadan  and  v.v.  Varadan, 
Time-Harmonic  Electromagnetic  Fields  in  Chiral  Media,  Heidelberg;  Springer-Verlag,  1989.] 

In  general,  when  a  wave  propagating  in  a  homogeneous  isotropic  chiral  medium  encounters  a 
homogeneous  isotropic  chiral  scatterer,  both  the  scattered  and  the  internally  induced  fields  contain 
LCP  and  RCP  components  regardless  of  the  state  of  polarization  of  the  incident  field.  In  the 
previous  Annual  Report,  we  had  shown  that  the  reflected  and  the  transmitted  plane  waves  across  the 
planar  interface  of  a  mirror-conjugated  chiral  media  have  the  same  handedness  as  the  incident  plane 
wave.  It  has  now  been  proved  that  the  scattered  as  well  as  the  internal  fields  have  the  same  state  of 
polarization  as  the  incident  field,  if  the  scatterer  and  the  surrounding  medium  are 
impedance-matched;  this  conclusion  holds  regardless  of  the  two  chiral  parameters  involved  as  well 
as  the  geometry  of  the  scatterer.  [A.  Lakhtakia,  V.K.  Varadan  &  V.V.  Varadan,  'Influence  of  impedance 
mismatch  between  a  chiral  scatterer  and  the  surrounding  chiral  medium,'  J.  Mod.  Opt.  36, 1989, 1385  -  1392.] 

A  fully  vector  treatment  of  the  problem  of  scattering  by  singly  periodic  interfaces  between 
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chiral  and  achiral  media  was  developed.  Using  an  improved  T-matrix  method  of  solution,  this 
treatment  is  able  to  bypass  the  restrictions  imposed  by  the  Rayleigh  hypothesis.  In  addition,  several 
other  related  problems  have  their  solutions  embedded  in  the  method  developed.  The  specific 
features  of  the  solution  procedure  allow  it  be  easily  converted  for  doubly  periodic  interfaces.  The 
analysis  may  be  useful  for  developing  new  surface-relief  gratings  made  using  chiral  materials.  [A. 
Lakhtakia,  VJC.  Varadan  &  V.V.  Varadan,  'Scattering  by  periodic  achiral-chiral  interfaces,'  accepted  for  publication  in 
J.  Opt.  Soc.  Am.  A].  Furthermore,  the  excitation  of  a  planar  achiral-chiral  interface  by  near  fields  was 
also  examined.  [A.  Lakhtakia,  V.K.  Varadan  &  V.V.  Varadan,  'Excitation  of  a  planar  achiral/chiral  interface  by 
near  fields,'/.  Wave-Mater.  Interact.  3, 1988,  231  -241.] 

Wave  propagation  in  unidirectionally  inhomogeneous  chiral  media  was  examined,  i.e.,  the 
constitutive  properties  of  the  chiral  media  vary  along  the  z  axis.  It  was  also  assumed  that  the 
electromagnetic  field  was  a  function  of  the  z  coordinate  alone.  It  was  shown  that  the  electromagnetic 
field,  in  this  arrangement,  can  be  decomposed  into  two  mutually-independent  circularly-polarized 
states.  Coupled  first-order  differential  equations  were  derived  to  describe  the  fields,  and  various 
solution  procedures  discussed.  Extension  was  made  to  the  case  when  the  medium  inhomogeneity  is 
periodic,  for  which  case  a  perturbational  solution  was  explicitly  given.[A.  Lakhtakia,  V.K.  Varadan  & 
V.V.  Varadan,  'Propagation  along  the  direction  of  inhomogeneity  in  an  inhomogeneous  chiral  medium,’  accepted  for 
publication  in  Int.  J.  Engg.  5c/.] 

In  the  previous  Annual  Report,  we  had  described  spherical  chiral  resonators.  Continuing  on  the 
same  lines,  propagation  of  electromagnetic  waves  in  a  parallel-plate  waveguide  wholly  filled  with  a 
chiral  medium  was  examined.  The  dispersion  equation  derived  leads  to  two  sets  of  modes,  and  the 
propagation  constants  for  the  two  sets  were  numerically  obtained.  [V.K.  Varadan,  A.  Lakhtakia  &  V.V. 
Varadan,  ’Propagation  in  a  parallel-plate  waveguide  wholly  filled  with  a  chiral  medium,’  /.  Wave-Mater.  Interact.  3, 
1988, 267  -  272.] 

In  another  effort  to  utilize  the  geometric  foundations  of  chirality,  the  reflection  and  the 
transmission  characteristics  of  a  structurally  chiral  slab  had  been  investigated  in  the  previous  year. 
This  effort  has  been  further  extended  this  year.  By  stacking  up  unit  cells  made  of  (identical) 
structurally  chiral  slabs,  a  periodically  inhomogeneous  medium  was  constructed.  Each  structurally 
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chiral  slab  was  made  up  of  a  certain  number  of  identical  uniaxially  anisotropic  plates,  the 
consecutive  optic  axes  describing  either  a  right-  or  a  left-handed  spiral.  When  the  unit  cell  thickness 
is  very  small  compared  to  the  principal  wavelengths  in  the  uniaxial  plates,  the  periodically 
inhomogeneous  uniaxial  medium  was  shown  to  be  equivalent  to  a  homogeneous  biaxial  medium, 
the  two  optic  axes  of  the  equivalent  medium  being  dependent  on  the  handedness  of  the  periodic 
medium.  [V.V.  Varadan,  A.  Lakhtakia  &  V.K.  Varadan,  'Propagation  through  a  periodic  chiral  arrangement  of 
identical  uniaxial  dielectric  layers  and  its  effective  properties, '  accepted  for  publication  in  Optik ]  A  similar  result 
was  also  obtained  when  the  uniaxial  plates  were  replaced  by  biaxial  plates.  [A.  Lakhtakia,  V.K.  Varadan 
&  V.V.  Varadan,  ’Effective  properties  of  a  periodic  chiral  arrangement  of  identical  biaxially  dielectric  plates,'  accepted 
for  publication  in  J.  Mater. 

When  particles  of  considerable  concentration  are  dispersed  in  a  host  medium  to  form  a 
composite,  multiple  scattering  dominates  the  scattered  energy  when  waves  are  impinging  upon 
them.  While  this  effect  is  known  and  shown  for  various  dense  systems,  it  has  never  been  examined 
for  chiral  composites.  Therefore,  preliminary  investigations  were  undertaken  this  year.  The 
scattering  response  of  a  single  chiral  particle  is  cast  in  terms  of  a  T-matrix.  The  circular 
birefringence  of  chiral  media  does  not  generate  any  difficulties  in  the  present  multiple  scattering 
formulation,  if  and  only  if  the  composite  medium  is  assumed  to  be  achiral  (or  weakly  chiral). 
Although  a  rigorous  multiple  scattering  formalism  needs  to  be  introduced  when  the  effective 
medium  is  chiral,  at  least  in  the  low  frequency  regime,  using  Maxwell-Gamett  Approximation 
(MGA),  we  have  derived  the  dispersion  equation  for  an  effectively  chiral  composite  medium.  [V.V. 
Varadan,  Y.  Ma  &  V.K.  Varadan,  'Effects  of  chiral  microstructure  on  EM  wave  propagation  in  a  lossy  dielectric 
composite  material,'  accepted  for  publication  in  Radio  Science ]  The  low  frequency  expressions  of  the 
dispersion  equation  using  the  Bruggeman  approximation  were  also  obtained.  [V.V.  Varadan,  A. 
Lakhtakia,  Y.  Ma  &  V.K.  Varadan,  'Long  and  intermediate  wavelength  models  for  the  effective  properties  of  a  chiral 
composite  material,’  URSI  International  Symposium  on  Electromagnetic  Theory ,  Stockholm,  Sweden,  August  14  - 
17, 1989].  Finally,  computed  results  of  microwave  properties  of  chiral  composite  materials,  were 
obtained  by  the  numerical  solution  of  the  dispersion  equations.  [V JK.  Varadan,  Y.  Ma,  V.V.  Varadan  & 
A.  Lakhtakia,  'Electromagnetic  wave  propagation  in  chiral  composite  materials,'  First  Progress  in  Electromagnetics 
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Research  Symposium,  Boston,  Massachusetts,  USA,  July  25  -  26, 1989] 

Acoustic  Chirality 

The  premise  for  acoustic  chirality  in  elastic  solids  has  been  investigated  by  us  in  the  past,  and 
was  reported  in  the  Previous  Annual  Reports.  In  an  effort  to  synthesize  structurally  chiral  media, 
layered  composites  were  examined.  These  are  made  up  by  stacking  identical  uniaxial  plates,  their 
consecutive  symmetric  axes  describing  either  a  right-  or  a  left-handed  spiral.  The  field  equations 
governing  the  harmonic  motions  of  the  layered  composite  were  written  as  a  first  order  matrix 
ordinary  differential  equation,  where  the  field  variable  is  a  6-vector  consisting  of  the  displacement 
and  traction  components.  A  matrix  representation  method  was  used  for  solving  the  reflection  and 
transmission  characteristics  of  the  layered  composites;  it  is  noteworthy  that  this  method  is  not 
affected  by  the  number  of  layers.  Numerical  results  of  the  plane  wave  reflection  and  transmission 
characteristics  were  obtained  for  these  chiral  arrangements.  It  was  observed  that  the  co-polarized 
characteristics  are  unaffected  by  the  structural  chirality,  while  the  cross-polarized  reflected  and 
transmitted  electromagnetic  fields  are  greatly  influenced  by  it.  These  investigations  will  shortly  be 
submitted  for  publication. 

Experimental  research  on  the  use  of  chiral  composites  for  acoustic  attenuation  continued  during 
the  present  period.  The  composites  were  constructed  by  embedding  piezoelectric  springs  in  a 
polymeric  material.  Each  chiral  inclusion  was  also  backed  by  an  air-filled  cavity,  while  echo 
reduction  and  insertion  loss  were  measured  in  an  acoustic  pulse  tube.  The  results  of  our 
investigations  demonstrated  that  underwater  sound  attenuation  could  be  enhanced  by  embedding 
chiral-shaped  inclusions  in  polymeric  materials.  Further  research  is  in  progress. 

Effective  Properties  of  Ferromagnetic  Composites 

Ferrite  composite  materials  as  well  as  other  ferromagnetic,  carbonyl  iron  composites  have 
recently  received  wide  attention  due  to  their  various  applications  over  quite  a  broad  microwave 
band.  Electromagnetic  shielding  ability  and  wave  absorbing  efficiency  are  just  two  of  their  popular 
characteristics.  In  order  to  use  these  materials  in  different  environments  and  to  optimize  their 
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capabilities,  we  must  carefully  tailor  their  properties.  In  this  context,  we  have  developed  a  direct 
formalism  which  can  be  used  to  predict  the  effective  electromagnetic  properties  of  ferrite  composite 
materials  [Y.  Ma,  V.K.  Varadan  &  V.V.  Varadan,  'Prediction  of  electromagnetic  properties  of  ferrite  composites,' 
submitted  for  publication  in  Progress  in  Electromagnetic  Research  ]. 

In  addition,  for  many  types  of  ferrites,  it  is  hard  to  directly  measure  the  frequency  dependent 
complex  permittivity  and  permeability  due  to  their  high  electric  conductivities.  We  have  developed  a 
new  technique  to  predict  these  complex  properties  using  the  effective  property  measurements  done 
for  a  high  (or  low)  porosity  ferrite  composite,  as  well  as  our  mutiple  scattering  formalism.  [Y.  Ma, 
V.K.  Varadan  &  D.K.  Ghodgaonkar,  'Inverse  problems  for  ferromagnetic  composites,'  91st  Annual  Meeting  of 
American  Ceramic  Society,  Indianapolis,  Indiana,  April  23-27, 1989]. 

Modified  Flux  Models  and  Their  Applications  to  Coating  Design 

In  two-  and  four-  flux  models  of  radiative  transfer  theory,  the  scattering  coefficients  or 
efficiencies  of  non-emitting  media  are  commonly  computed  using  the  single-scattering  albedo,  while 
the  interactions  among  particles  are  neglected.  The  flux  models  are  modified  by  considering  multiple 
scattering.  One  application  of  the  modified  flux  model  is  to  predict  the  opacity  of  high  pigment 
volume  concentration  paint  coatings.  Therefore,  we  introduced  the  multiple  scattering  formalism 
and  have  shown  from  the  resulting  equations  how  the  multiple  scattering  albedo  can  be  calculated 
from  the  effective  propagation  constant  K.  The  multiple  scattering  albedo  so  obtained  can  be  used  in 
the  flux  model  to  compute  the  diffuse  reflectance  of  the  corresponding  paint  film. 

Because  the  refractive  index  of  pigment  particles  used  in  the  measurements  is  direction- 
dependent,  the  modeling  of  an  anisotropic  scatterer  wass  introduced.  For  an  unpolarized  light 
source,  illuminated  anisotropic  particles  of  random  orientation  can  be  thought  of  as  having  an 
average  refractive  index  over  the  optical  axes  within  the  particles  in  most  of  the  visible  spectrum. 
This  information  is  used  to  compute  scattering  cross  section  or  the  T-matrix. 

We  also  reviewed  the  existing  formalisms  of  the  flux  models.  The  coupled  differential 
equations  for  the  diffuse  reflectance  and  transmittance  are  solvable  in  general.  However,  the 
coefficients  in  the  equations  were  modified  by  incorporating  multiple  scattering.  The  boundary 
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conditions  for  the  real  paint  films  cannot  be  ignored  and  we  considered  two  different  diffuse 
substrates  in  order  to  justify  the  opacity  through  the  contrast  ratio.  Though  in  most  cases  theoretical 
results  compare  favorably  well  with  those  measured,  the  discrepancies  in  some  case  due  to  either 
the  anisotropic  modeling  or  inaccuracy  in  the  values  used  for  the  substrate  are  discussed  [Y.  Ma, 
V.V.  Varadan  &  V.K.  Varadan,  'Modified  flux  models  considering  multiple  scattering  -  application  to  opacity  of 
coatings,'  submitted.  Applied  Optics ]. 

Enhanced  Absorption  due  to  Dependent  Scattering 

In  a  recent  series  of  studies  on  light  scattering  from  pigmented  surface  coatings  or  by  radiative 
heat  transfer  in  powder  insulators,  it  has  been  shown  that  for  densely  populated  heterogeneous 
media  dependent  scattering  as  well  as  dependent  absorption  have  to  be  considered.  As  a  result,  for 
paint  opacity,  the  dependent  scattering-absorption  can  greatly  affect  the  diffuse  reflectance  as  it  does 
the  overall  thermal  radiation  resistance  of  the  powder  insulation  composites.  Analogous  to  the 
aforementioned  applications,  chemical  and  nuclear  reactors,  fuel  combustors,  cryogenic  insulation, 
microwave  and  laser  coatings,  artificial  obscuration  materials  and  many  other  commercial  and 
military  systems  usually  involve  a  high  concentration  of  particles.  In  order  to  analyze  the  energy 
transport  characteristics,  the  dependent  scattering  and  absorption  properties  cannot  be  ignored. 

The  role  of  multiple  scattering  on  dependent  scattering  as  well  as  on  dependent  absorption  was 
investigated  for  heterogeneous  media  containing  a  high  concentration  of  particles.  The  decrease  of 
scattering  and  increase  of  absorption  for  lossy  (with  intrinsic  absorption)  particles  in  a  lossless 
matrix  was  quantitatively  described  by  the  use  of  a  Rayleigh  region  solution  derived  from  a  multiple 
scattering  formalism.  For  smaller  wavelengths,  dependent  scattering  and  absorption  was  obtained 
through  nymerically  solving  the  resulting  dispersion  equation.  Numerical  results  were  computed  for 
lossless  particles  in  a  lossy  matrix  which  models  an  optical  coating  system,  [(i)  Y.  Ma,  VX  Varadan  & 
V.V.  Varadan,  'Enhanced  absorption  due  to  dependent  scattering,'  accepted  for  publication  in  ASME  Journal  of  Heat 
Transfer;  (ii)  Y.  Ma,  VX.  Varadan  &  V.V.  Varadan,  'Enhanced  absorption  due  to  dependent  scattering,'  1989 
CRDEC  Scientific  Conference  on  Obscuration  and  Aerosol  Research ,  Aberdeen  Proving  Ground,  Maryland,  June 
26-30, 1989] 
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Modeling  of  Dielectric  Properties  of  Porous  Ceramic  Materials 

In  our  experimental  studies  on  the  microwave  sintering  of  ceramics,  it  has  been  observed  that 
the  (initial)  porosity  of  the  green  sample  predetermines  the  subsequent  sintering  rate  [V.K.  Varadan,  Y. 
Ma,  A.  Lakhtakia  and  V.V.  Varadan,  ‘Microwave  sintering  of  ceramics,'  in  Microwave  Processing  of  Materials  (Eds.: 
W.H.  Sutton,  M.H.  Brooks  &  IJ.  Chabinski),  Pittsburgh:  MRS,  l988].The  change  in  porosity  during 
sintering  is  a  complex  thermodynamic  process,  rooted  in  the  absorbed  microwave  energy;  the 
instantaneous  absorption  depends  on  the  temperature-dependent  dielectric  loss  of  the  porous 
sample. 

We  have  proposed  a  model  to  quantitatively  describe  the  change  in  the  dielectric  properties 
using  the  effective  two-body  inter-particle  potentials.  These  potentials  depend  on  the  porosity  as 
well  as  on  the  temperature.  At  the  critical  porosity  (or  ,the  critical  temperature),  a  phase  switch  can 
thus  be  predicted  [V.K.  Varadan,  Y.  Ma  and  V.V.  Varadan,  Modeling  of  dielectric  properties  of  porous  ceramics 
during  microwave  sintering,'  91st  Annual  Meeting  of  American  Ceramic  Society,  Indianapolis,  Indiana,  April 
23-27,  1989].  In  addition,  using  analytical  expressions  [Y.  Ma,  V.K.  Varadan  &  V.V.  Varadan,  'Effective 
properties  of  microwave  composites,' /.  Wave-Mater .  Interact.  3, 1988, 243  -  248]  we  have  predicted  dielectric 
properties  of  some  sintered  ceramic  materials.  Hie  computed  results  compare  quite  well  with  those 
measured  which  were  provided  by  Georges  Roussy  from  Universite'  de  Nancy  1,  France. 
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The  subject  of  power  sonic  and  ultrasonic  transducers  has  been  rapidly 
expanding  on  an  international  scale  over  the  last  several  years.  This  book  is 
a  collection  of  papers  resulting  from  the  Proceedings  of  the  International 
Workshop  on  the  Design  of  Power  Sonic  and  Ultrasonic  Transducers,  as 
held  in  Lille,  France,  on  26  and  27  May  1 987.  A  presentation  of  some  of  the 
latest  transducer  design  techniques,  applications,  and  available  materials 
are  given.  The  complete  text  includes  222  figures  and  is  written  ih  English. 

The  contents  of  the  book  begin  with  an  introductory  lecture  by  B. 
Tocquet  and  follows  with  the  13  papers  written  by  each  of  the  invited  Work¬ 
shop  contributors.  Each  chapter  is  an  individual  paper,  and,  as  such,  they 
can  stand  on  their  own. 

The  first  chapter  is  entitled  “Power  Limitations  of  Piezoelectric 
Length  Expander  Transducers,”  by  Oscar  B.  Wilson  of  the  U.S.  Naval 
Postgraduate  School  in  Monterey,  California.  This  paper  discusses  some  of 
the  internal  and  external  power  limiting  factors  of  electroacoustic  trans¬ 
ducers.  Effective  use  of  an  equivalent  network  is  provided  with  the  example 
discussion  of  a  longitudinal  vibrator. 

L.  Eyraud  of  the  Laboratoire  de  Genie  Electrique  et  Ferroelectricite  in 
Villeurbanne,  France  authors  the  second  paper  “The  Material  for  Piezo¬ 
electric  Power  Transducers”  and  coauthors  the  third  paper  “Characteriza¬ 
tion  of  Piezoelectric  Ceramics  for  High  Power  Transducers"  with  P.  Gon- 
nard  and  P.  Champ  of  the  same  address.  Chapter  2  discusses  the 
mechanisms  of  piezoceramic  stability  as  a  function  of  aging  and  polariza¬ 
tion  through  a  discussion  of  chemical  compositions,  while  Chap.  3  charac¬ 
terizes  various  ferroelectric  parameters.  An  analysis  of  piezoceramics  under 
cw  and  pulsed  conditions  is  presented  with  the  limiting  factors  of  each 
detailed. 

The  fourth  chapter  is  titled  "Highly  Magnetostrictive  Rare  Earth 
Compounds  for  High  Power  Acoustic  Projectors,”  by  Arthur  E.  Clark  from 
the  Naval  Surface  Weapons  Center  (NSWC)  in  Silver  Spring.  Maryland. 
This  paper  is  a  reproduction  of  a  chapter  (7)  written  in  1980  in  the  book 
Ferromagnetic  Materials,  Volume  1,  published  by  North-Holland,  Amster¬ 
dam.  It  presents  a  background  and  detailed  discussion  of  magnetostrictive 
lanthanide  compounds.  Although  the  discussion  appears  complete,  much 
of  the  literature  is  dated  and  does  not  include  the  excellent  recent  work  done 
by  Clark  and  his  colleagues  at  NSWC. 

D.  Boucher  of  Groupe  d'Etude  et  de  Recherche  en  Detection  Sous- 
Marine  in  Toulon,  France  penned  Chap.  5  “Trends  and  Problems  in  Low 
Frequency  Sonar  Projectors  Design."  A  discussion  on  various  classical  so¬ 
nar  transducer  designs  is  presented  along  with  an  effective  briefing  on  array 
interaction  effects.  Unfortunately  several  of  his  comments  on  specific  trans¬ 
ducer  designs  are  quite  misleading.  For  example,  when  he  states  “...the 
work  in  this  Aeld  (Class  IV  flextensional  design)  has  been  pursued  as  far  as 
possible,”  he  ignores  the  recent  resurgence  of  the  Class  IV  flextensional 
design  utilizing  lanthanide  drivers,  high  strength  epoxy  shell  composites, 
and  new  multidimensional  array  designs.  On  the  other  hand,  his  summary  is 


an  accurate  portrayal  of  the  present  problems  being  encountered  in  low- 
frequency  sonar  transducers. 

“Frequency,  Power  and  Depth  Performance  of  Class  IV  Flextensional 
T ransducers"  is  the  title  of  Chap.  6  as  cowritten  by  J.  Oswin  and  1.  Dunn  of 
British  Aerospace,  Underwater  Research  and  Engineering  Unit  in  Wey¬ 
mouth,  United  Kingdom.  This  paper  gives  an  excellent  overview  of  the 
Class  IV  flextensional  transducers  from  the  original  patents  to  present  Dis¬ 
cussions  on  depth,  power,  and  frequency  performance  are  given  along  with 
simple  design  guidelines.  A  correction  on  the  dates  of  the  two  Toulis  patents 
( references  6. 1  and  6. 1 1 )  must  be  clarified  though.  The  authors  incorrectly 
listed  the  patents  as  being  granted  in  1963  but  in  actuality  the  patents  were 
first  filed  in  1963  and  granted  in  1966.  A  minor  point  but  one  that  should  be 
quantified  to  establish  an  accurate  historical  setting. 

Chapter  7  is  titled  “Opportunities  and  Challenges  in  the  Use  of  Ter- 
fenol  for  Sonar  Transducers,”  by  J.  M.  Sewell  and  P.  M.  Kuhn  of  Martin 
Marietta  Aerospace  in  Baltimore,  Maryland.  This  paper  begins  with  a  brief 
comparison  between  piezoceramic  and  terfenol,  and  is  then  followed  with  a 
discussion  of  terfenol  transducer  performance  characteristics.  Even  though 
the  authors  found  several  positive  points  in  the  designs  presented,  the  tech¬ 
nology  shown  here  is  antiquated.  Their  methods  for  applying  both  the  dc 
(magnetic)  and  the  mechanical  bias  are  ineffective  and  will  result  in  low- 
efficiency  devices.  The  educated  reader  will  readily  notice  that  the  latest 
reference  given  is  1977,  despite  the  literature  and  patents  published  by  the 
Naval  Underwater  Systems  Center  (NUSC),  Raytheon  Company,  and  Im¬ 
age  Acoustics,  among  others,  that  have  advanced  this  technology  much 
further  along. 

“Application  of  the  Finite  Element  Method  to  the  Design  of  Power 
Piezoelectric  Sonar  Transducers”  and  “Determination  of  the  Power  limits 
of  a  High  Frequency  Transducer  Using  the  Finite  Element  Method”  are  the 
titles  to  Chaps.  8  and  9,  respectively.  The  first  one,  authored  by  B.  Harmon¬ 
ic  (one  of  the  two  book  editors )  of  the  Institut  Superieur  d'Electronique  du 
Nord,  Laboratoire  d’Acoustique,  in  Lille,  France,  presents  FEM  analysis 
(with  ATI  LA  code)  on  several  sonar  transducer  designs.  Chapter  9,  co¬ 
written  by  W.  Steichen,  G.  Vanderborck,  and  Y.  Lagicr  of  Thomson-Sintra 
in  Valbonne,  France,  explains  a  multiprocedure  FEM  analysis  that  includes 
acoustic  radiation  and  thermal  and  dielectric  losses.  Although  this  method 
appears  rather  tedious,  the  experimental  and  computed  examples  are  in 
general  agreement. 

The  next  two  chapters  present  discussions  on  novel  new  ultrasonic 
transducer  designs  for  industrial  applications.  Chapter  10  is  entitled  “High 
Power  Ultrasonic  Transducers  for  Use  in  Gases  and  Interphases,”  by  J.  A. 
Gallego-Juarez  of  the  Instituto  de  Acustica  in  Madrid,  Spain,  while  Chap. 

1 1  is  entitled  “Design  of  High  Power  Ultrasonic  Transducers  for  Use  in 
Macrosonics,”  by  P.  Tierce  of  Sinaptec  Sari  and  J.  N.  Decarpigny  (the 
second  editor  of  the  book)  of  the  Institut  Superieur  d'Electronique  du 
Nord,  Laboratoire  d’Acoustique,  both  in  Lille,  France. 

The  functions  and  architecture  of  a  sonar  emitting  system  are  present¬ 
ed  in  Chap.  12  under  the  heading  “Power  Electronic  Devices  for  Sonar 
Systems,"  by  C.  Pohlenz  of  Groupe  d'Etude  et  de  Recherche  en  Detection 
Sous-Marine  (GERDSM)  in  Toulon,  France.  This  paper  offers  several 
methods  for  solving  difficult  impedance  matching  problems  in  a  sonar  sys¬ 
tem.  Disappointingly,  it  is  the  only  paper  that  did  not  offer  any  references 
for  further  study.  The  final  paper  in  the  book.  Chap.  13,  is  "The  Electronic 
Driving  Sources  for  Ultrasonic  Machining,”  by  B.  Thirion  of  Slice  Sari  in 
Lille,  France. 
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In  summary,  although  I  found  several  papers  reporting  on  “past  glor¬ 
ies,”  I  feel  that  most  transducer  designers  and  material  scientists  will  benefit 
from  this  book.  It  offers  a  rare  glimpseof  the  international  transducer  world 
with  a  wide  variety  of  useful  global  references.  It  is  not  limited  to  the  under¬ 
water  sonar  designer  since  it  additionally  offers  innovative  ultrasonic  de¬ 
vices  that  may  be  found  in  use  in  the  private  sector. 

THOMAS  R.  HOWARTH 

Department  of  Engineering  Science  and  Mechanics 

Research  Center  for  the  Engineering  of  Electronic  and  Acoustic  Materials 

The  Pennsylvania  State  University 

University  Park,  Pennsylvania  16802 


The  Journal  of  the  Acoustical  Society  of  America 
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TECHNICAL  NOTES  AND  RESEARCH  BRIEFS 

Advanced-degree  dissertations  in  acoustics 


Experimental  studies  using  chiral  composites  as  acoustic  energy 

attenuators  [43.30.Ky] — Thomas  R.  Howarth,  The  Research  Center  for 
Engineering  of Electronic  and  Acoustic  Materials.  Department  of  Engineer¬ 
ing  Science  and  Mechanics.  The  Pennsylvania  State  University.  University 
Park.  PA  16802.  May  1988 (M.S).  The  area  of  underwater  acoustic  attenu¬ 
ation  has  been  a  subject  of  much  research  since  the  1940s.  The  Germans  first 
investigated  the  use  of  air-rubber  baffles  to  camouflage  their  submarines  for 
active  sonar.  In  their  study  they  also  developed  a  water-filled  acoustic  pulse 
tube  device  for  determining  air-rubber  performance  characteristics  under 
various  hydrostatic  pressures  and  temperatures.  In  the  late  1970s  and  into 
the  1980s  other  types  of  designs  have  been  sought.  It  was  found  by  using  the 
concept  of  a  T  matrix  in  conjunction  with  multiple  scattering  theory  that 
the  use  of  piezoelectric  chiral-shaped  elements  could  replace  the  air  pockets 
of  the  previous  designs  and  offer  new  mechanisms  for  acoustic  energy  at¬ 
tenuation.  This  thesis  investigates  experimental  studies  using  piezoelectric 
chiral  compositics  for  acoustic  attenuation.  A  discussion  of  the  theoretical 
concepts  is  given  followed  by  the  techniques  used  for  fabrication  of  test 
samples.  Background  and  operation  of  the  pulse  tube  measurement  system 
are  provided  followed  by  a  discussion  of  measured  results.  Recommenda¬ 
tions  for  future  investigations  conclude  the  presentation. 


Thesis  Advisor.  Vijay  K.  Varadan. 


Incommensurate  Numbers,  Continued  Fractions, 
and  Fractal  Immittances 


A.  Lakhtakia*,  R.  Messier*1**,  V.  V.  Varadan*,  and  V.  K.  Varadan* 

The  Pennsylvania  State  University,  University  Park,  PA  16802  (USA) 

Z.  Naturforsch.  43a,  943-955  (1988);  received  June  10,  1988 

Continued  fractions  have  a  rich  tradition  in  the  theory  of  numbers;  c.g.,  non-terminating  con¬ 
tinued  fractions  represent  irrational  numbers.  It  will  be  shown  that  a  class  of  continued  fractions 
possess  the  property  of  self-referential  decomposition,  and  their  interpretation  in  the  form  of  non- 
terminating  ladder  circuits  gives  rise  to  fractal  immittances  with  potential  analogies  to  rough 
surfaces,  thin  cermet  films,  as  well  as  to  the  internal  void  network  structure  of  thick  films. 


Introduction 

The  motivation  for  this  work  comes  from  number 
theory:  the  representation  of  irrational  numbers  in  the 
Stieltjes  continued  fraction  form,  and  which  has  al¬ 
ready  inspired  the  development  of  ladder  circuits  [1], 
Continued  fractions  have  recently  been  used  in  under¬ 
standing  fractal  quantization  of  particles  in  one¬ 
dimensional  potentials  with  incommensurate  periods 
[2],  as  well  as  in  two-dimensional  electron  gases  [3]:  a 
particularly  appealing  and  simple  interpretation  of 
the  relevant  Hamiltonian  has  been  given  by  Chao  [4], 
Continued  fractions  have  also  been  applied  in  examin¬ 
ing  the  frustrated  instabilities  of  active  optical  resona¬ 
tors  [5]. 

Although  continued  fractions  have  seen  some  use  in 
the  characterisation  of  the  rough  surfaces  of  real  mate¬ 
rials  [6,7],  they  have  not  been  applied  yet  to  the 
inhomogeneous  internal  structure  of  materials.  It  will 
be  shown  that  specific  examples  of  the  resulting  ladder 
circuits  have  fractal  immittances  with  potential  ana¬ 
logies  to  both  cermet  thin  films  and  films  with  internal 
void  networks;  and,  hence,  to  the  related  film  proper¬ 
ties. 


*  Department  of  Engineering  Science  and  Mechanics. 

**  Materials  Research  Laboratory.  , 

Reprint  requests  to  Prof.  R.  Messier,  Engineering  Science 
and  Mechanics,  265  Materials  Research  Labroalory,  The 
Pennsylvania  State  University,  University  Park,  PA  16802, 
USA 


Continued  Fractions  and  Quadratic  Irrational  Numbers 

Any  rational  or  irrational  number  can  be  written 
down  in  the  continued  fraction  form  as  [8] 

1 

<a0,a1,a2,a3,...>=  a0  + - j - 

a,  H - - - 

a2  +  — — — 

«»  +  •••  .  0) 

with  a0,  a,,  etc.  being  positive  integers  [but  see  Ap¬ 
pendix],  The  continued  fraction  is  finite  in  size  when 
it  represents  any  rational  number  p/q,  with  a„  being 
the  integral  part  of  the  ratio  p/q;  it  turns  out  that  every 
rational  number  has  exactly  two  such  representations 
[9].  On  the  other  hand,  irrational  numbers  have  con¬ 
tinued  fraction  representations  which  are  infinite  in 
size.  Of  these  irrational  numbers,  there  is  a  class  of 
quadratic  irrational  numbers  which  are  solutions  of  a 
quadratic  equation.  A  theorem  due  to  Lagrange  [10] 
states  that  the  continued  fraction  expansion  of  any 
quadratic  irrational  is  periodic  after  a  certain  state, 
eg., 

v/15  =  <3,1.6, 1,6, 1,6,...), 
n/31  =  <5, 1, 1,  3, 5,  3, 1, 1, 10, 1, 1,  3, 5,  3, 1, 1, 
10,1,1,3,  5,  3, 1,1, 10,... >, 

{24  -  yi5}/l  7  =  <1, 5, 2, 3, 2, 3, 2, 3, . .  .> . 

It  is  the  representation  of  quadratic  irrational  num¬ 
bers  which  is  of  particular  interest  in  the  present  con¬ 
text.  The  simplest  example  is  the  number 

Q(m)  =  <0,p,m,p,...>;  p>0,  (2) 
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ABSTRACT 

The  refraction  of  near  fields  by  a  planar  achirallchiral  interface  has  been  examined.  A  realistic  source,  a 
constant  current  line  source,  has  been  used  here,  as  apart  from  the  usual  analyses  involving  incident 
planewaves  only.  Provided  a  planewave  spectral  decomposition  of  the  incident  field  is  possible,  this 
procedure  can  be  extended  to  include  sources  of  other  configurations  and  polarizations  as  well.  Maps  of  the 
refracted  field  are  drawn  to  elucidate  the  effect  of  both  the  near-zone  irradiation  as  well  as  of  the  handedness  of 
the  chiral  medium.  Such  an  analysis  would  be  of  use  to  various  researchers  in  the  area  so  that  the  effect  of 
the  near fields  cf  radiating  sources  may  not  be  ignored. 

1.  INTRODUCTION 

Ever  since  the  discovery  of  optical  activity  early  in  the  last  century,  there  has  been  a  great  interest  in 
measuring  the  circular  dichroic  (CD)  and  the  optical  rotatory  dispersion  (ORD)  spectra  of  molecular 
aggregations  [1].  Such  measurements  are  routinely  made,  at  frequencies  down  into  the  infra-red  regime 
nowadays.  The  measurement  procedures  usually  are  variants  of  a  simple  technique:  A  planar  slab  of  the 
optically  active,  or  chiral,  material  is  irradiated  by  a  source,  the  polarization  state  of  the  irradiating  field  being 
known.  On  the  other  side  of  the  slab,  the  rotation  of  the  plane  of  polarization  of  the  transmitted  field  with 
respect  to  that  of  the  incident  field  is  measured,  and  the  CD  and  the  ORD  of  the  material  determined.  A  good 
description  of  such  an  experiment  conducted  at  frequencies  around  10  GHz  is  available  in  [2],  where  Tinoco 
and  Freeman  describe  the  procedure  for  measuring  the  optical  rotatory  activity  of  a  collection  of  oriented 
copper  helices,  each  approximately  0.S  cm  in  diam  and  about  1  cm  long. 

However,  this  and  other  such  measurement  techniques  usually  ignore  the  proximity  of  the  chiral  slab  to 
the  source.  If  the  source  is  far  away  from  the  exposed  face  of  the  slab,  the  illuminating  field  may  be 
conveniently  taken  to  be  a  planewave.  But  when  that  is  not  so,  the  slab  lies  in  the  near-zone  of  the  source, 
where  the  irradiating  field  has  reactive  components  of  large  magnitudes.  In  such  circumstances,  the  planewave 
approximation  of  the  source  field  can  be  quite  gross.  In  this  communication,  the  nature  of  the  fields  excited 
in  a  chiral  half-space  are  examined.  The  source  is  taken  to  be  a  constant  current  line  source,  the  field  radiated 
by  which  source  is  decomposed  into  an  infinite  set  of  planewaves,  some  of  which  are  propagating  and  the 
remaining  ones  are  evanescent  It  is  to  be  emphasized  here  that  any  other  source  can  be  accommodated  in  this 
theoretical  procedure,  provided  its  field  can  be  expressed  by  a  planewave  spectral  (PWS)  representation  [3].  By 
mapping  the  field  refracted  into  the  chiral  half-space,  some  understanding  of  the  effect  of  near-field  irradiation 
of  achiral/chiral  interfaces  is  obtained. 

2.  REFRACTION  OF  A  TE-POLARIZED  PLANEWAVE  BY  A  PLANE 
ACHIRAL/CHIRAL  INTERFACE 

Let  the  space  z  £  0  be  filled  by  a  non-dissipative  achiral  medium  (here  taken  to  be  free-space)  in  which 
the  constitutive  relations 

D  =  eE  ,  B  =  hoH  (la,b) 


hold,  and  from  which  a  TE-poiarized  planewave 
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Ej(k)  =  j  exp  □  {kx  +  {3o(k)z}]  ,  H.(k)  =  (l/jcofxo)  V  x  E.(tc)  ,  (2a,b) 

with 

pi>(K)=V0-icJ]'n  ,  ko  =  <0ltiA]'B  ,  <3*.b) 

is  incident  on  the  interface  z  =  0,  the  hannonic  time-dependence  exp[-jcot]  being  suppressed  here  and  hereafter, 
and  the  unit  vectors  i,  j,  k  having  their  usual  meanings  in  a  rectangular  co-ordinate  system. 

The  region  above  the  z  =  0  interface,  z  £  0,  is  occupied  by  a  chiral  medium  in  which  the  constitutive 


equations  are  given  as  [4,5] 

D  =  e{E  +  aVxE}  ,  B  =  p{H  +  aVxH}  (4a,b) 

and  in  which  it  would  be  proper,  because  of  Snell's  law,  to  express  the  transmitted  field  as  [6]: 

Et(K)=TL(K)QL(K)  +  aRTR(K)QR(K)  ;  (5a) 

H(k)  =  Tl(k)  Ql(k)  +  Tr(k)  Qr(k)  (5b) 

in  which  T^rOc)  are  to  be  determined  by  the  boundary  conditions  on  the  interface  z  =  0,  and  the  left-  and 
right-circulariy  polarized  (LCP  and  RCP)  wavefunctions,  respectively,  are  given  as 

Ql(k)  =  (l/kL)  {-pL(K)i  -  jkLj  +  Kk }  exptJ  {  kx  +  Pl(k>}]  ,  (6a) 

and 

Qr  (k)  =  (l/kR){-PR  (K)i  -jkRj+  Kk}  exp[j  {  kx  +  pR  (k)z}  ]  .  (6b) 

In  Eqs.  (5)  and  (6)  the  following  definitions  hold: 

aL=-j(e/p]1/2  ,  aR  =  -j[n/e],/2  ,  (7a,b) 

kL  =  k[l-kaf1  ,  kR  =  k[l  +ka]_1  ,  (7c,d) 

pL(K)  =  +  [kL-K2]1/2  .  PR(K)  =  +  [kR-K2]I/2  ,  (7e,f) 

k  =  o[n£]1/2  .  (7g) 


The  field  reflected  back  into  the  achiral  medium  also  has  LCP  and  RCP  components,  and  it  can  be 
expressed  in  the  form  [6]: 

Er(K)  =  Rr(k) Pr(k)  +  Rl(k) PL(K)  ,  Hr(K)  =  {\rm0)  V X  E(k)  ,  (8a, b) 

where  R^rW  are  to  be  determined  and  the  wavefunctions, 

PR(K)  =  (lyko)  {Po(K)i+jkJ  +  Kk}cxp[j{Kx-Po(K)z)]  ,  (9a) 

and 

-PL (k)  =  (iyko)  (Po(K)i -jkj  +  Kk}  exp(j  {kx  - P0(k)z} ]  .  (9b) 

By  invoking  the  usual  boundary  conditions  across  the  z  s  0  interface,  i.e., 

kx[E.(K)  +  Er(K)-Et(K)]  =  0  ,  kx[Hi(K)  +  H/K)-Ht(K)]  =  0  ,  (10a, b) 

and  using  Eqs.  (1),  (5)  and  (8),  it  can  be  shown  that 

Tl(k)  =  2j  [S  +  aR  Q(k)]  /  A(k)  ,  (1  la) 
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TR(K)  =  2j[aLS-P(K)]/A(K)  ,  (lib) 

RL(K)  =  -(l/4)[TL(K){P(K)+l}{aLS  +  l}+TR(K){Q(K)-l){S+aR}]  ,  (11c) 
RR(K)=(l/4)[TL(K){P(K)-l}{aLS-l}+TR(K){Q(K)+l){S-aR}]  ,  (lid) 
A(k)  =  [S  +  ^ Q(k)]  [1  -  ^ P(k)  S]  -  \  S  - P(k)]  [aR  +  Q(k)  S]  ,  (lie) 

aid 

1/2 

P(K)  =  (ko/kL)[pL(K)/|5o(K)]  ,  Q(K)  =  (k<ykR)[PR(K)/po(K)]  ,  S=-j[R0/e/  . 

(llf.g.h) 

It  should  be  noted  that  if  a  =  0,  i.e.,  the  medium  above  the  z  s  0  interface  also  becomes  achiral,  then 


kR  =  kL=k  ,  pR(K)  =  pL(K)=p(K)  =  +  [k2-K2]1/2  ,  (12a,b) 

P(K)=Q(K)  =  (k{/k)[p(KVpo(K)]  ,  (12C) 

TL(K)  =  -aRTR(K)  =  jaR[aR  +  Q(K)Sr1  ,  (12d) 

and 

Rr(k)  =  -  Rl(k)  =  (l/2)j  -  Q(K)  S]  [a*  +  Q(k)  S]*1  ;  (12e) 

consequently,  the  familiar  expressions  available  in  any  standard  textbook  [e.g.,  7]  on  EM  theory  emerge: 

Et(K)=j2aR[aR  +  Q(K)Sf1exp[j  {kx  +  P(k)z}]  ,  (13a) 

and 

E/K)  =  -j  [aj^-QOO  S]  [aR+  Q(k)  S]"1  exp[j  {kx-  Po(k)z}]  .  (13b) 


3 .  REFRACTION  OF  THE  FIELD  OF  A  LINE  SOURCE 

Let  now  an  isotropically  radiating,  y-directed  current  line  source  be  located  at  Tp  =  -kdp  in  the  achiral 
half-space,  and  its  radiated  field  be  set  as  [8]: 

EjsjjiHQ  (klr — rpl)  .  (14) 

A  planewave  spectral  (PWS)  decomposition  [3]  of  this  field  is  possible,  and  is  given  by 

•• 

Ej  =  j J[dK/po(K)]exp(jpo(K)dp]exp[j{Kx  +  po(K)z)]  ,  z>-dp  ,  (15a) 

and 


Ej  =  j  J  [dtc/Po(K)] exp[-jPo(K)dp] exp[j  {kx -  P0(k)z)]  ,  z<-dp  .  (15b) 

In  view  of  the  PWS  decomposition,  Eq.  (15a)  and  the  development  of  the  preceding  section,  it  is 
possible  to  write  down  the  total  field  existing  in  the  chiral  half-space  due  to  the  line  source  as 
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V  J [dK/po(K)l cxp(jp0(K)dp] •  (Tl(k) Ql(k)  +  a* Tr(k) Qr(k)}  ,  z>0  .  (16) 

Likewise,  provided  it  is  assumed  that  the  field  reflected  from  the  achiral/chiral  interface  does  not  influence  the 
line  source  current,  the  total  field  existing  in  the  space  z  <  -dp  can  be  set  down  as 

M 

ER  =  J  [die/po(K)]  exp(jpo(K)dp]  •  (R^k)  Pr(k)  +  Rl(k)  Pl(k)}  + 

M 

+j  J[dJC/po(K)]exp[-jpo(K)dp]exp|j  {kx-Po(k)z}]  ,  z<-dp  .  (17) 


4.  DISCUSSION 

It  is  the  function  Ep  which  has  to  be  mapped  in  the  xz  plane  in  order  to  understand  the  fields  generated 
in  the  chiral  half-space.  It  is  obvious  that  the  K-domain  poles  can  be  extracted  simply  by  solving  the 
equation  A(k)  *  0.  However,  such  poles  are  of  little  interest  in  the  present  analysis.  Instead,  the  only  recourse 
is  to  actually  map  Ep  over  an  xz  grid  for  a  given  set  of  parameters.  It  would  certainly  be  useful  then  to 
exploit  the  x-  and  the  K-symmetries  of  the  various  parts  of  the  integrands  in  Eq.  (16),  and  it  turns  out  that 

eo 

ET=2  J  [dK/kLPo(K)]  exp[jPo(K)dp]  exp(jPL(K)z]  Tl(k)  {-iPL(K)cosKx-jjkLcosKx+kjtcsinKx )  + 

0 


+  2aR  J  (dK/kRpo(K)]  exp(jPo(K)dp]  exp|jpR(K)z]  Tr(k)  { -iPR(K)cosKx+jjkRcosKx+kjKsin  kx  }  , 

0  z>0 ,  (18) 

whose  cartesian  components  have  either  even  or  odd  symmetries  with  respect  to  the  x  coordinate. 

A  computer  program  to  calculate  Ep  was  implemented  on  a  DEC  VAX  1 1/730  minicomputer  via  Eq. 
(18).  The  infinite  K-integral  in  Eq.  (18)  was  truncated  to  hold  over  the  range  0  £  K/k0  5  20.0,  it  being 
observed  that  this  truncation  did  not  give  rise  to  any  errors  for  the  selected  values  e/Eg  £  5.0,  p/Po  =  1.0,  otk 
£  0.1  and  kgdp  =  2.0.  The  magnitudes  of  the  x-,  y-,  and  z-  directed  components  of  the  LCP  and  the  RCP 
parts  of  Ep,  as  well  as  of  the  total  Ep,  are  illustrated  in  Figs.  1  -  4  for  several  cases  over  the  rectangular 
region  0  £  kx/dp,  kz/dp  £  10.0. 

In  Fig.  1,  Ep  =  Ej  since  e/e0  =  1.0  and  ok  »  0.0.  This  illustration,  therefore,  simply  shows  the  field 
radiated  by  the  line  source  in  the  specified  xz  domain  when  all  space  is  covered  by  the  non-dissipative  achiral 
medium.  As  was  expected  and  can  be  observed  from  this  figure,  the  magnitudes  of  the  RCP  and  the  LCP 
parts  of  Ep  are  identical,  but  Ep  is  only  y-directed.  In  Fig.  2,  the  ratio  e/e0  was  increased  to  5.0  with  ak 
still  equal  to  zero.  Again,  Ep  turns  out  to  have  only  a  y-directed  component,  the  x-  and  the  z-  directed 
components  of  its  LCP  and  RCP  parts  having  cancelled  themselves  out  These  two  figures  are  then  simply 
symbolic  of  what  can  be  expected  when  the  refracting  half-space  is  also  achiral  [8]. 

This  cancellation  of  the  x-  and  the  z-  directed  components  of  the  LCP  and  the  RCP  parts  of  Ep  does  not 
occur  in  Fig.  3  where  e/Cq  ■  5.0  and  ak  =  0.01.  The  magnitudes  of  the  LCP  and  the  RCP  parts  of  Ep  are 
still  equal  by  virtue  of  the  fact  that  all  possible  values  of  k  ( positive  as  well  as  negative)  are  included  in  Eq. 
(18).  However,  the  sum  of  these  parts  still  shows  the  presence  of  the  x-  and  the  z-  directed  components, 
which  had  cancelled  out  in  Figs.  1  and  2  where  ak  =  0.0.  The  effect  of  handedness  in  the  half-space  z  2  0  is 
betrayed,  thus,  by  the  presence  of  a  transmitted  electric  field  Ep  which  is  not  TE  polarised  like  the  source 
field  E|.  This  tendency  is  even  more  marked  in  Fig.  4  where  e/e0  =  5.0  and  ak  =  0.1.  Furthermore,  in 
Fig.  4,  Ep  records  well-defined  extrema,  which  are  not  as  prominent  in  Fig.  3.  It  becomes  possible  to  state. 
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Figure  1  Magnitudes  of  the  various  components  of  Ej  of  Eq.  (18)  when  Ej  is  given  by  Eq.  (14).  From  left 
to  right  in  each  row:  magnitudes  of  the  x-,  y-  and  z-  directed  components.  From  top  to  bottom: 
LCP  part  of  Ej,  RCP  part  of  Ej,  and  Ej  itself.  The  parameters  k0dp  =  2.0,  e/e0  =  1.0, 
ji/jiQ  =  1.0  and  ak  =  0.0.  The  Field  maps  are  drawn  on  the  xz  domain  0  5  kx/dp,  kz/dp  <  10.0. 
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Figure  3  Same  as  Fig.  1.  except  =  5.0  and  ok  =  0.01 . 
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Figure  4  Same  as  Fig.  1,  except  e/e0  =  5.0  and  ak  =  0.1. 
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Figure  5  Magnitudes,  from  left  to  right,  of  the  x-  and  y-  directed  components  of  E-p  when  Ej  is  given  by 
Eq.  (19).  The  parameters  kodp  =  2.0,  e/e 0  =  5.0,  p/Po  =  1.0  and  ak  =  0.1.  The  field  maps  are 
drawn  on  the  xz  domain  0  S  kx/dn,  kz/dn  £  10.0. 
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Figure  6  Magnitudes  of  the  various  components  of  Ej  when  Ei  is  given  by  Eq.  (21).  From  left  to  right  in 
each  row:  magnitudes  of  the  x-  and  y-  directed  components.  From  top  to  bottom:  LCP  part  of  Ej, 
RCP  part  of  Ej,  and  E-p  itself.  The  parameters  l^dp  =  2.0,  e/e0  =  5.0,  p/p0  =  1.0  and  ak  =  0.1. 
The  field  maps  are  drawn  on  the  xz  domain  0  S  kx/dp,  kz/dp  S  10.0. 
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therefore,  that  if  ak  «  1.0,  then  Icr  ■■  kp  «  k  and  the  effect  of  the  chirality  parameter  a  is  slight  on  the 
field  Ep  as  can  be  seen  from  comparing  Fig.  3  with  Fig.  2.  But  as  ak  increases,  then  kg  and  differ 
widely  from  k  [9],  and  the  interference  of  the  LCP  and  the  RCP  parts  of  Ep  is  all  too  visible. 

Some  idea  of  this  interference  can  be  drawn  from  considering  only  specific  values  of  k  in  Eq.  (18])  or, 
equivalently,  by  modifying  the  line  source  field  Ej  of  Eq.  (14)  to 

EI  =  jjiH0(klr-rpl)8(K)  ,  (19) 

where  S(*)  is  the  Dirac  delta  function.  Shown  in  Fig.  S  are  the  magnitudes  of  the  x-  and  the  y-  directed 
components  of  Ep,  with  £/e0  =  5.0  and  ak  =  0.1,  there  being  no  z-directed  component  of  Ep.  If  Eq.  (19])is 
used  as  the  incident  field,  then,  from  the  preceding  theoretical  analysis,  it  is  easy  to  see  that 

Ej.  =  [2/ko(S  +  a^]  exp(j(kz-kodp)]  {-i  sin  (k2az)  +  j  cos  (k2az) }  ,  (20a) 

provided  ak  <  1  and  the  approximations 

kL * k [1  +  ak]  ;  kR»k[l- ak]  (20b, c) 

are  valid.  Fig.  5  vindicates  Eq.  (20a)  very  well  and  it  also  becomes  easy  to  observe  the  mutual  interference  of 
the  LCP  and  the  RCP  parts  of  Ep  in  deriving  this  equation. 

Finally,  in  Fig.  6,  the  incident  Held  is  modified  to 

Er=jrtH0(klr-rpl)5(K-ko/4)  ,  (21) 

and  E-p  is  computed;  e/e0  =  5.0  and  ak  =  0.1.  It  turns  out  that  P(ko/4)  **  QOcq/4),  while  Pl^o/4)/^  » 
PR(V4)/kR.  Consequently,  TlOcqM)  *■  ~aRT^(ko/4).  Hence,  the  RCP  and  the  LCP  parts  of  E-p  are 
virtually  identical  in  magnitude,  their  z-directed  components  are  negligible,  and  their  x-  and  y-  directed 
components  have  a  cos(k0x/4)  dependence  on  the  x  coordinate.  When,  the  LCP  and  the  RCP  parts  are  added 
up,  the  cartesian  components  of  E-p  have  well-defined  extrema  in  the  xz  plane.  In  view  of  Figs.  5  and  6, 
then,  the  variations  of  Ep  in  Fig.  4  can  be  easily  explained. 

It  should  be  noted  that  Ep(x,y)  in  Figs.  1  -  4  tends  to  decay  away  as  one  goes  farther  from  the  source  at 
rp.  This  is  natural  since  Ej,  in  the  vicinity  of  the  source,  has  reactive  components  with  large  magnitudes.  As 
one  moves  farther  and  farther  away  from  the  source,  these  reactive,  near-zone  field  components  tend  to  die 
out.  This  behavior  of  Ej  is  replicated  qualitatively  by  Ep  as  well,  a  phenomenon  which  has  been  noted  in 
other  electromagnetic  scattering  problems  also  [10]. 

In  summary,  it  has  been  shown  here  that  the  fields  refracted  into  a  chiral  half-space  are  markedly  different 
when  ak  is  substantially  different  from  zero,  i.e.,  when  ak  is  of  the  order  of  0.1.  A  realistic  source  has  been 
used  here,  as  apart  from  the  usual  analyses  involving  incident  planewaves  only  [6].  Provided  a  PWS 
decomposition  [3]  of  the  incident  field  is  possible,  this  procedure  can  be  extended  to  include  sources  of  other 
configurations  and  polarizations  as  well.  Such  an  analysis  would  be  of  use  to  various  researchers  in  the  area 
so  that  the  effect  of  the  near  fields  [10]  of  radiating  sources  may  not  be  ignored. 
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The  plane-wave  reflection  and  transmission  characteristics  of  bimaterial  interfaces  between  chiral  and  chiral- 
achiral  interfaces  have  been  extensively  explored.  We  report  on  the  curious  characteristics  of  the  interface  formed 
by  two  chiral  half -spaces,  one  of  which  is  the  mirror  image  of  the  other;  this  is  referred  to  as  problem  1.  It  is  shown 
that  these  characteristics  are  related  to  the  reflection  of  plane  waves  on  the  interface  of  a  chiral  half -space  and  a 
perfectly  conducting  one,  which  constitutes  problem  2. 


INTRODUCTION 

The  lack  of  geometric  symmetry  between  an  object  and  its 
mirror  image  is  referred  to  as  chirality,1-2  and  the  mirror 
image  of  a  chiral  object  cannot  be  made  to  coincide  with  the 
object  itself  by  any  operation  involving  only  rotations  and/or 
translations.  The  most  commonly  investigated  chiral  ob¬ 
jects  are  the  L-  and  the  D-type  stereoisomers,  which  are 
familiar  to  organic  chemists.  The  basis  for  the  difference  in 
the  physical  properties  of  the  mirror  conjugates  lies  in  the 
handedness,  or  the  chirality,  possessed  by  their  molecular 
configurations.  When  an  electromagnetic  disturbance  trav¬ 
els  through  a  medium  consisting  of  chiral  molecules,  it  is 
forced  to  adapt  to  the  handedness  of  the  molecules.  In  other 
words,  linearly  polarized  plane  waves  cannot  be  made  to 
propagate  through  such  a  medium,  whereas  left-circularly 
polarized  (LCP)  and  right-circularly  polarized  (RCP)  plane 
waves,  traveling  with  different  phase  velocities,  are  perfectly 
acceptable  solutions  of  the  vector  wave  equation  for  this 
class  of  medium. 

The  usual  constitutive  equations  D  =  «E  and  B  =  nH  do 
not  hold  for  chiral  media;  instead,  the  equations1-2 

D  =  tE  +  0«V  x  E,  B  »  *.H  +  X  H  (1) 

are  deemed  applicable,  with  /S  being  he  chirality  parameter. 
The  regular  time-harmonic  Maxwell  equations  (exp(-;wt)] 
are  now  utilized  along  with  Eqs.  (1),  and,  following  Bohren,3 
the  electric  and  the  magnetic  fields  are  transformed  into 

E-Qi  +  afiQs,  H  -  Qfl  +  aLQL,  (2) 

where  the  LCP  and  the  RCP  fields,  Q*,  and  Q«,  respectively, 
must  satisfy  the  conditions 

(V2  +  *,2)Qi.  -  0,  (V2  +  k22)QR  =>  0, 

VXQi-fcjQi,  V.Ql«0, 
VXQr--A2Qr,  V  ■  Qr  -  0.  (3) 

In  these  equations,  fci  -  fe/(l  -  k0)  and  k2 «  k/(l  +  kfi),  while 
k  -  «(«fi)1/2  is  a  convenient  abbreviation;  aR  ■»  - ;{m/<)i/2  is  an 
impedance,  and  ai  -  —j{t/n)xn  is  an  admittance. 
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The  plane-wave  reflection  and  transmission  characteris¬ 
tics  of  bimaterial  interfaces  between  chiral  and  chiral- 
achiral  interfaces  have  been  extensively  explored  recently.4-6 
Here  we  report  on  the  curious  characteristics  of  the  interface 
formed  by  two  chiral  half-spaces,  one  of  which  is  the  mirror 
image  of  the  other;  this  is  referred  to  as  problem  1.  We  show 
that  these  characteristics  are  related  to  the  reflection  of 
plane  waves  on  the  interface  of  a  chiral  half-space  and  a 
perfectly  conducting  one;  this  constitutes  problem  2. 

PROBLEM  1 

Let  the  space  z  <  0  be  occupied  by  the  chiral  medium  (e,  #t, 
0),  while  the  half-space  z  >  0  is  occupied  by  the  mirror- 
conjugate  medium  («,  n,  -0).  Either  a  LCP  or  a  RCP  plane 
wave  is  incident  upon  the  interface  z  -  0  from  the  zone  z  <  0. 
It  is  appropriate  to  express  the  fields  in  the  zone  z  <  0  by  the 
fields 

Q i  -  AL[ey  +  j(—eza1  +  ejxj/kJexpb’Ux  +  a,z)] 

+  Bd*y  +  /(e*«i  +  e^l/kJexpL/Ax  -  aqz)], 

z  <  0,  (4a) 

Q R  ~  AR[ey  +  j(eza2  -  ejxJ/kJexpb'Ocx  +  a2z)] 

+  #R[ey  ~  j(eza2  +  e^J/kJexpL/Ux  -  ar^z)], 

z  <  0.  (4b) 

In  these  equations,  Al  and  Ar  represent  the  incident  plane 
waves,  while  Bl  and  Br  are  the  amplitudes  of  the  reflected 
ones;  x  is  the  horizontal  wave  number  required  by  Snell’s 
laws  to  satisfy  the  phase-matching  condition  oh  the  inter¬ 
face  z  =»  0;  «i  —  +(ki2  —  k2)1/2  and  a2  m  +(k22  —  k2)l/i;  and  ex, 
etc.  are  the  unit  Cartesian  vectors. 

The  half-space  z  i  0  is  occupied  by  the  mirror-conjugate 
medium;  this  means  that  phase  velocities  of  the  LCP  and  the 
RCP  plane  waves  here  are,  respectively,  those  of  the  RCP 
and  the  LCP  plane  waves  in  the  medium  of  incidence.  Con¬ 
sequently,  an  acceptable  representation  of  the  fields  in  the 
medium  of  transmission  is  given  by 
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Q  L  m  Cdey  +  j(-exa2  +  e,K)/fe2]exp[/'(<cz  + 

z>0,  (5a) 

Q R  =  CR[ey+  jie^  -  erf/k^explHKX  +  oqs)], 

z>  0,  (5b) 

with  Ci  and  CR  being  the  transmission  coefficients. 

The  solution  of  the  boundary-value  problem  is  sought  in 
the  convenient  matrix  forms 


pi]  pLL 
[Br\  =  [Rri 

rlr~\  r  a<i 

(6a) 

[c:H£ 

Tulpil 

(6b) 

with  the  R’s  constituting  the  reflection  matrix  and  the  Ts 
constituting  the  transmission  matrix.  The  utilization  of 
Eqs.  (2),  (4),  and  (5)  in  ensuring  that  there  are  no  discontin¬ 
uities  in  the  tangential  E  fields  and  the  tangential  H  fields 
across  the  interface  z  -  0  yields 

Rll  —  ~Rrr  —  («1^2  —  ®2^i)/(«i^2  a2^l)> 

Tn  =  2alk2/(alk2  + 

Trr  m  2a2kl/(alk2  +  ot2k i), 

rlr  =  rrl  m  Tlr  «  Trl  =  0.  (7) 


This  curious  result  should  be  noted:  If  the  incident  plane 
wave  is  LCP  (RCP),  then  the  reflected  and  the  transmitted 
waves  are  also  LCP  (RCP).  There  are  no  waves  of  the 
opposite  handedness  generated  at  the  planar  boundaries 
between  mirror-conjugated  chiral  media.  Thus,  the  ar¬ 
rangement  of  problem  1  acts  somewhat  like  a  beam  splitter; 
an  incident  LCP  (RCP)  plane  wave  is  broken  into  two  LCP 
(RCP)  plane  waves,  which  leave  the  interface  in  opposite  z 
directions  with  different  amplitudes  and  phase  velocities. 

PROBLEM  2 

Let  the  medium  in  the  zone  z  >  0  be  perfectly  conducting. 
Then  a  solution  of  the  form 

is  to  be  sought.  For  this  purpose,  Eqs.  (2)  and  (4)  are 
utilized  to  ensure  the  nulling  of  the  tangential  E  field  at  the 
impenetrable  surface  z  *  0,  and  the  result  obtained  is 

rn=  -rRR  *  (ajfej  -  a2k1)/(«1k2  +  a2*i)» 

rflL=  2aLa,A2/(«1A2  +  ofjfcj), 

rLR=-  -2 aRa2kJ(axk2  +  a2kx).  (9) 

We  observe,  therefore,  that  when  a  LCP  or  a  RCP  plane 
wave  hits  a  perfectly  conducting  surface,  the  reflected  field 
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Fig.  1.  Illustration  of  the  correspondenM  between  problems  1  and  2  in  accordance  with  Eqs.  (10).  Problem  1  considers  the  interface  of  a  chin 
medium  and  a  perfect  conductor,  while  in  problem  2  the  two  chiral  media  on  either  side  of  the  interface  are  mirror  conjugates. 
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consists  of  components  of  both  circular  polarization  states. 
However,  when  *  -  ku  then  rRL  =  0,  rLL  -  -1,  rLR  -  -2 aR, 
and  rRR  —  1,  which  gives  rise  to  the  trivial  case  of  no  reflec¬ 
tion  for  the  grazing  LCP  incidence.  Likewise,  when  k  =  k2, 
then  rLR  -  0,  rRR  -  -1,  rRL  *  2 aL,  and  rLL  =  1. 

DISCUSSION 

The  correspondence  between  the  solutions  [Eqs.  (7)  and  (9)] 
of  problems  1  and  2  should  be  noted.  That  is, 

R-LL  m  r LI’  ® HR  m  r flfi" 

Tll  =  rnJaL>  Trr  *  -rLR/aR.  (10) 

For  illustration,  let  us  consider  the  case  of  LCP  incidence. 
In  problems  1  and  2,  ru  is  the  amplitude  of  the  reflected 
LCP  wave  that  travels  with  a  phase  velocity  wiki.  In  prob¬ 
lem  2,  rRi  is  the  amplitude  of  the  reflected  RCP  wave  with  a 
phase  velocity  u>/k2,  but  in  problem  1,  rRJai  is  the  amplitude 
of  the  transmitted  LCP  wave,  which  also  travels  with  a 
phase  velocity  w/fe2  because  the  medium  of  transmission  is 
the  mirror  conjugate  of  the  medium  of  incidence  and  reflec¬ 
tion.  Analogous  comments  also  apply  to  the  case  of  RCP 
plane-wave  incidence,  and  both  cases  are  schematically  il¬ 
lustrated  in  Fig.  1.  It  should  be  noted  that  Fermat’s  princi¬ 
ple  is  equally  well  satisfied  in  problems  1  and  2,  and  in 
identical  fashion.  -  _  ~ 

The  correspondences  given  in  Eqs.  (10),  coupled  with  the 
complementary  relations 


r ll(~P)  m  rRR(f)),  rRR(~fi)  ”  r u,(0) > 

“  aRrRL(P),  -  aL2rLR(0)  (11) 

point  to  the  imaging  concept.  Consider  problem  2  again 
with  an  incident  LCP  plane  wave  of  unit  amplitude.  Rele¬ 
vant  to  the  zone  z  <  0,  this  problem  is  equivalent  to  the 
superposition  of  two  problems,  each  of  which  is  a  specializa¬ 
tion  of  problem  1.  These  two  problems  are  (i)  a  problem  1  in 
which  a  LCP  plane  wave  is  incident  upon  the  interface  from 
the  zone  z  <  0  with  unit  amplitude  and  (ii)  a  problem  1  in 
which  a  RCP  plane  wave  is  incident  upon  the  interface  from 
the  zone  z  S  0  with  an  amplitude  equal  to  a£.  The  case  of  an 
incident  RCP  plane  wave  in  problem  2  can  also  be  handled  in 
this  way,  and  both  cases  are  schematically  illustrated  in  Fig. 
2.  However,  as  becomes  clear  from  Fig.  2,  the  use  of  an 
imaging  theory  for  chiral  media  is  complicated  for  scattering 
problems  in  general:  not  only  do  the  sources  get  imaged  but 
the  medium  does  also. 

The  authors  are  also  with  the  Center  for  the  Engineering 
of  Electronic  and  Acoustic  Materials,  The  Pennsylvania 
State  University,  University  Park,  Pennsylvania. 
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Two  approaches  involving  the  spatial  and  temporal  Fourier  transforms  have  been  used  to 
derive  time-  and  space-dependent  Green’s  functions  pertinent  to  the  propagation  of  sound 
waves  in  a  fluid  that  is  moving  with  a  constant  velocity  v.  The  two  approaches  give  rise  to 
differing  interpretations  of  the  observations  made  by  a  stationary  observer  vis-a-vis  those  made 
by  an  observer  moving  with  the  fluid.  The  properties  of  the  causal  and  the  noncausal  Green’s 
functions  are  analyzed,  and  are  shown  to  be  equivalent. 

PACS  numbers:  43.20.Rz,  43.30.Es 


INTRODUCTION 

It  is  readily  apparent  that,  at  the  point  of  observation, 
the  acoustic  frequency  spectrum  of  a  radiating  source  is  al¬ 
tered  by  the  flow  of  the  fluid  in  which  the  source  and  the 
observer  are  embedded;  this  phenomenon  has  been  the  sub¬ 
ject  of  intense  investigation.1-*  In  spite  of  the  widespread 
interest,  particularly  in  the  area  of  fluid  mechanics,5-7  we 
have  not  been  able  to  come  across  a  systematic  derivation  of 
the  pertinent  Green’s  function;  although,  it  should  be  noted 
that,  for  the  corresponding  problem  in  electromagnetic  field 
theory.  Green’s  functions  are  commonly  available.®-'0 

Here,  we  will  utilize  spatial  as  well  as  temporal  Fourier 
transforms  in  order  to  derive  two  separate  expressions  for 
the  Green’s  function  pertinent  to  the  radiation  and  propaga¬ 
tion  of  sound  waves  in  an  ideal  fluid  flowing  with  a  constant 
velocity  v.  These  two  approaches  involve  interchanging  the 
order  in  which  the  spatial  and  the  temporal  transforms,  and 
correspond  to  Doppler  shifts  in  either  spatial  or  temporal 
frequencies.  Thus  the  two  approaches  give  rise  to  differing 
-interpretations  of  the  observations  made  by  a  stationary  ob¬ 
server  vis-d-vis  those  made  by  an  observer  moving  with  the 
fluid.  It  is  shown  that  the  causal  Green’s  functions  derived 
from  the  two  approaches  are  identical. 

I.  PRELIMINARIES 

As  per  Morse  and  Ingard,11  in  a  fluid  moving  with  a 
velocity  v,  the  acoustic  pressure  p(r,t)  satisfies  the  wave 
equation 


[v2-c-2^  +  t*vJ  p{ r,/)=0 


within  a  source-free  region,  c  being  the  speed  of  sound  in  the 
medium  at  rest  (i.e.,  when  ▼  =  0);  this  equation  is  valid  in  a 
coordinate  system  that  is  moving  with  a  velocity  —  v  with 
respect  to  a  parallel  coordinate  frame  affixed  to  the  fluid. 
The  first  coordinate  system  (xj/p;,t)  is  called  stationary, 
while  another  coordinate  system  {xyj,t ' ) ,  to  be  discussed 
later,  will  be  referred  to  as  the  auxiliary  system;  the  velocity 
▼  ■*  oe,  in  the  stationary  system,  e,  being  the  unit  vector. 


The  corresponding  Green’s  function,  therefore,  must  be  the 
solution  of  the  inhomogeneous  wave  equation 

[v2  -  c-2  (iL  +  v-v)2]  g(r,t)  =  -  d(r)<5(r),  (2) 

6(‘)  being  the  Dirac  delta  function. 

As  Jackson12  has  pointed  out,  ( 1 )  is  not  invariant  under 
Galilean  transformations,  and  no  kinematic  transformation 
of p(r,t)  can  eliminate  the  term  v»V  in  ( l).  The  solution  of 
(2)  can,  however,  be  very  easily  attempted  if  the  term  (d  / 
dt  -f  vV)  can  be  simplified,  and  Fourier  analysis  comes  in 
very  useful  for  that  purpose.  We  define  two  sets  of  Fourier 
transform  pairs:  one  in  the  time  o  temporal  frequency  do¬ 
mains  specified  by  the  twin  sets  of  relations13 

a(r)  =  (2ir)-1  J  deocxpi  —  iat)a(a),  (3a) 
a{co)  =  J  d/exp(Kuf)a(r),  (3b) 

co  being  the  temporal  frequency;  and  the  second  in  the  space 
O  spatial  frequency  domains  given  as14 

a(r)  =  (2rr)-3  J  jj  d3q  exp(/qr)a(q),  (4a) 

a(q)=j"  J"J  d3rexp(  —  iqT)a(r),  (4b) 

in  which  q  is  the  (vector)  three-dimensional  spatial  frequen¬ 
cy  variable.  Armed  with  these  two  sets  of  Fourier  transform 
pairs,  we  can  solve  (2)  via  any  one  of  the  two  approaches 
possible. 


II.  DOPPLER  SHIFT  IN  THE  TEMPORAL  FREQUENCY 

The  first  approach  involves  taking  the  temporal  Fourier 
transform  of  (2),  which  yields  the  relation 

[V2  — c-2(  —  ua  +  vV)2]g(r,a>)  =  -5(  r),  (5) 

after  using  the  integral  representation  of  S(t)  given  as 
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<5(0  =  (2ir)~l  J  d<uexp(  -  itot).  (6) 

Next,  we  employ  the  factorization4-15"17 

g(  r,o))  =  A(r,<u)exp[  —  i(a/c)My2z],  (7) 

which  is  nothing  but  an  expression  of  the  Doppler  shift  in  the 
temporal  frequency  domain.  In  (7),  the  Mach  number 
M  —v/c,  and  y2!  1  —  M2)  =  1;  we  restrict  M <  1. 

As  a  consequence  of  (7),  (5)  simplifies  to 

[V1  —  (A fepV)2  +  (y<y/c)2]A(r,<a) 

=*  —  ^(rJexpt/fo/cJA/y^z].  (8) 

Utilization  of  (4a)  and  (4b),  along  with  the  representation 

<5(r)  —  (2ir)~3  J*  JJ  d3q  exp(iq*r)  (9) 

for  the  Dirac  delta  function  <5  ( r ) ,  in  ( 8 )  then  gives  rise  to  the 
algebraic  relation 

A(q,<y)  =  —  (c/y)2/W-  (c/y)2 [q*q  -  Af^e^q)2]}. 

(10) 

As  a  first  step  towards  obtaining  h(r,t),  we  take  the 
inverse  temporal  Fourier  transform  of  ( 10);  the  application 
of  (3a)  to  (10)  yields13 


A(q,f)  =  (2 ir)~'(c/y)  [irsgn(r)] 

X  [q*q  —  A/2(er*q)2]~,/2 
Xsin{(cr/7')V[q-q  -  Af  ■i(e1-q)'i]  >,  (11) 

whence,  on  taking  the  inverse  spatial  Fourier  transform  also, 
we  get 

h(r,t)  =  (2v)~4(c/y)  [?r  sgn(f)]/(Af;r,t),  (12) 

in  which  the  signum  function  sgn(|)  =  2U(£)  —  1,  whereas 
the  Heaviside  function  17(£)  =  1V|>0  and 
U(£)=0V£<0. 

On  expanding  q  in  a  cylindrical  coordinate  space 
(tfpWVtfx).  =  q*e*.  the  integral  I(M;r,t)  in  ( 12)  can  be 
given  as 


J(M;r,t) 


dqt  exp (iq,z) 


X  J  dtp,  exp(  iqpp  cos  <pq ) . 


(13) 


The  integral  on  the  azimuthal  angles,  is  easily  implemented 
via  an  integral  representation  of  the  cylindrical  Bessel  func¬ 
tion />(£>•• 


rlw 

MS)  =  (2ir)-1  J  da  exp  (if  cos  a) ; 


and  the  consequent  integral  on  can  be  calculated  through 

Sonine’s  discontinuous  integral15: 


J“  dni9{cg)m(«ir+F)(S*  +  b2rl/2 

=  cos (by]al-ci)(a2-c1)-'nU(a-c),  c>0. 

(14b) 

As  a  result, 

=  (2jt)  U  (-1)  J"  dq, 

X  exp  ( /$,z )  cos  ( y” 1  f0) ,  (15a) 


xK‘+t)+,('"t)]-  (isi” 

on  utilizing  (6);  in  these  equations, 

t0  =  sKct/y)1  - p 2,  (16a) 

p  =  +JP~+~jr.  (16b) 

Further  simplification  is  afforded  by  the  fact  that 

2|a|<5(£2-a2)=<5(|-+a)-|-<5(£-a),  (17) 

and  that  ct/y>p  due  to  the  Heaviside  function;  consequent¬ 
ly. 

I(M;r,t)  =  (4ir2/y)U(cr/y— p)S(z1  —  y~2t\). 


h(r ,/)  =  (c/4y2tr2)[trsgn(t)]U(ct/y—p) 

X<5  (z2-y"2r5),  (18b) 

From  (3)  and  (7),  it  is  easy  to  see  that 
g(r,/)  =  A(r,r  +  Me_yz);  (19) 

hence,  a  solution  of  (2),  to  be  subscripted  t  for  later  use,  is 
given  as 

g,(r,t)  =  (c/4a2)[vsgn(t  +  Mc~iy2z)] 

XU\c(t  +  Afc~,y2z)/y — p] 

XSir2  —  2zvt  —  c*y~2t2).  (20) 

To  be  noted  is  the  fact  that  the  inversion  (13)  of  (12)  had 
not  assumed  causality;  if  causality  is  to  be  incorporated 
when  the  inverse  temporal  Fourier  transform  of  ( 1 1 )  is  tak¬ 
en,13-19  then  the  signum  function  sgn(/  +  Mc~  ly2z)  in  (20) 
should  be  replaced  by  the  2U(  t  -f  Mc~  xyLz).  Here  and  here¬ 
after,  the  term  causal  should  be  interpreted  as  causal  will, 
respect  to  time. 


III.  DOPPLER  SHIFT  IN  THE  SPATIAL  FREQUENCY 

The  second  approach  mirrors  the  first  one  in  that  tht 
order  of  taking  the  temporal  and  the  spatial  Fourier  trans 
forms  is  reversed.  Thus,  on  using  (4)  and  (9)  in  (2),  w* 
obtain 

=  —  S(t).  (21 

The  introduction  of  the  factorization 

£(<LO  =/(q>t)exp(  -  te,*qi>0  (22 

in  (21)  results  in 


1853 


J.  Acoust  Soc.  Am.,  Vol.  85,  No.  5,  May  1 989 


Lakfitakia  etei.:  Sound  In  moving  fluid  185C 


^  _  q-q_  C-J  =  -  $(r)exp(iex*qi/r). 

(23) 

It  should  be  noted  that  (22)  tantamounts  to  a  Doppler  shift  in 
the  spatial  frequency  domain. 

The  temporal  Fourier  transform  of  (23 )  gives  rise  to  the 
algebraic  expression  for/tq,a>)  given  by 

flq,a)  =  —  <f/{o2  —  q-qc1).  (24) 

As  a  first  step  towards  deriving/(r,r),  the  inverse  temporal 
Fourier  transform  is  taken;  hence,13 

f(q,t)  —  (2jr)-,c2[jrsgn(/)](qc)_l  sin(?cf),  (25) 

from  where,  on  taking  the  inverse  spatial  Fourier  transform 
also,  we  get 

f(r,t)  —  (2ir)  “■*<:[«•  sgn(r)  ]/(M;r,r),  (26) 

in  which  the  integral  J(M;r,t)  is  given  by 

J(M;r,t)  =  J  dqtexp(iqIz)  dqp  q„  (<fp  +  qj ) " 1/2 

Xsin (ct^fi  +  q£ ) 

X  f  dq>q  exp(iqpp  cos  q>4).  (27) 

Jo 

The  use  of  ( 14a),  ( 14b)  and  ( 17)  in  (27)  finally  yields 
J(M;r,t)  =  4ir2U(.ct-p)S(rs  -  (ft2).  (28a) 

whence 

f{r,t)  =  (c/Aifl)[irs^n{t)]U{ct  —  p)8{f  —  <ft2). 

(28b) 

Finally,  by  applying  (22)  in  (4),  it  can  be  shown  that 
g(r,t)  =flr-it,t),  (29) 

and  the  use  of  this  relationship  in  (28b)  gives  the  Green’s 
function 

g,(r,f)  =  (c/Air2) 

X  [ir  sgn(f)  ]  U{ct  —  p)8{fl  -  2 zvt  —  (fy~zt z); 

(30) 

the  subscript  r  on  g  in  (30)  is  for  later  use.  Again,  we  note 
that  gr  derived  above  is  noncausal;  for  causality  to  be  incor¬ 
porated,  sgn(f)  should  be  replaced  by  2U(t). 

IV.  DISCUSSION 

The  two  Green’s  functions  derived  above  share  one  im¬ 
portant  feature,  i.e.,  they  possess  the  same  singularities.  The 
singularities  of  both  g,  (r.t)  and  g,  (r ,t)  are  embedded  in  the 
delta  function  that  can  be  expanded  as 

6(t*  —  2zvt  —  cfy~2t 2) 

=  {2yjR  2  +  Izut  4-  <fy~2t z )  “ 1 

X  l<5(z 4-  y-'JR*  +  2zvt  +  c1y-iti) 

+  8{z  —  y~  lyJR  +  2 zvt  +  fy~lt  *)  ],  (31a) 

where 

R1^M1fzl-p2.  (31b) 

In  this  respect,  therefore,  there  is  nothing  to  distinguish 


between  g,  {r.t)  and  gr  (r ,t).  In  addition,  it  should  be  noted 
that  if  v  =*  0,  then  both  g,  and  gr  reduce  to  the  familiar 
expression 

g,(r,r)|M_0  =&(r,f)|,_o 

=  ( l/Aifr)  [rr  sgn(t)  ]  U(ct  —  p) 
X[6(t-rc-l)+8(t  +  rc-1)],  (32) 

which  too  can  be  made  causal  if  sgn(  t)  is  replaced  by  2(7(0 
and  the  second  delta  function  is  omitted.  However,  the 
switching  functions,  i.e.,  the  signum  and  the  Heaviside  func¬ 
tions,  included  in  (20)  and  (30)  have  different  arguments. 

The  major  reason  for  the  difference  between  gr  ( r.t)  and 
g,(r,t)  appears  to  be  due  to  the  specific  factorizations  (7) 
and  (22),  respectively.  As  stated  earlier,  while  (7)  expresses 
a  Doppler  shift  for  time-harmonic  waves,  (22)  gives  rise  to  a 
Doppler  shift  for  space-harmonic  fields.  Howevef,  it  is  easy 
to  show  that 


g,  {r.t)  =  (2 v)~*  |"  da  exp(  —  iat)ex p  ^ 


■ ur 


d3q  exp  ( /q*r )  A  ( q,a ) 


(33) 


is  a  proper  solution  of  (2)  by  substituting  (33)  in  (2)  and 
noting  that  the  differential  operators  of  the  left  side  of  (2) 
will  act  only  on  the  exponential  terms  of  (33).  Moreover,  it 
is  also  possible  to  show  similarly  that 

h(r.t)  =  (2ir)~4  J  da  exp(  —  iat) 

xj  JJ  d3qexp(/q-r)A(q^>)  (34a) 
satisfies  the  differential  equation 

(34b) 

It  should  be  noted  that  by  using  ( 19)  in  (2),  it  is  observed 
that 


[(S)+(SM3)  . 

-  (-^)J  (~4)]  Hx,yS,n  =  ~8{r)8{t),  (34c) 


where  t'  =  t  +  M-fz/c,  f  =  z,  and  it  is  assumed  that  the  sta¬ 
tionary  and  the  auxiliary  frames  coincide  at  f  =  f '  =  0.  Like¬ 
wise,  it  can  also  be  shown  that 


gr(r,t)  -  (2ir) 


- ur - 


q  exp(/q«r)exp(  —  /q*vt) 


xj*  do  cxp(  —  iat)f(qjD)  (35) 


is  also  a  proper  solution  of  (2)  by  substituting  (35)  in  (2).  In 
addition. 


/tr,r)  =  (2r)-4J  JJ" 


d3qexp(/qf) 


da  exp(  —  iat)f{q,a) 


(36a) 
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satisfies  the  differential  equation 

(v2  -  c"2  j^/tr.r)  =  -  <5(05(0.  (36b) 

Again,  it  is  to  be  noted  that  by  substituting  (29)  in  (2),  and 
assuming  that  the  stationary  and  the  auxiliary  frames  coin* 
dde  at  z  =  z'  —  0,  that 


Kdy2 

-5(05(0. 


(36c) 


where  z*  —  z  —  vt  and  t  '  =  t.  Thus  the  proprieties  of  the  two 
procedures  are  in  good  standing. 

The  equivalence  of  g,(r,t)  and  g,(r,t)  is  now  demon¬ 
strated.  Using  the  relations  (4b)  and  (7),  it  can  be  shown 
that 

g,  (q,<y)  =  A(q  4-  toy/y/tr.cv),  (37a) 

while,  from  (3b)  and  (22), 

gr(.q,co)  =ytqIo  —  q-v).  (37b) 

However,  the  right  sides  of  (37a)  and  (37b)  are  identical  as 
can  be  seen  by  utilizing  (10)  and  (24).  Therefore,  in  the 
temporal-spatial-frequency  domain,  gr(q,<u)  and  g,( q,a>) 
are  identical,  as  should  be  expected.20  As  a  result  of  the 
uniqueness  theorem  for  the  inverse  Fourier  transforms,21  it 
follows  that  gr  (r ,/)  =  g,  (r,r)  except  perhaps  at  their  singu¬ 
larities.  It  is  obvious  that  everywhere  except  at  their  (com¬ 
mon)  singularities,  gr  (r,t)  —  g,  (r.r)  =  0  as  a  consequence 
of  the  delta  functions  in  (20)  and  (30). 

Finally,  we  show  that  the  two  causal  Green's  functions 
have  identical  values  at  their  singularities.  Analyzing,  first, 
for  the  causal  g,  (r ,r) ,  it  should  be  noted  that 

r*  —  2 ror  —  <?Y~xt 2  =p2  +  yPz2  —  (c/y)2(t  +  M-fz/c)2 

(38) 

and />> 0.  Therefore,  at  the  singularities  p2  -f  y^z2  =  (c/y)2 
X  (t  +  My^z/c)2;  together  withp>0,  this  implies  that  (c/y) 
X  (t  -F  MyPz/c)  ~>p.  As  a  result,  at  its  singularities,  the  caus¬ 
al  g,  (r.O  =  c/2ir.  Next,  for  the  causal  g,  (r ,t)  it  should  be 
noted  that 

r1  —  2zvt  —  c2y~2t2=p2  +  (z  —  vt)2  —  c2*2.  (39) 

Therefore,  at  thesingularities,/?2  +  ( z  —  vt )2  =  c?t 2;  togeth¬ 
er  with /?>0,  this  implies  that  cr>p.  Consequently,  at  its  sin¬ 
gularities,  the  causal  g,  (r.r)  =  c/lv.  Hence,  the  causal 
g,(r,r)  andg,(r,r)  are  equal  for  all  rand  r>  0.  Similar  reason¬ 
ing  also  gives  the  noncausal  g,(r,r)  —gT(r,t)  =c/4ir  at 
their  singularities,  consequently,  the  noncausal  g,  (r,r)  and 
j;  (r ,r)  are  also  equal  everywhere. 

From  the  available  literature,  it  is  observed  that 


t,(T/a)  -  (^jexp^-^VI*4  +  /  +  Y1*1) 


Xexp£  —  Afy*zj  (x2  +jp  +  yV)"1'2 


(40) 


has  been  commonly  utilized  for  solving  time-harmonic  prob- 
lona.4-1*~17  On  the  other  hand,  g,  (r ,/)  is  in  a  form  more 
amenable  to  interpretation  for  time-dependent  problems.  In 
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order  to  illustrate  the  utility  of  these  Green’s  functions,  we 
consider  a  source  pressure  which  acts  on  a  ring  of  radius  a. 

First,  let  the  time-harmonic  problem  be  considered,  i.e., 
the  source  pressure 

p,(t,a)  =  (P/a)S(p  —  a)5(z)ex  p(  —Uot),  (41) 

P  being  some  constant.  The  radiated  pressure  will  be  isotrop¬ 
ic  in  the  x-y  plane;  hence,,  in  the  x-z  plane  it  will  be  given, 
using  (40)  and  (41 ),  as 

Pnd  (*,CU,<u) 

=  (-^-)«P(  -  ^t) exp |  -  i My>z\  J  d<p0 

Xexp [/(jWcJVx2  +  a2  —  lax  cos  <p0  +  yV] 

X  (x2  +  a2  —  lax  cos  <pQ  +  yV)  ~ l/!,  x^a.  (42) 

The  evaluation  of  this  integral  must  be  done  numerically; 
however  Fraunhoffer-type  approximations14  can  be  made  to 
obtain  the  far-zone  pressure  quite  easily.  If 
(x2  -(-  jrz2) 1/2 ->  oo,  then  the  integrand  can  be  simplified  to 

exp  [  /( ytu/c)  (x2  +  y 2z2 ) 1  /2  ]  (x2  +  /z2 ) " 1  /2 

Xexp[  —  i (yeu/c) ax  cos  tp^x2  -f  yV)“,/2]. 

In  which  case,  (42)  simplifies  to 

Pnd  (x,0,z,<y) 

=  (Py/  2)exp(  —  /<yf)exp[  —  i(a/c)MyLz\ 

Xexp[;(y<u/c)(x2  +  y2z2)1/2](x2  +  y2z2)~ul 

'X.JQ[(ya>/c)ax(x1  +  y3z2)~'12].  (43) 

To  illustrate  the  use  of  gr  (r,r),  we  consider  the  same 
ring  source  with  an  impulse  excitation;  i.e. 

P,(r,t)  =  (P/a)S(p  -  a)S(z)S(t).  (44) 

In  the  equatorial  plane  z  =  0,  the  radiation  is  isotropic. 
Hence,  from  the  causal  version  of  ( 30 ) ,  the  radiated  pressure 
will  be  given  by 

^U0.0,0  =  (^)f^o 

X6(ct  —  yJ?~+lF^lAxco$p^) 

Xix2  +  a*  —  lax  cos  <p0)~ 1,2 ,  t> 0.  (45) 

Thus  p^  (x,0,0,r)  exists  only  for  y(x  —  a)<kct<y(x  +  a) 
and  has  a  magnitude  of  (PyP/2ir)/t  during  that  interval.  On 
the  other  hand,  along  the  z  axis,  the  radiated  pressure  is 

Pnd  (0,0 ,z,r)  =  Pcy 2  5(ct  -(-  zMy2 

—  Vz2 M 2y*  +  yV  +  yV).  (46) 

It  should  be  noted  that  the  time  of  arrival  at  the  origin  from 
(46)  equals  t  =  ( y/c)a ;  thus  the  radiated  pulse  arrives  there 
later  in  a  moving  fluid  than  in  a  stationary  fluid.  Also, 
t  =  a/c  is  time  of  arrival  at  z  =  Ma,  while  t  —  (a/ 
c)y2(  1  4-  M 2)  is  the  time  of  arrival  at  z  =  —  Ma  on  the  z 
axis,  which  conclusions  correspond,  respectively,  to  the  up¬ 
stream  and  the  downstream  velocities  of  z-directed  plane 
waves  derived  by  Ingard  and  Singbal  for  moving  fluids.20 
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Time-harmonic  and  time-dependent  Green’s  functions 
are  derived  for  a  lossless,  uniaxial  gyroelectromagnetic  me¬ 
dium  whose  permeability  tensor  is  a  scalar  multiple  of  its 
permittivity  tensor,  and  their  properties  are  nvestigated. 
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The  derived.  Green’s  functions  can  be  used  for  the  solution 
of  initial  and  boundary  value  problems,  as  well  as  for  ob¬ 
taining  the  electromagnetic  fields  radiated  by  electric  and 
magnetic  sources. 

Researchers  have  begun  to  focus  attention  on  the  electro¬ 
magnetic  theory  of  anisotropic  media  because  of  the  recent 
proliferation  in  the  use  of  such  materials  for  a  variety  of 
applications.  Post1  laid  the  foundation  of  the  constitutive 
equations  of  a  general,  linear,  bianisotropic,  homogeneous 
medium  around  a  quarter  of  a  century  ago;  and  his  work  was 
later  extended  by  several  workers  (see,  e.g.,  Refs.  2-8).  In 
that  vein,  the  dyadic  formalism  promoted  by  Chen9  in  con¬ 
nection  with  his  work  on  the  time-harmonic  electric  Green's 
function10  for  uniaxial  dielectrics  has  recently  been  utilized 
to  compute  the  corresponding  magnetic  Green's  function 
also11;  furthermore,  time-dependent  Green’s  functions  have 
been  derived  for  the  same  media.1'2 

Gyroelectromagnetic  media  have  not  been  given  much 
attention.  Ray  optics  of  crystals  such  as  calcite  and  rutile  is 
well  known,13  and  some  work  on  radiation  due  to  Chow14  is 
also  available.  Some  results  can  be  obtained  by  appropriate 
simplifications  of  the  extant  results  pertaining  to  bianisotro¬ 
pic  media.2"8  In  general,  however,  a  systematic  electrody¬ 
namic  theory  of  such  media  is  not  well-developed,  particular¬ 
ly  with  respect  to  the  infinite  medium  Green’s  functions.  In 
the  sequel  we  have  derived  the  time-harmonic  as  well  as  the 
time-dependent  Green’s  functions  for  a  homogeneous,  uni¬ 
axial  gyroelectromagnetic  medium  whose  permeability  ten¬ 
sor  is  a  scalar  multiple  of  its  permittivity  tensor,  and  we  have 
explored  the  characteristics  of  the  derived  Green’s  functions. 
Following  a  proposal  due  to  Rumsey15  for  constructing  artifi¬ 
cial  uniaxial  gyromagnetic  materials,  it  is  conjectured  that 
materials  of  the  type  studied  here  can  be  constructed  by 
embedding  parallel  ferrite  and  dielectric  fibers  in  some  host 
medium.  It  is  to  be  noted  that  while  boldface  English  letters 
represent  vectors,  the  German  letters  denote  tensors  or  dy¬ 
ads;  in  addition,  all  vector  operations  have  been  ordered  to 
proceed  from  right  to  left. 

Uniaxial  gyroelectromagnetic  media  possess  a  single  optic 
axis  for  both  their  electric  and  magnetic  properties;  in  addi¬ 
tion,  we  assume  that  the  permeability  tensor  is  a  scalar 
multiple  of  the  permittivity  tensor.  Thus,  the  appropriate 
constitutive  equations  are  given  as 

D  =  c0a-E;  B  =  n0qa-B,  (la,b) 

in  which  the  uniaxial  tensor  a  is  specified  by9 

a  =  axa  +  (a,  -ax)e,ec,  (2) 

ec  being  the  unit  vector  parallel  to  the  optic  axis,  and  21,  the 
idempotent;  a  j_,  a\\,  and  q  are  scalar  constitutive  parameters 
assumed  constant;  and  co  and  mo  refer  to  free  space.  Substi¬ 
tution  into  the  time-dependent  Maxwell’s  equations  yields 
the  field  equations 

(V  X  a)  •  a"1  •  (V  X  8)  •  E  +  •  |d*E !dt2\ 

-  -MofllW/atl  -  V  X  (a'1  •  K),  (3a) 
(VXH).«-'.(VXJ).H  +  eoMo<?a  •  I^H/at2) 

-  -'0\dVJdt\  +  V  x  (a"1  •  J).  (3b) 

Likewise,  substitution  into  the  time-harmonic  Maxwell’s 
equations,  with  an  exp{— imt j  dependence,  gives  rise  to  the 
differential  equations 


(V  x  8)  •  a-1  •  (V  x  H)  •  E  -  k2a  •  E 

»  ivfiQqJ  -  V  X  (a-1  •  K),  (4a) 

(V  x  8)  •  a'1  •  (V  X  8)  ■  H  —  fc2a  •  H 

m  iut0 K  +  V  X  (a-1  •  J),  (4b) 

in  which  k  =  ajy(coMo<?)  has  been  used  for  convenience.  It 
should  be  noted  that  in  Eqs.  (3a, b)  and  (4a, b),  as  well  as 
hereafter,  J  and  K,  respectively,  denote  the  electric  and  the 
magnetic  source  current  densities. 

The  solutions  of  Eqs.  (3a, b)  and  (4a, b)  are  now  sought  in 
the  forms  of  the  infinite-domain  Green’s  functions,  begin¬ 
ning  with  the  time-harmonic  solutions. 

For  the  time-harmonic  case,  by  direct  substitution  in  Eqs. 
(4a, b)  it  can  be  easily  verified  that  the  solutions 

E(rl<i>)  =  S  jj  d3r,[i«n092S,(Rlu)  •  J(rju) 

~  93m(Rli*>)  •  la"1  •  K(r,|w)|],  (5a) 

H(rlui)  =  SSS  dVji^oSMRlu)  -  K(r,U) 

+  »„,(Rlw)  •  |<r‘  •  J(rju>)|]  (5b) 

hold,  provided  the  Green’s  functions  $«,  and  23m  satisfy, 
respectively,  the  equations 

[(V  X  8)  •  a"1  •  (V  x  8)  -  *2a]  •  ®„(RU)  *  85(R),  (6a) 

[(V  X  »)  •  a-1  •  (V  x  8)  -  fc2n)  •  ®m(Rl«)  =  V  X  ai(R),  (6b) 

S(-)  being  the  Dirac  delta  function.  It  should  be  noted  that 
the  appearance  of  R  =  r  -  r,  in  Eqs.  (6a,b)  signifies  the 
spatial  invariance  of  the  uniaxial  gyroelectromagnetic  medi¬ 
um  considered  here,  with  r,  being  the  source  point  and  r,  the 
field  point. 

To  solve  Eq.  (6a),  its  3-D  spatial  Fourier  transform  is 
taken,  as  has  been  done  elsewhere  for  uniaxial  dielectrics.9"11 
It  turns  out  that 

S8,(R|a>)  ■  (2t)-3  jfj  d3p  exp[ip  •  RlSB^'tolw),  (7) 
in  which  the  dyadic 

SSUpIw)  =  -(p  X  a)  •  a"1  •  (p  X  8)  -  k2a,  (8) 

and  p  is  the  3-D  spatial  frequency  vector.  Using  dyadic 
algebra,9  the  inverse  ©^(plw)  can  be  derived  as 

2BJl(pl<*j)  =  [a^cT1  -  pp/fe2]/[p  •  a  •  p  -  feV^a,].  (9) 

By  substituting  Eq.  (9)  into  Eq.  (7)  and  evaluating  the  result¬ 
ing  integral,  the  Green’s  dyadic  ®c(Rlu>)  is  calculated  to  be 

SB,(RU)  =  Ko.a"1  +  VV/k2}  exp[ikR']/4rR',  (10) 

in  which 

R,  -  v'[a^ajR  -  a"1  •  R).  (11) 

Similarly,  by  taking  the  Fourier  transform  of  Eq.  (6b),  it 
can  be  shown  that 

®m(Riu>)  »  (2t)"3  jj!  d3pexp[ip-R][ffl;'l(p|a)).ipxa], 

(12) 

whence, 

®m(R|o»)  -  [o2±a,a-‘l  •  V  X  8  eip[i'fcR,)/4rR, 

-  a\a,o''  X  V|exp|ifeR,]/4Tfi,|.  (13) 
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From  Eqs.  (10)  and  (13),  it  is  quite  clear  that  the  uniaxial 
gyroelectromagnetic  media  considered  in  this  paper  are  sin¬ 
gly  refringent,  as  opposed  to  uniaxial  dielectrics  which  are 
birefringent.  Additionally,  all  waves  in  these  uniaxial  gyro- 
electromagnetic  media  are  of  the  extraordinary  type,  where¬ 
as  both  ordinary  and  extraordinary  waves  can  propagate  in 
uniaxial  dielectrics.*-12 

Coordinate-free  forms9  of  33r(R|u)  and  33m(R|oj)  can  be 
easily  obtained  from  Eqs.  (10)  and  (13),  and  are,  respectively, 
given  as 

33,(Rlw)  -  (aaxo,  exp(ifeflj/4rfi,)((l  -  (ikRerl  +  (iWy-2|<rl 
-  [1  -  3(ikR')-'  +  3(ikRer2] 

(a-1  •  R)(a-1  •  R)/(R  •  a-1  •  R]),  (14a) 

®„(Rlw)  =  (*  Vx«?  explikR,]/4*Re) 

( UikS,)- 2  -  (ikRfrl]a-1  X  [a-1  •  R]).  (14b) 

These  forms  may  be  of  use  when  solving  radiation  and  scat¬ 
tering  problems  using  numerical  techniques  which  call  for 
the  discretization  of  surfaces  and/or  volumes. 

The  reciprocal  nature4  of  these  uniaxial  gyroelectromag¬ 
netic  media  is  reflected  by  the  symmetry  properties 

»,(RU)  -  a,(-Rl«);  »m(Rlu)  =  -a^-Rlu);  '  (15a, b) 

in  addition,  the  following  transpose  properties  also  hold: 

[a/Rltt))]"  =  S«(Rl«):  (16a) 

[SBJRU  •  a~T  -  -®m(RU>)  -  o-‘.  (16b) 

These  relationships  can  be  verified  very  easily  from  the 
coordinate-free  forms  (14a,b).  Furthermore,  using  Eq.  (10) 
it  can  be  shown  that  93e(Rloj)  and  93m(R|aj)  are  connected 
through  the  relationship 

®m(Rl«)  «  a'1  •  (V  x  ®((R|«)1  •  a;  (17a) 

and  the  combination  of  Eqs.  (6a)  and  (17a)  gives 

V  X  ®m(Rlu)  -  k2a  ■  ®,(Rlw)  -  a  =  a«(R).  (17b) 

Correct  to  the  order  R~3,  the  near-zone  approximations  of 
the  Green’s  dyadics  can  be  obtaineo  from  Eqs.  (14a,b)  as 

S8,(R|«)  -  (a2xa,  exp[ifefij/4rft,)(ifefl,)-2 

X  (a-1  -  3(<r‘  •  R)(a_1  •  R)/[R  ■  a-1  •  R]),  (18a) 

SS^Rlo,)  -  0.  (18b) 

The  far-zone  approximations,  correct  to  order  R~l,  are  de¬ 
rived  from  Eqs.  (14a,b)  as 

8,(R|«)  -  (aio,  exp[»ARJ/4*Re) 

X  (a-1  -  (a-1  •  R)(a-‘  •  R)/[R  -  a'1  •  RJ),  (19a) 
8„(Rl«)  -  -(A2a4xa}exp[»ARJ/4xR,) 

X  (i*Re)-‘(a-1  x  (a’1  ■  R]).  (19b) 

As  a  by-product  of  the  foregoing  analysis,  the  propagation 
characteristics  of  plane  waves  in  the  uniaxial  gyroelectro- 
magnetic  media  considered  here  can  be  obtained.  Consider 
a  plane  wave  traveling  in  the  direction  ep  with  the  phase 
velocity  w/p.  The  wavenumber  p(ep)  can  be  calculated  from 
the  relation 


P2(ep)  =  *2a5.aj/[ax (ep  X  ef)  •  (ep  x  ec)  +  a,(ep  •  e,)2).  (20) 

When  ep  is  parallel  to  the  optic  axis,  the  associated  E  and  H 
fields  are  mutually  orthogonal,  and  are  orthogonal  to  ep  as 
well,  and  any  polarization  can  be  chosen  subject  to  the  re¬ 
strictions  imposed  by  the  time-harmonic  Maxwell’s  equa¬ 
tions.  On  the  other  hand,  if  ep  is  not  collinear  with  ec,  two 
well-defined  polarizations  are  possible.  Respectively,  these 
are  given  as 

(i)  E:  =  ep  X  ec;  H,  »  -[A2ox/wmo<7P]E2; 

(21a, b) 

(ii)  E2  =>  ec  -  ep(ep  •  er)(p/feax)2;  H2  =  -{wtopWaJEj. 

(22a, b) 

It  should  be  noted  that  the  first  polarization,  Eqs.  (21a,b),  is 
TE-to-ep;  while  the  second  one,  Eqs.  (22a, b),  is  TH-to-ep;  in 
either  case,  the  B  and  D  fields  are  orthogonal  to  ep. 

The  solutions  of  Eqs.  (3a, b)  are  also  sought  for  in  the  forms 

E(rlt)  =  SSS  d\  S  <tt,|po98i(Klr)  •  J(r,lt,> 

+  ®2(R|r)  •  a-1  •  K(r,|t,)],  (23a) 

H(rlt)  =  ///  d\  f  dtJ[-®2(R|r)  •  a"1  •  J(rjt,) 

+  H)8I(RlT).K(r,M,  (23b) 

in  which  the  use  of  r  =  t  —  t,  is  in  accordance  with  the 
temporal  invariance  of  the  medium,  with  r  >  0  demanded  by 
the  consideration  of  causality9;  t  is  the  time  associated  with 
the  field  point,  while  t,  is  associated  with  the  source  point. 
By  substituting  Eqs.  (23a, b)  into  Eqs.  (3a, b),  and  then  taking 
the  (causal)  temporal  Fourier  transform  of  the  resulting 
equations,  it  can  be  shown  that 

S^Rlr)  “  (2x)-1  |  du(iu >)  exp[-i(i)t]®e(RU),  (24a) 

J—+& 

SB2(Rlr)  =•  -(2x)-1  [  du  exp[— ia>£]®m(Rl«).  (24b) 

J—+& 

Causality  is  obtained  by  invoking  analytic  continuation  from 
real  to  complex  w  by  moving  the  integration  path  above  the 
real  axis  in  the  complex  o>  plane;  therefore,  A  >  0  in  Eqs. 
(23a, b)  so  as  to  ensure  avoiding  the  pole  singularities.1617 
Finally,  by  using  the  unit  step  function  u(t )  defined  as18 

u(t)  *  +  (1/t)  j  (du/a)  sina>£,  f  2:  0 

•0,  £  <  0,  (25) 

as  well  as  the  Dirac  delta  function,  the  integrals  involved  on 
the  right-hand  sides  of  Eqs.  (23a, b)  are  evaluated  out  to  yield 

93,(R|t)  -  (a2xa,/4xRe)([-i'(rt)  +  (cJR,)i( r,) 

-  (C'/R,)2u(T,))a-'  -  [— i'(rt)  +  3(cJRe)S(rt) 

-  Z(c JRJhj. (r.JKa-1  •  R)(«-‘  •  R)/[R  -  a"1  •  RJ),  (26) 
®2(R|t)  -  (o4xaf/4TC^J)([S'(r,) 

+  (cy«.W(re)Ja-1X[«-1-Rl),  (27) 

in  which  5'(*)  is  the  unit  doublet  impulse17, 

c, - W?c0]"l/2-  Tt“T~  Rjcg,  (28a,b) 
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r,  being  the  retarded  time  in  the  direction  R/R.  The  time- 
domain  Green’s  dyadics,  Eqs.  (26)  and  (27),  once  again  show 
that  the  media  being  considered  here  are  singly  refringent,  as 
opposed  to  uniaxial  dielectric  or  uniaxial  magnetic  materials 
which  are  birefringent.12 

The  reciprocal  nature  of  these  uniaxial  gyroelectromag- 
netic  media  is  reaffirmed  by  the  symmetry  properties 

SBj(Blr)  -  a,(-R|r);  »2(R|t)  -  -®2(-R|r);  (29a, b) 

and  the  following  transpose  properties  also  hold: 

[SitRlrir  -  SS^Rlr);  (SB2(R1t)  •  a*1]"  =  -®2(R|t)  •  a'1. 

(30a,  b) 

The  wavehead  is  established  by  the  condition  re  =  0  and  is 
dependent  on  the  direction  relative  to  the  optic  axis;  while  at 
large  times  such  that  re  >  0,  the  Green's  functions  of  Eqs. 
(26)  and  (27)  can  be  reduced  to 

®i(Rlr)  -  (a\ai>/i-KRl)(cJRl)2{-a~1  +  3(a'‘  •  R) 

(a_1-R)/[R-a"I-RJ);  r,>0,  (31) 

®2(R|t)  -  0;  t e  >  0.  (32) 

In  summary,  time-harmonic  and  time-dependent  dyadic 
Green’s  functions  have  been  derived  for  a  lossless,  uniaxial 
gyroelectromagnetic  medium  whose  permeability  tensor  is  a 
scalar  multiple  of  its  permittivity  tensor,  and  their  proper¬ 
ties  are  investigated.  The  derived  Green’s  functions  can  be 
used  for  the  solution  of  initial  and  boundary  value  problems, 
as  well  as  for  obtaining  the  electromagnetic  fields  radiated  by 
electric  and  magnetic  sources. 

This  work  was  supported  in  part  by  the  Ben  Franklin 
Partnership  Program  of  the  Commonwealth  of  Pennsylva¬ 
nia.  Additional  support  from  the  industrial  sponsors  of  the 
Center  for  the  Engineering  of  Electronic  and  Acoustic  Mate¬ 
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Abstract.  The  eigenmodes  of  a  spherical  resonator,  with  a  perfectly  conducting 
wall  and  filled  with  a  homogeneous  isotropic  chiral  medium,  have  been  identified 
and  studied.  In  particular,  the  roots  of  the  dispersion  equation  have  been  given  as 
functions  of  /3/a.  where  a  is  the  cavity  radius  and  ft  is  the  chirality  parameter. 
Application  as  a  microwave  circuit  element  is  suggested. 


1.  Introduction  and  preliminaries 

Though  the  phenomenon  of  chirality  is  known  chiefly 
at  the  molecular  level,  and  therefore,  at  frequencies  in 
or  above  the  ultraviolet  range,  it  has  been  suggested 
(Lakhtakia  et  al  1988)  that  particles  endowed  with  chi¬ 
rality  can  exist  at  even  lower  frequencies,  say,  in  the 
GHz  range.  This  is  because  chirality,  or  handedness,  is 
a  geometric  property;  for  example,  the  electromagnetic 
(em)  response  of  a  right-handed  helix  is  different  from 
that  of  a  left-handed  one  (Varadan  et  al  1988).  Fur¬ 
thermore,  by  embedding  such  chiral  particles  in  a  low- 
loss  dielectric  medium,  the  resulting  composite  medium 
too  will  possess  handedness.  With  advances  in  polymer 
science,  it  is  becoming  increasingly  possible  that  such 
artificial  materials  can  be  manufactured  with  ease,  and 
their  properties  tailored  by  altering  the  sizes  and  con¬ 
centration  of  the  embedded  chiral  particles. 

Significant  advances  have  taken  place  recently  in 
the  formulation  of  a  frequency-domain  electromagnetic 
theory  for  chiral  media;  we  have  summarised  these 
elsewhere  (Lakhtakia  et  al  1988).  With  these  devel¬ 
opments,  it  is  time  that  aspects  relating  to  the  appli¬ 
cation  of  chiral  media  for  practical  problems  be 
explored.  To  that  end,  we  consider  here  the  eigen¬ 
modes  of  a  perfectly  conducting  sphere  filled  with  a 
homogeneous,  isotropic,  chiral  medium.  Such  an 
arrangement  constitutes  a  microwave  resonator,  an 
important  circuit  element  (Harrington  1964). 

Consider  a  source-free  region  occupied  by  an  iso¬ 
tropic  chiral  medium  in  which  the  usual  constitutive 
relations  D  =  eE  and  B  =  fxH  are  not  adequate  because 
of  their  incompatibility  with  the  handedness  of  the 
medium.  Instead,  the  relations 

D  =  eE  +  fieV  x  E  B  =  fiH  +  ftiV  x  H  (1) 

hold,  and  satisfy  the  requirements  of  time-reversal  sym- 
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metry  and  reciprocity;  here  /?  is  the  chirality  parameter 
measured  in  units  of  length.  Following  Bohren  (1974), 
the  EM  field  is  transformed  to 


E  =  Q  i  +aRQ  2 

H  =  Q2  +aLQ{ 

(2) 

where  the  left-  and  the  right-circularly  polarised  (lcp 
and  rcp)  fields,  Q{  and  Q2,  respectively,  must  satisfy 
the  Helmholtz  equations 

(V2  +  Yt)Qi  —  0 

(V2  +  H)0:  =  0 

(3) 

along  with  the  rotational  conditionst 

v  x  Q\  -  YiQi 

V  x  Qz  -  ~y2Qi- 

(4) 

Needless  to  say,  these  fields 
i.e. 

are  also  divergence-free 

V-Ci«0 

VQ2  =0. 

(5) 

In  these  equations,  the  two  wavenumbers  are  given  by 

$ 

1 

r—4 

II 

£ 

Y:  =  k/(  1  +  m 

(6) 

and 

=  -i(£/>)1/2 

aR  =  ~i  (ti/e),/2. 

(7) 

An  exp(-iart)  time  dependence  has  been  assumed 
throughout  this  work,  while  k  =  a)( iie)i/2  is  simply  a 
short-hand  notation. 


2.  Fields  inside  the  spherical  resonator 

Consider  now  the  spherical  resonator  r  =  a  which  is 
bounded  by  perfectly  conducting  walls  and  is  wholly 
filled  with  a  chiral  medium.  Using  a  representation 
due  to  Bohren  (1974)  and  justified  by  us  elsewhere 

t  We  are  indebted  to  an  anonymous  referee  for  pointing  out  that 
Qi  and  Q2  are  Beltrami  fields  (Caraman  1974). 
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Figure  1 .  Roots  (x  =  ka)  of  the  dispersion  equation 
F(x;  n)  =  0  as  functions  of  fi/a  for  n-  1.  Computations  for 
figures  1-4  were  performed  on  a  MAC  II  minicomputer 
using  Absoft  macfortran/020. 


Logip/a) 

Figure  2.  Roots  (x  =  ka)  of  the  dispersion  equation 
F(x;n)  =  0  as  functions  of  p/a  for  n  =  2. 


(Lakhtakia  et  al  1987),  the  field  inside  the  cavity  can 
be  adequately  expressed  by 

Q\(r)=  2  AtL(vl)(Yir) 

jm/i 

fi2(r)=  2  ByR^(y2r). 

V 

In  these  equations,  the  functions 

U,)(yir)  =  M(J\Yir)  +  Ni'KYir) 

R^{Yir)=M^i.Yir)-NW{Y2r) 

are  left-  and  right-circularly  polarised,  respectively. 
The  subscript  v  is  a  triple  index — the  index  s  can  be 
even  or  odd,  the  index  n  varies  from  1  to  and  the 
index  m  assumes  values  from  0  to  n;  the  well  known 
vector  spherical  wavefunctions  M{^(ar)  and  N^\or) 


(8) 

(9a) 

(9b) 


U. 


3  o  A 


20  21 


x  =  ka 


Figure  3.  Function  F(x;n  =  1)  for  /3/a  =  0  and  /3/a  =  0.04. 


have  been  defined  by  Morse  and  Feshbach  (1953), 
while  A„  and  B„  are  the  unknown  expansion  coef¬ 
ficients. 

The  determination  of  the  expansion  coefficients  can 
be  accomplished  using  the  boundary  conditions 
e,  x  E  =  0  on  the  surface  r  =  a,  where  E  can  be  syn¬ 
thesised  using  (2)  and  (8),  and  eT  is  the  unit  radial 
vector.  On  taking  the  inner  productt  of  Eg  with 
(m'i r?'(cos9)  sinm»  and  that  of  Ev  with 
(f£'(cos0)cosm»  on  r  =  a,  and  adding  the  two 
results,  it  is  easy  to  see  that 

Aemn)n(X\)  +  ttRBrmnJn(x2)  —  0.  (10a) 

In  these  equations  /„(or)  is  the  spherical  Bessel  function 
of  order  n ,  ;C(cos0)=  P?(cos  0)/sin  0, 
r”(cos  0)  =  dP™(cos  0)/d0,  P™(cos  0)  is  the  associ¬ 
ated  Legendre  function  of  order  n  and  degree  m,  and 
xp  =  YPa  f°r  p  —  1,2.  Likewise,  after  taking  the  inner 
product  of  Eg  with  (m'ji%'( cos  0)  cosm'tp)  and  that  of 
Ey  with  (-r™  (cos  0)  sin  m'<p)  on  r  =  a,  the  equation 

A  omn  jn  (x  1 )  +  OrB  omnj  n(x2)  —  0  ( 10Z?) 

can  be  obtained.  Furthermore,  from  the  inner  products 
of  Ee  with  (r£(cos0)  cos  m'<p)  and  of  -Ev  with 
(m'jr™  (cos  0)  sin m'q>)  on  r  =  a,  one  can  obtain 

Aemnx2d%pn(xl)  -  aRBemHxxdyn(x2)  =  0  (10c) 

in  which  rpn(%)  =  | /„(§)  is  the  Riccatti-Bessel  function 
and  dV>„(£)  *  dt//„(§)/d§.  Finally,  adding  the  inner 
product  of  Ee  with  (r™'(cos  9)  sin  m'<p)  to  that  of  Ev 
with  ( m'jt "  (cos  9)  cos  m'cp)  on  r  =  a  yields 

Aomnx2dxpn(x i)  -  aRBomnxxd y„(x2)  =  0.  (lOd) 
Is  it  obvious  from  (10)  that  azimuthal  parity  in  terms 
of  the  indices  s  and  m  is  conserved  for  the  present 
problem,  which  contributes  to  the  degeneracy  of  the 
eigenvalues. 

t  The  inner  product  of  two  functions  f(9,  <p)  and  g(9,  <p),  as  used 
here,  is  defined  as  the  integral 

f  d<P 
J0 


I  dtfsin  Of (9,  <p)g(9,  <p). 
Jo 
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>wH 


x  =  ka 

Figure  4.  Function  F(x;  n  =  2)  for  /3/a  =  0  and  /3/a  =  0.04. 

The  dispersion  equation  for  the  chiral-filled  res¬ 
onator  is  simply  obtained  by  the  simultaneous  solution 
of  (10a)-(10d);  consequently,  each  solution  of  the 
equation 

F(x;  n )  =  H>n(x i)  drpn(x2)  +  d*pn(xx )ip„(x2)  =  0 

x  -  ka  (11a) 

adds  yet  another  index  q  (for  its  gth  zero),  represented 
jointly  by  xnq,  x,„q  and  x2m}.  Corresponding  to  this  qth 
solution  of  (lla),  the  ratio 

^smnq/^R^smnq  ~  '  j  n(.X  2nq)  / j n(.X  Inq)  •  0  ^ ) 

Thus,  mode-ordering  for  the  chiral-filled  spherical 
resonator  requires  four  separate  indices,  of  which  s,  m 
and  n  are  as  before,  and  q  denotes  the  <jth  zero  of 
(11a)  for  a  given  n.  The  eigenmodes  of  this  cavity  are 

Esmnq  —  jn  i*2nq  v  '  (V  lnqr)  "b  jn  (-*•  Inq  )R  v  (  Y  2nq  fj 

(12) 

aR  Hsmnq  =  jnix^L^iYlnqr)  +  jniXlnqW^iYTnqr) 

(13) 

and  the  field  inside  the  chiral-filled  resonator  can  be 
expressed  in  terms  of  these  modes: 

E  =  2  FsmnqEsmnq  H  =  2  FsmnqHsmnq  (14) 

smnq  smnq 

with  the  usual  restriction  that  the  m  =  0,  s  =  odd  mode 
is  non-existent.  In  these  equations  xpnq  =  ypnqa  for  p  = 
1,  2;  and  (o^  =  (ep)l/2/k„q  are  the  resonance  fre¬ 
quencies,  with 

x ^  =  (k^a)-1  =  i(*r4  +  x^l,).  (15) 

The  stored  energy  in  a  given  eigenmode  can  be 
computed  in  either  of  two  ways  (Lakhtakia  et  al  1988): 
using  the  relation 

W}mnq  =  JR c(Esmnq  •  D*mna  -  Hsmnq  •  Btmnq)  (16a) 


or  equivalently  through  the  expression 

^  smnq  =  JRc(f  * D smnq  *  D smnq  ~  [A  1 B lmnq  •  Bsmnq). 

(16  6) 

In  either  case,  it  can  be  shown  that 

^ smnq  ~j n{x  inq)j  n{x2nq)^s{y  [„qY  2nq  /  k^) 

x  f  dq>  f  d  0  sin  0 

•'o  ■'o 

*  f  drr2L<»n  (ylnqr)  ■  Bi^ty^r)  (lla) 

provided  e,  \i  and  /?  are  all  assumed  real.  The  angular 
integrals  can  be  evaluated  analytically  to  yield 

^ smnq  J n  Xnq  )/n  (*^2 nq  )^£(  /  Xnq  Y 2nqf^nq  ) 

x  [(2 n  +  l)(n  -  m)!]_1[2jr(l  +  6^) 
x  n(n  +  l)(n  +  m)!]/^  (176) 

in  which  <5W  is  the  Dirac  delta,  and  the  integral 

hq  =  f  dr  r2{jn(Ylnqr)jn(Y2nqr) 

X  [1  -n(n+  1  )/Y\nqY2nqr2] 

-  dV’n(Yinqr)drl>n(Yinqr)/YlnqY2„qr2}  (17 c) 

may  have  to  be  evaluated  numerically  due  to  its  signifi¬ 
cant  complexity. 

It  is  well  known  that  although  lcp  and  rcp  waves 
can  propagate  independently  in  an  unbounded  chiral 
region,  at  a  bimaterial  interface  mode-conversion 
occurs  (Lakhtakia  et  al  1985).  Therefore,  pure  LCP  or 
RCP  modes  cannot  exist  within  a  bound  chiral  volume, 
as  is  also  apparent  from  (12)  and  (13).  Moreover,  in 
view  of  the  definitions  (6),  the  solutions  x^  of  (lla) 
do  not  depend  on  the  sign  of  the  chirality  parameter 
>3;  if  /3  changes  its  sign  then  xlnq  and  x^  are  simply 
interchanged.  In  the  sequel  it  is  assumed  that  /3  is  real 
and  positive. 

3.  Results  and  discussion 

Returning  to  (lla),  it  is  clear  that  with  /3/a  =  0,  i.e., 
when  the  medium  is  non-chiral,  this  equation  reduces 
to 

Vn(x)dipn(x)  =  0  x  =  ka  (18a) 

in  which  case  the  TE-to-r  eigenvalues  come  from  the 
zeros  of  the  equation  (Harrington  1964) 

jn(x)  =  0  x  =  ka  (186) 

and  yield  the  transverse  electric  modes  with  er  •  E  =  0. 
On  the  other  hand,  the  TM-to-r  eigenvalues  are  the 
solutions  of  the  equation  (Harrington  1964) 

dy>n(x)  =  0  x  =  ka  (18c) 

and  give  the  transverse  magnetic  modes  with  er  •  H  = 
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0.  Furthermore,  when  xfi/a  <  1,  then  the  Riccatti- 
Bessel  functions  in  (11a)  can  be  expanded  around  x  to 
first  order  in  x2fi/a  using  a  Taylor  series  expansion; 
with  this  approximation  (11a)  consequently  reduces  to 

y>n(x)dy>„(x)  -  (x2fi/a)2  dipn(x)  d2V„(x)  =  0. 

(19a) 

But 

x2d2' l>n(x)  =  [n(n  +  1)  -x2]ipn(x) 
so  that  (19a)  further  simplifies  to 
V>H(x)  d \l>n (x)[l  -  x2fi2a~2n(n  +  1)  +  x*fi2a~2]  =  0 

x  =  ka  (19b) 

suggesting  that  the  eigenvalues  for  xfi/a  <  1  do  not 
appreciably  differ  from  the  case  /3/a  =  0.  The  third 
factor  on  the  right  side  of  (196),  however,  suggests 
additional  roots  (due  to  chirality)  which  will  be  signifi¬ 
cant  when  xfi/a  is  larger. 

The  roots  (*)  of  (11a)  are  plotted  in  figures  1  and 
2  for  n  =  1  and  n  =  2,  respectively,  as  functions  of  the 
parameter  fi/a,  and  these  figures  confirm  the  deduc¬ 
tions  of  the  previous  paragraph.  Only  those  roots  have 
been  given  for  which  both  xXnq  and  are  non-nega¬ 
tive.  There  is  considerable  deviation  from  the  achiral 
values,  however,  when  fi/a  becomes  larger  than  about 
10~2.  In  fact,  the  roots  are  extremely  crowded  around 
fi/a  =  10~3/2.  In  order  to  bring  out  this  crowding 
further,  shown  in  figures  3  and  4,  respectively,  are  the 
plots  of  F(x\  n)  as  functions  of  x  for  n  =  1  and  n  =  2; 
the  normalised  chirality  parameter  fi/a  =  0  and  fi/a  * 
0.04  in  these  two  diagrams.  The  increase  in  the  num¬ 
bers  of  zeros  of  F(x\  n )  for  fi/a  =  0.04  and  should  be 
noted.  As  xfi/a-*  1,  jcj  — *• 00  very  rapidly;  consequently 
the  cavity  becomes  electrically  large  for  the  lcp  com¬ 
ponents  (and  constant  but  small  for  the  rcp 
components),  creating  conditions  thereby  for  the 


accommodation  for  an  increasingly  larger  number  of 
eigenmodes. 

It  is  well  known  that  resonators  find  use  as  the 
microwave  analogues  of  lc  circuits.  For  fixed  n  and 
fi  =  0,  the  resonance  frequencies  are  spaced  far  apart, 
as  is  obvious  also  from  figures  1  and  2;  thus,  non-chiral 
resonators  are  essentially  narrow-band.  On  the  other 
hand,  when  fi/a  is  large  in  the  present  context,  the 
resonance  frequencies  are  closely  spaced.  Therefore, 
chiral  resonators  may  be  utilised  for  ‘wide-band’  pur¬ 
poses,  and,  in  particular,  to  replace  the  multicavity 
staggered  turning  arrangements  (Gewartowski  and 
Watson  1965)  in  klystron  amplifiers. 
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Abstract.  In  a  homogeneous  isotropic  chiral  medium,  the  electromagnetic  field 
has  to  be  described  in  terms  of  left-circularly  and  right-circularly  polarized  (LCP 
and  RCP)  components,  each  of  which  has  a  different  wavenumber  associated  with 
it.  In  general,  when  a  wave  propagating  in  a  homogeneous  isotropic  chiral 
medium  encounters  a  homogeneous  isotropic  chiral  scatterer,  both  the  scattered 
and  the  internally  induced  fields  contain  LCP  and  RCP  components  regardless  of 
the  state  of  polarization  of  the  incident  field.  It  is  proved  here  that  the  scattered  as 
well  as  the  internal  fields  have  the  same  state  of  polarization  as  the  incident  field,  if 
the  scatterer  and  the  surrounding  medium  are  impedance-matched;  this  conclu¬ 
sion  holds  regardless  of  the  two  chiral  parameters  involved  as  well  as  the  geometry 
of  the  scatterer.  Two  examples  are  also  given. 


1.  Introduction  and  preliminaries 

Homogeneous  isotropic  chiral  media  are  circularly  birefringent  [1,2]:  the 
electromagnetic  field  in  such  a  medium  has  to  be  described  in  terms  of  left-circularly 
and  right-circularly  polarized  components,  each  of  which  has  a  different  wave- 
number  associated  with  it.  Either  of  the  two  components  is  capable  of  independent 
propagation  in  an  unbounded  homogeneous  region.  However,  when  either  a  left- 
circularly  polarized  (LCP)  or  a  right-circularly  polarized  (RCP)  wave  hits  an 
isotropic  obstacle,  then  the  scattered  field  is  composed  of  both  LCP  and  RCP 
components;  in  addition,  if  any  field  can  be  induced  inside  the  obstacle,  it  too  is 
composed  of  both  components. 

Chiral  media  can  come  in  two  forms,  one  being  the  mirror-conjugate  of  the  other; 
this  is  commonly  exemplified  by  L-type  and  D- type  enantiomorphs  of  many  organic 
molecules  [3],  A  homogeneous,  isotropic,  reciprocal  chiral  medium  can  be  described 
by  the  constitutive  equations 

D  =  e[E+/JV  x  E],  B=/z[H  +  0VxH],  (la,  b) 

in  which  e  and  /i  have  their  usual  meaning,  while  /?  is  the  chirality  parameter.  Two 
media  having  the  same  e,  n  and  |/J|,  but  with  their  chirality  parameters  differing  in 
sign,  are  mirror-conjugates  of  each  other.  In  a  recent  paper  [4],  it  was  shown  that 
when  a  plane  wave  is  scattered  by  the  planar  interface  between  two  mirror- 
conjugated  chiral  media,  the  reflected  as  well  as  the  refracted  plane  waves  are  of  the 
same  circular  polarization  state  as  the  incident  plane  wave.  This  phenomenon  will  be 
generalized  here,  and  shown  to  exist  at  the  interface  of  any  two  homogeneous 
isotropic  chiral  media,  provided  that  they  are  impedance-matched. 
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Following  the  pioneering  work  of  Bohren  [2],  the  electromagnetic  field  in  a  chiral 
medium  is  transformed  by  the  prescription  of  a  LCP  field  Qt  and  a  RCP  field  Q2> 
defined  as  per  the  following  relations: 

E=0j+aR02,  H  =  aL0j+02,  (2a,6) 

with 


It  should  be  noted  that  aR  has  the  units  of  an  impedance  and  aL>  that  of  an  admittance. 
In  an  unbounded  chiral  medium  LCP  and  RCP  fields  can  exist  independent  of  each 
other,  while  linearly  polarized  fields  are  not  allowed  [5].  From  the  aforementioned 
transformation,  it  follows  that  E  =  i(/i/e)1/2H  for  a  LCP  field,  and  E=  —  i(/i/e)1/2H  for 
a  RCP  field;  consequently,  the  ratio  (file)112  may  be  considered  as  the  intrinsic  wave 
impedance  of  a  chiral  medium,  just  as  in  the  case  of  achiral  media  [6]. 

Substitution  of  equations  (2  a,  6)  into  source-free  Maxwell’s  equations,  with  a 
harmonic  time-dependence  exp[  —  icot],  leads  to  the  circulation  equations 

VxQ^yjQi,  VxQ2=-y2Q2,  (4a,  b) 

both  Qi  and  Q2  are  solenoidal.  The  wavenumbers  corresponding  to  the  LCP  and  the 
RCP  fields  are  given  as 


Vl  1  —kp  ’  l+kfi’  (5) 

in  which  k  —  co(en)112  is  simply  a  shorthand  notation  and  does  not  represent  any 
wavenumber  except  when  fl—O. 


2.  Analysis 

In  an  unbounded  source-free  region  occupied  by  an  isotropic  chiral  medium,  it 
can  be  shown  that  E  and  H  satisfy  the  same  homogeneous  governing  differential 
equations, 


V2E  +  2yly2^VxE  +  yly2E=0,  V-E=0,  (6) 

V2H  +  2y1y2^VxH  +  y,y2H=0,  V-H  =  0.  (7) 

For  this  isotropic  chiral  medium,  the  infinite  medium  Green  function  ®(r,  r0)  satisfies 
the  dyadic  inhomogeneous  relation  [5] 

[VV  -  V23  -  y  ty23  -  2y  jy20V  x  3]  •  ®(r,  r0)  =  3  S(r— r0),  (8) 

3  being  the  unit  dyadic.  Its  solution  is  given  by  the  decomposition 

®(r,r0)  =  ®1(r,ro)  +  ®2(r,ro),  (9  a) 


Influence  of  impedance  mismatch 


1387 


s 

Quantities  of  relevance  to  the  scattering  analysis. 


with 


©i(r,r0)  = 

©2(r.ro)  = 


k 


k 

8nyiy2 


[yi3  +  yrlVV  +  Vx3]^(yi;|r-r0|), 
[y +  y 2  u W  - V  x  3]g(y  2;  |r  -  r J), 


(96) 
(9  c) 


while  g(o;  i?)  =  exp  [iaRjjR  is  a  scalar  Green’s  function. 

Consider  the  region  V+  shown  in  the  figure,  which  is  occupied  by  a  chiral 
medium  (e,p,P).  The  obstacle  volume  F_  is  occupied  by  another  chiral  medium 
(s',  p',  ft'),  no  restriction  being  placed  on  the  chirality  parameters  of  either  medium. 
Let  an  electromagnetic  field,  {Einc,  H,nc},  be  incident  on  the  obstacle,  it  being 
understood  that  the  source  of  the  incident  field  lies  outside  the  circumscribing  sphere 
of  V _ .  As  per  Huygens’s  principle  and  the  extinction  theorem  [5, 7]  applied  to  V+ ,  it 
turns  out  that 


—  Ein'(r)  =  (yJy2A  2)jJ  d2r0©1(r,r0)*[ico/t«1,x H+(r0)+*«nx E+(r0)] 

+  (yi72*"2)  |  j^d2r0©2(r,  r») '  x  H+(r.)-*eB  x  E+(r„)], 


II 


reF_,  (10a) 


—  Hinc(r)  =  (y1y2A~2)  d2^©^,  r0)  •  [*•„  x  H+(r0)-icoee0  x  E+(r0)] 


— (Vi72*~z)  J*  J  d2r0©2(r,  r„) •  [6ea  x  H+(r0)  +  ia)ee0  x  E+(r0)], 

reF_.  (106) 

In  these  equations,  •„  is  the  unit  outward  normal  to  the  bimaterial  interface  S\  and 
{E  +  .H+}  denote  the  electromagnetic  field  on  the  exterior  side  of  S.  The  Green 
functions  used  above  also  refer  to  the  medium  in  V+ . 


■  »  .V»-i  a  * 
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Equations  (10  a,  6)  form  the  bedrock  of  the  T’-matrix  method  [7, 8].  Once  these 
have  been  solved,  the  scattered  field,  {Fc,  H**},  can  be  easily  determined  via  the 
relations 

E“(0  =  (y iVi* " 2)  |  £  d2r0 © , (r, r0) •  [i£ojten xH+(r0)  +  *»„xEt (r„)] 

+  (yiV2*"2)  |  |s  d2r0®2(r, r„) •  [i<o/i«n  xHt(rJ-l»„  x  E+(r0)], 

reF+,  (11a) 

H’c(r)  =  (yIy2A“2)  J*  £  d2^®^!,  r0)-[£«n  x  H+(r0)— ia»een  x  E+(r0)] 

—  (Viyx*-2)  j*  |  d2r0®2(r,r0)'[^enxH  +  (r0)  +  ia;eenxE+(r0)], 

reF+.  (116) 

Let  the  transformation  (2  a,  6)  be  applied  to  both  the  incident  and  the  scattered 
fields.  As  a  consequence,  equations  (10a,  6)  and  (11  a,  6)  can  be  respectively 
rewritten  in  the  forms 


“>JI 


-Oi“c(r)  =  (y1y2*“l)  d2r0®1(r,r0)-[*nxE+(r0)-altenxH+(r0)],  reV_, 


Q2C(0  =  (y1y2*  d2r0®2(r,r0)-[enxH+(r0)-aLenxE+(r0)],  reF_, 


0“(r)=(y1y2* 


-1)  J  j*^d2r0®2(r 

4>  T  is  d2r°®^r’ 
■‘){{f  d2r0®2(r, 


r.)  •  [en  x  E+(r0)— aRen  x  H+(re)],  reV+, 


~02(r)=(yly2*  l)  d2r0®2(r,  r0)- [en  x  H+(r0)-aLen  x  E+(r0)],  reV+. 


The  boundary  conditions  on  the  interface  S’  are  the  continuity  of  the  tangential 
components  of  the  E  and  H  fields  across  S.  Consequently,  equations  (12  a,  6)  and 
(13  a,  6)  can  be  re-expressed  to  bring  in  the  field,  {Pn‘,  H1"1},  induced  inside  F_  [5, 8] 
as 

—  Qirc(r)  =  (y1y2A_1)  J  £  d2r0  S^r,  r„)  •  [•„  x  Ein,(r0) -  aRen  x  H,m(r0)],  reV., 


-n 


Qz‘(r)=(yiyz*  )  d2r0 ®2(r ,  O •  [•„  x  Hio,(r0)  —  aLe„  x  Eint(r0)],  «=F_, 


046) 
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Or(r)  =  (yl72A 


l)|  £d2r.®i(r. 


r„)  *  [•„  x  ^"'(r,,)  —  a*«n  x  Hin,(r0)],  ref\, 


-Q\%T)  =  (yiy2k 


■>I1 


d2r0ffi2(r,  r0)  •  [en  x  Hin,(ro)-aLe0  x  Eim(r0)],  reV+. 


But,  Bohren’s  decomposition  (2  a,  b)  also  applies  to  the  medium  inside  V _ ,  albeit 
with  the  primed  constitutive  parameters.  In  strict  analogy  with  expression  (3),  the 
quantities 

1 W’-  ?  "6I 

can  also  be  defined.  Consequently,  the  foregoing  sets  of  integral  equations  are 
transformed  to 

-Qi“c(r)=7iy2A-1[l  +  aR/a'R ]  J  £  d’r,®^,  rj-  [•„  x  Qf  (r0)] 

+yiy2*“lK-aR]||  d2r0©1(r,r0)-[onxQi2,,,(r0)],  r eF_,  (17 a) 


-II 


Q’c(t)  =  7i72*  HI+Or/or]  d2r0©l(r,r0)-[enxOin,(r0)] 


+yiy2* 


d2r0©1(r,r0)-[enxQ‘2«(r)],  reF+,  (176) 


— 02l*(r) = y  jy2*  ” 1  t°L — J  Jsd2r0©2(r,r0)-  [e„  x  0^(r0)] 

-yiy2*_l[l+aR/aR]|  |  d2r0 ©2(r,  r0)  •  [eB  x  Q2m(r0)],  reF_, 


0‘/(r)  =  yty2k~  l[aL  -  a|J  J  £  d\ ®2(r,  r0) •  [e„  x  Q‘»‘(r0)] 

-y^'Hl+aR/aRjJ  |  d2r0©2(r,  r0)  * [•„  x  Q2M(r0)],  reV+. 


As  was  expected,  equations  ( 1 7  a,  b)  and  (1 8  a,  6)  show  that  both  the  scattered  and 
the  internally  induced  fields  contain  LCP  and  RCP  components  regardless  of  the 
state  of  polarization  of  the  incident  field.  It  should  be  noted  that  in  deriving 
equations  (17  a,  6)  and  (18  a,  6),  the  chiral  parameters  of  either  medium  were  not 
explicitly  used;  nor  do  equations  (2  a,  6)  need  them.  In  fact,  these  equations  are 
applicable  whether  both,  any  one  or  neither  of  the  two  materials  are  chiral. 
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0‘ic(r)  =  (Vly2A-1)  J  £ 
-0‘/(D  =  (yly^-1)  J  Js 


d2r0©i(r,  r0)  •  [en  x  E‘B,(r0)  - aRe„  x  Hio'(r#)], 


d2r0  ©2(r,  r0)  •  [en  x  Hin*(r0)  -  aLen  x  Eint(r0)], 


reV+, 

(IS  a) 

reV+. 

(15  b) 


But,  Bohren’s  decomposition  (2  a,  6)  also  applies  to  the  medium  inside  V _ ,  albeit 
with  the  primed  constitutive  parameters.  In  strict  analogy  with  expression  (3),  the 
quantities 


(16) 


can  also  be  defined.  Consequently,  the  foregoing  sets  of  integral  equations  are 
transformed  to 


-II 


_Or(r)=y1y2*-1[l+aR/4]  d2r0©1(r,r0)-[enxQHr0)] 


+Viy2* 


-II 


d2r0©1(r,r0)-[enxQ,2Dl(r0)],  reF_,  (17  a) 


Qic(r)  =  V172*  ‘[l+aR/aU  d2r0©1(r,r0)-[enxQy’t(r0)] 


+  yi72* 


‘K-«j||s 


d!r.®,(r,r„)-[.„xQi”(r)l,  reV,,  <174) 


~OT(r)=yIy2A 


||. 


d2r„  ©2(r,  r0)  •  [•„  x  Qin,(r0)] 


-Viyik-'V+a'M  J  |s  d2r0©2(r,r0)-[« 


•nXQnO],  reV., 


(18  a) 


Oi(r)=yly2k 


"K-4U|[ 


d2r0  ®2(r,  r0)  •  [en  x  Q\n‘(r0)] 


-yiy2*'l[l  +  «»/««]  J  Js  d2r0 ®2(r, r0) •  [en x Q‘2  *(r0)],  reF+. 


(18  b) 


As  was  expected,  equations  (17  a,  b)  and  (18  a,  b)  show  that  both  the  scattered  and 
the  internally  induced  fields  contain  LCP  and  RCP  components  regardless  of  the 
state  of  polarization  of  the  incident  field.  It  should  be  noted  that  in  deriving 
equations  (17  a,  b)  and  (18  a,  6),  the  chiral  parameters  of  either  medium  were  not 
explicitly  used;  nor  do  equations  (2  a,b)  need  them.  In  fact,  these  equations  are 
applicable  whether  both,  any  one  or  neither  of  the  two  materials  are  chiral. 


1390 


A.  Lakhtakia  et  al. 


One  notable  case  stands  out:  let  the  scatterer  be  impedance-matched  with  the 
surrounding  medium;  fi'le'  =  nle.  In  that  case,  equations  (17  a,  6)  and  (18  a,  6) 
simplify  to 

—  Oinc(r)  =  2y,y26-1 j  jsd2r0©1(r,r0)-[enxQ^(r<>)], 

reF_, 

(19  a) 

Q‘ie(r)  =  2yly26‘ 1  j  d2r0 ©t(r,  r0)  •  [en  x  Q^rJ], 

reF+, 

(196) 

Qi2"e(r)  =  2y1y2ft-1  J  £  d2r9S2(r,  rB)  •  [en  x  Qf(r#)], 

reF_, 

(20  a) 

-Q*2‘(r)  =  2y1y26-1 1  £  d2r0®2(r,ro)-[e,xQ“'(r0)], 

reF+. 

(206) 

Equations  (19  a,  6)  and  (20  a,  6)  clearly  show  that  the  scattered  as  well  as  the  internal 
fields  have  the  same  state  of  polarization  as  the  incident  held,  if  the  scatterer  and  the 
surrounding  medium  are  impedance-matched.  This  conclusion  is  not  affected  by  the 
value  of  the  chirality  parameter  in  either  of  the  two  media,  and  forms  the  major  result 
of  this  paper. 

3.  Scattering  by  planar  interface 

As  an  example  of  the  analysis  given  above,  let  the  reflection  and  refraction  of 
plane  waves  at  bimaterial  interfaces  be  considered.  Let  F+  be  the  half  space  z^O; 
and  F_,  the  half  space  z>0.  The  interface  S  is  the  plane  z  =  0. 

In  V+ ,  the  unprimed  constitutive  parameters  are  used,  and  the  incident  fields  can 


be  represented  by  [9] 

Q\“c(r)  =  4 ,  C®,  +  K  -  “i®* + »cez)/y1]  exp  [i(icx  +  o^z)],  reV+,  (21  a) 

Qife{r)  =  A2[9,+i(.<i1Bx-Kez)ly2]exp[i(Kx  +  <x2z)],  reF+,  (216) 

while  the  scattered  (i.e.  reflected)  fields  are  given  by 

01*(r)=sf?1[*j»+i(«iex+Kex)/yl]exp  (i(icx— atz)],  reF+,  (22a) 

Q?(r)=R2[By-i(a2ex+Kaz)ly2]exp[i(Kx-tz2z),  reV+.  (22  b) 

The  fields  induced  in  V_  (i.e.  the  refracted  fields)  are  to  be  given  as 

Q'1m(r)  =  r1[e,-(-i(— a'1eI+x»I)//1]exp[i(K*+a',z)],  reF_,  (23a) 

Qj’(r)  =  T2[»y+ i(a'2«x — »cer)//2]  exp  [i(ic* + a^z)],  reF_,  (236) 


in  which  the  primed  constitutive  parameters  have  been  used.  In  these  equations,  k  is 
the  (common)  horizontal  wavenumber  required  to  satisfy  Snell’s  laws  at  the  interface 
5;  •„ . . .  are  the  unit  Cartesian  vectors;  while  =  +  (y2  —  x2)1/2,  a2  =*  +  (y§  —  k2)1/2, 
aj  =  -(y'i—K2)112  and  a'2  =  +(y'22-K:2)1/2. 

The  x-directed  and  y-directed  components  of  the  E  and  the  H  fields  have 
to  be  continuous  across  S ;  in  addition,  the  two  media  are  impedance-matched, 
namely  n’l& =/i/e.  Use  of  these  conditions  leads  to  the  solutions 

=  +  1Al,  Ri^Tx  —  Au  (24) 

T2  =  2y’2a2\y'2<i2+a!2y^\  iA2,  R2=T2  —  A2.  (25) 
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This  illustrates  the  observation  that  the  incident,  the  reflected,  and  the  refracted 
plane  waves  are  of  the  same  circular-polarization  state  whenever  the  two  chiral  media 
on  either  side  of  a  planar  interface  are  impedance-matched.  Given  that  fi'/e'  =  pl&,  it 
should  also  be  noted  from  equations  (24)  and  (25)  that  Ti/A1  =  T2/A2  =  1  for  the  case 
of  normal  incidence  (rc  =  0);  further,  if  j3  =  0,  then  the  case  of  normal  incidence 
trivially  satisfies  the  Brewster  condition  [10]. 


4.  Scattering  by  a  sphere 

The  second,  and  the  last,  example  involves  the  scattering  characteristics  of  a 
sphere.  Let  V _  be  a  sphere  of  radius  a.  The  incident  field  can  be  expanded  as 

Or«=  £  AtM^O  +  Nt^r)],  (26 a) 

v~smn 

Qf(r)=  £  5v[M<l>(y2r)-N«l>(y2r)].  (26  6) 

v*5m* 

The  subscript  v  is  a  triple-index — the  index  s  can  be  even  or  odd,  the  index  n  varies 
from  1  to  X),  and  the  index  m  assumes  values  from  0  to  n;  the  well  known  vector 
spherical  wavefunctions,  M^'fcr)  and  H^\ar,j=  1, 3,  have  been  defined  by  Morse 
and  Feshbach  [11],  while  A „  and  B,  are  the  expansion  coefficients. 

The  scattered  field  is  likewise  expressed  in  terms  of  the  vector  wavefunctions  as 


Qic(r) = Z  Fy[  M<3>(y ,  r)  +  N‘3>(y  l  r)] , 

V 

reV+, 

(27  a) 

Q?(r) = I  Gy[M<3>(y2r)  -  N<3>(y2r)], 

V 

reF+, 

(276) 

while  the  field  induced  inside  F_  has  the  representation 

Q‘r(r)=Z  C¥[M‘I)(y'ir)  +  N^r)], 

reF_, 

(28  a) 

Q‘2m(r) = Z  ^.[W<.,,(y'2r)  -  N["(y'2r)], 

reF_. 

(286) 

After  enforcing  the  continuity  of  the  tangential  components  of  the  E  and  the  H 
fields  across  the  sphere  r=a,  utilizing  the  orthogonalities  of  the  spherical  harmonics 
[11],  and  recognizing  that  the  two  media  are  impedance-matched  (i.e.  p'le'  =  ple),  one 
obtains  the  following  solution: 

jFv  =  -  Ay[cil/n(y1a)\l/a(y\a)  -  ^.(y  ia)#„(/ia)][5<5n(y  ,a)^n(y',a)  -  ^n(yia)d\f/n(y  ',a)] " 1 , 

(29  a) 

C,  =  |0[^B(y  i^.(y'.a)  -  <J„(yla)#n(/1a)]  - 1 ,  (29  6) 

G,  =  -  B,[d<l/a(y2a)>l/B(y2a)  -  ^B(y  2a)3^a(y'2a)][d^B(y2o)^n(y'2a)  -  ZB(y2a)d\l/B(y'2a)] " 1 , 

(30  a) 

£>,=B^^[^n(y2a)^(/2a)-^n(y2a)#0(y'2a)]-1. 


(306) 
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In  these  expressions  <l/n(x)  =  xja(x)  and  £n(x)  =  xhif\x)  are  the  Riccati-Bessel 
functions  [12];  and  di^n(x)  =  di^n(jc)/djc,  and  so  on.  Again,  the  major  result  of  this 
paper  has  been  exemplified  by  equations  (29  a,  b)  and  (30  a,b). 
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ABSTRACT 

We  employ  the  multiple  scattering  formalism  to  obtain  two  independent  expressions  for  the  effective 
permittivity  <  e>  and  permeability  <  p  >  of  microwave  composite  materials.  These  two  expressions  c'e 
required,  but  not  provided  in  the  literature,  in  the  investigation  of  the  reflection  and  transmission 
characteristics  for  a  layer  of  microwave  composite  materials  in  response  to  incident  millimeter-waves  and 
microwaves.  In  order  to  check  the  validity  of  the  expressions,  several  cases  involving  lossy/ loss’  ess  matrix 
and  inclusion  materials  are  discussed. 

1 .  INTRODUCTION 

It  is  well  known  that  when  waves  propagate  through  a  medium  containing  scatterers,  the  entrained 
energy  (intensity)  is  going  to  be  either  redistributed  in  various  directions  by  scattering  or  absorbed  by 
intrinsic  absorption  mechanisms  This  suggests  the  well  known  idea  that  to  achieve  a  high  attenuation  of  the 
incident  wave  energy,  the  geometry  and  material  properties  of  the  scatterers  as  well  as  physical  mechanisms 
for  scattering  and  absorption  should  be  carefully  selected.  The  selection  is  not  an  easy  task  and  it  depends  on 
different  applications,  civilian  as  well  as  military  considerations  and  advancements  in  material  technology. 

In  many  applications  involving  microwaves,  it  is  desirable  to  design  materials  that  will  have  prescribed 
reflection  and  transmission  characteristics  as  a  function  of  frequency  (narrow  or  wide  band)  and  at  the  same 
time  conform  to  restrictions  on  weight,  structural  properties,  thickness,  etc.  Microwave  composite  materials 
that  contain  a  distribution  of  inclusions  of  specific  concentration,  distribution,  geometry  and  material 
properties  can  often  achieve  such  goals.  Since  the  actual  response  of  microwave  composite  materials  to  the 
incident  wave  of  low  to  high  gigahertz  frequency  is  quite  diverse,  it  becomes  very  expensive  and  time 
consuming  to  actually  prepare  samples  of  such  materials  and  test  them  experimentally.  However,  an  optimal 
design  through  theoretical  analysis  of  such  materials  is  relatively  efficient  and  the  parameters  involved  are 
tractable. 

When  particles  are  dispersed  in  a  host  medium  to  form  a  composite,  depending  on  the  volume  fraction  of 
particles,  single  or  multiple  scattering  dominates  the  scattered  energy  when  waves  are  incident  on  them. 
Multiple  scattering  effect  cannot  be  ignored  when  the  concentration  of  particles  is  considerable  as  in  the  case 
of  microwave  composites.  In  this  theoretical  investigation,  we  employ  the  multiple  scattering  formalism, 
which  has  been  previously  applied  to  get  the  effective  propagation  constant  for  the  ferrite  composites  [1],  to 
obtain  two  independent  expressions  for  the  effective  permittivity  <  e  >  and  permeability  <  p  >  of  microwave 
composite  materials.  These  two  expressions  are  required,  but  not  provided  in  the  literature,  in  the 
investigation  of  the  reflection  and  transmission  characteristics  for  a  layer  of  microwave  composite  materials 
in  response  to  incident  millimeter- waves  and  microwaves.  In  order  to  check  the  validity  of  the  expressions, 
several  cases  involving  lossy/lossless  matrix  and  inclusion  materials  are  discussed. 

2 .  MULTIPLE  SCATTERING  FORMALISM 

Consider  the  propagation  of  a  plane  harmonic  electromagnetic  wave  in  a  medium  referred  to  as  the 
matrix  characterized  by  complex  values  of  the  dielectric  function  £2  =  £2+1  e2" an(*  magnetic  permeability 
P2  *  M-2  +  j  M-2"-  Embedded  in  the  matrix  is  a  random  distribution  of  uniform  size  spherical  scatterers  whose 
permittivity  -  ej’  + j  ej"  and  permeability  pj  =  pj’  +  j  Pj”  .  The  matrix  of  volume  V  contains  N 
scatterers  such  that  N  -»=*>,  V  -»  <*>  but  n0  (=  N/V)  the  number  density  is  finite. 
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Let  E,  E°,  Eje  and  Ejs  be  respectively  the  total  electromagnetic  field,  the  incident  wave  field,  the  field 
exciting  the  ith  scatterer  and  the  field  scattered  by  the  ith  scatterer.  All  the  fields  have  implicit  exp(-icot) 
time  dependence  and  satisfy  vector  Helmholtz  equation.  Let  Re  <)>n  and  Ou  4>n  denote  the  basis  of  orthogonal 
functions  which  are  eigenfunctions  of  the  vector  Helmholtz  equation.  The  qualifiers  Re  and  Ou  denote 
functions  which  are  regular  at  the  origin  and  outgoing  at  infinity  [2].  Thus,  we  can  write  the  following  set 
of  self-consistent  equations: 

E=E°  +  fEf  (1) 

i=l 


E°  =  £  a1  Re  0  (r  -  r.) 

n  n  Tnv  v 

(2) 

Ee=£  a‘ Re<j)  (r-r.)  ;  a<lr-r.l<2a 

i  n  n  Tn'  i/  *  i 

(3) 

ES  =  2  £  Ou(j)  (r-r.)  ;  lr-r.l>a 

i  n  n  Tnv  v  i 

(4) 

where  anl  and  fn*  are  unknown  expansion  coefficients  and  an‘  are  known  incident  field  coefficients.  We 
observe  in  (3)  and  (4)  that  a  is  the  radius  of  the  sphere  and  that  all  expansions  are  with  respect  to  a  coordinate 
origin  located  in  a  particular  scatterer. 

The  T-matrix  by  definition  simply  relates  the  expansion  coefficients  of  Eje  and  Ejs  provided  the  sum  of 
them  is  the  total  field.  Thus  [2], 


and  the  following  addition  theorem  for  the  basis  functions  is  invoked 


(5) 


°u  <t>n(r  -  r.)  =  £  c  (r.  -  r )  Re  *  (r  -  r )  .  (6) 

11  J  n  nn  1  J  n  * 

Substituting  (2)  -  (6)  in  (1),  and  using  the  orthogonality  of  the  basis  functions  we  obtain 

N 

a  =  a1  +  X  T*  °(ri  -  rj)  aJ  .  (7) 

j*» 

This  is  a  set  of  coupled  algebraic  equations  for  the  exciting  field  coefficients  which  can  be  iterated  and  leads 
to  a  multiple  scattering  series. 

For  randomly  distributed  scatterers,  an  ensemble  average  can  be  performed  on  (7)  leading  to 

<a>i  =  a1  +  <  a  (r  -  r.)  T*  <oci>.j> .  (8) 

where  the  angle  brackets  and  their  subscripts  i  and  ij  denote  conditional  averages.  Equation  (8)  when  iterated 
is  an  infinite  hierarchy  and  usually  the  hierarchy  is  truncated  by  use  of  the  quasi-crystalline  approximation 
(QCA)  [3].  The  QCA  states  as 

<  cc*>..  =  <  a*>.  .  (9) 

ij  j 

Then,  (8)  simplifies  to 

<  a‘>.  =  j  +  <  a  (r.  -  r.)  T*  <o^>j>i  (10) 

an  integral  equation  for  <  a1  >j  which  in  principle  can  be  solved.  We  ob serve  that  the  ensemble  average  in 
(10)  only  requires  P(rjlrj),  the  joint  probability  distribution  function.  In  particular,  the  homogeneous 
solution  of  (10)  leads  to  a  dispersion  equation  for  the  effective  medium  in  the  quasi-crystalline 
approximation.  Defining  the  spatial  Fourier  transform  of  <  a1  >j  as 

<a>.  =  I  e  X  (K)dK  (11) 
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and  substituting  in  (10),  we  obtain  for  the  homogeneous  solution 


X‘  (K)  =  X  f  <J(r.  -  r)  T*  pfclr.)  x  '  'j)  Xj  (K)  dr 


If  the  scattered  arc  identical 


Xi(K)  =  XJ(K)=X(K) 

and  thus  for  a  nontrivial  solution  to  <  a1  >j,  we  require 

I  £  f  j  >K-(r ■  -  r-) 

T  \  I  ^  tJ  _  I.  \  „  1  J 


I~X  |o(r.-rj)Tip(rjlr.)i 


dr.  |  =0 


In  (14),  P(rjlrj)  is  the  joint  probability  distribution  function  and  I  is  an  identity  matrix.  For  isotropic  or 
spherical  statistics, 

(  0 ;  I  r.  -  r.  |  <2a 


P(rjlriM  ,  ,  ,  ,  (15) 

J  ^ g ( |  r. - r. | )/V  ;  |r.-r.|>2a 

where  we  have  assumed  that  the  scatterers  are  impenetrable  with  a  minimum  separation  between  the  centers, 
and  in  (15),  2a  could  be  the  diameter  of  the  circumscribing  sphere.  For  spherical  scatterers,  the  joint 
probability  distribution  depends  only  on  the  interparticle  distance  and  not  on  the  orientation  of  the  vector 
joining  the  centers;  and  the  function  g(l  rj  -  rj  I)  is  called  the  radial  distribution  function. 

3.  DISPERSION  EQUATION  IN  THE  LONG  WAVELENGTH  LIMIT 

The  objective  of  the  self-consistent  multiple  scattering  formalism  derived  for  the  problem  is  to  obtain  an 
analytical  expression  for  the  dispersion  equation  (14)  of  a  wave  propagating  in  the  random  medium.  By 
solving  the  dispersion  equation,  we  can  obtain  an  effective  propagation  constant  of  the  random  medium. 
However,  this  is  a  forbidden  task  and  can  only  be  done  numerically  when  the  size  of  the  scatterers  is  large 
compared  with  the  incident  wavelength. 

Fortunately,  for  most  microwave  composite  materials,  the  scatterers  are  much  smaller  than  a  millimeter 
in  size.  Therefore,  at  frequencies  in  the  gigahertz  range,  the  wavelength  in  the  matrix  material  is  much  larger 
than  the  size  a  of  particles  and  make  the  nondimensional  frequency  ka  fall  in  the  Rayleigh  region  .  Thus,  one 
can  solve  the  dispersion  equation  in  the  long  wavelength  regime  and  obtain  a  closed  form  solution. 

Retaining  only  the  dipole  terms,  we  have  the  following  dispersion  equation  [1] 


T1  (JHn  +  JH-  /  2)  —  1 


3  JHj  r/2 


3JH,T  /2 


TZ(JH0  +  JH2/2)-l 


JH0  =  3jQ/X2  +  w 
JH,  =  3  j  Q  (K/k2)/X^ 


JH,  =  3jQ(K/k  .)2IX] 


Q=c/{l-(K/k  2n 
T,=j(2x5/3)[(ji1-H2)/(2ii2  +  Hi)j 
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12=j(2X^/3)[(e1-e2)/(2£2  +  £1)] 

w  =  (l-c)4/(l  +2c)2 
Xj  =  Re  (k2a) 


i 


£,=e;+j£i 


n2=  VJ42 


£ 

2 


k2  =  tO(^£2)1/9-/Co 


K  =  Effective  propagation  constant  =  CD  (<p>  <£>)1/2  /  cq 
cq  =  Light  speed  in  vacuum 


c  =  Volume  fraction  of  scattered 

The  effective  dielectric  constant  <  e  >  and  permeability  <  jx  >  can  be  easily  obtained  from  Eq.  (16)  due 
the  fact  that  T1  and  T2  are  uncoupled  and  the  effective  propagation  constant  K  is  polarization  independent 
the  Rayleigh  region.  This  approach  can  also  be  applied  to  the  study  of  electromagnetic  scattering  by  a 
periodic  array  of  particles. 

For  randomly  distributed  lossy/lossless  scatterers  with  a  high  concentration  in  a  lossy/lossless  matrix, 
the  effective  permeability  and  permittivity  are,  respectively, 

<  H  >/  H2  =  {  1  -  [pw  -  (cy— cj)/(l  +  y2)]  (B  +  jA)}  /  {  1  -  [pw  +  (cy-cj) /2(1  -  y2)]  (B  +  jA)} 


<£>/£  ={ 
2 

where. 


1  -  [Pw  -  (cy-  cj)/(l  +  y2)]  (C  +  jD)}  /  {  1  -  [pw  +  (cy-  cj)  /  2(1  -  y2)]  (C  +  jD) } 

(17) 


A=2(By+ A)  B  =  2  (Ay- B) 

C  =  2  (Cy-  D)  D  =  2(Dy+  C) 

A=  [Oi'  - if) (2\l'2  +  |i' )  +  Oi'  -  if) (2 n'  +  M-j)] / A 

B  =  [(if  -  if) (2\i’2  +  ji')  +  Of  -  if )  (2 if  +  if)]  /  A 

C  =  [(£;  -f)(2£2  + c' )  +  (£j  “  £j)  (2£2  +  e[)J 1 A 

D  =  [(f  -  f)  (2e^  +  e  ) + (e  -  f)  (2e"  +  e')]  /  A' 
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A  =  (2^  +  n'i)U(2tl%|i;)2 
A'-Oej  +  e'^  +  Gc'  +  e')2 
T=(xJ-3XjX;)/(Xj-3XjX^ 
P  =  (Xj  -  3X2  Xj)  /  3 


/ 

=  Im  (k2a) 

4.  RESULTS  AND  DISCUSSION 

Different  from  various  equations,  which  are  quite  popularly  used,  such  as  Maxwell-Gamett  equation, 
Clausius-Mosotti  equation,  and  other  equations  like  Lichtenecker  logarithmic  equation,  our  equations  for 
independent  permittivity  and  permeability  can  be  used  in  either  purely  scattering  domain  or  absorption 
domain.  In  other  words,  the  matrix  as  well  as  the  inclusion  materials  can  both  be  lossy  and  lossless  or 
combinations  of  the  two.  Furthermore,  a  rigorous  multiple  scattering  theory  has  been  employed  in  the 
derivations  accounted  for  high  concentrations  of  the  inclusion  phase  while  none  of  the  mentioned  equations 
has  considered  multiple  scattering. 

In  the  following,  we  discuss  some  limiting  cases  when  (17)  is  applied  to  compute  the  effective 
properties  of  microwave  composite  materials.  Because  the  effective  permittivity  and  permeability  bear  the 
same  form.  For  simplicity,  we  use  only  the  expressions  for  the  effective  permittivity.  However,  all  the 
following  discussion  are  exactly  those  for  the  effective  permeability  when  all  the  e’s  are  replaced  by  the  |Ts 
in  the  statements. 

Case  1.  Pure  Scattering 

If  neither  the  scauerers  nor  the  matrix  material  are  lossy,  i.e.  Ej"  =  £2"  =  0,  only  scattering  contributes 
to  the  imaginary  part  of  the  effective  dielectric  constant  <  e  >. 

<e>/e2  =  (l  +2cT)/(l  -cT) +j  6(JwcT2/(l -cT)2  (18) 

where  T  -  ( t\  -  62'  )  /  ( +  2  £2'  ).  The  effective  dielectric  constant  <  e  >  becomes  £2'  and  Ej’, 
respectively,  when  c  equals  to  0  (no  scatterer)  and  1  (  matrix  material  totally  occupied  by  scatterers). 
Meanwhile,  the  imaginary  part  of  <  e  >  vanishes  due  to  the  disappearance  of  scatterers. 

Case  2.  Lossy  scatterers  in  a  lossless  matrix 

This  is  a  common  case  happened  in  various  problems,  e.g.  suspensions  in  fluids.  For  this  case,  £2"  =  0 
but  Ej"  is  not.  Equation  (17)  can  be  reduced  to  the  form  as  follows 

<e>/e'={  1+2 (Pw+cj)(D-jC)}/ {1  +  2  [pw-cj/2](D-jO)  •  (19) 

Again,  <  e  >  becomes  when  c  =  0.  However,  when  c  =  1, 

<£>/£'  =(l+2C  +  2jD)/(l-C-jD)  =  (£  +)£*)/£  .  (20) 

2  112 

It  is  easy  to  see  that  the  effective  dielectric  constant  becomes  that  of  the  lossy  scatterers. 

Case  3.  Lossless  scatterers  in  a  lossy  matrix 

Although  this  problem  has  seldom  been  considered  for  electromagnetic  case  ,  e.g.  particles  in  a  viscous 
fluid,  it  is  useful  to  check  the  validity  of  the  formalism.  For  this  case,  Ei”  =  0  and  £2"  is  finite.  It  is  easy  to 

Journal  of  Wave-Material  Interaction,  Vol.  3,  No.  3,  July  1988 

247 


Ma,  Varadan  and  Varadan 


Effective  Properties  of  Microwave  Composites 


show  using  (17),  for  c  =  0,  <  e  >  =  E2  =  e2  +  J  e2"  •  For  c  =  1,  (17)  reduces  to 

<£>/£2  =  (l+2C  +  2jD)/(l-C-jD)  =  £j/(e^+j£')  ,  (21) 

which  means  the  effective  dielectric  constant  becomes  that  of  the  scatterers  and  is  real. 

Case  4.  Lossy  scatterers  in  a  lossy  matrix 

For  this  case,  the  complete  form  of  (17)  must  be  considered.  As  for  the  relative  contribution  form 
scattering  and  absorption  towards  <  e  >,  one  can  easily  show  by  using  case  2  that 

<£>/£' =[(1+ 2c  C)  (1  -  c  C)  -  2  c2  D2]  /  [(1  -  c  C)2  +  c2  D2] 

+j[3cD  +  6p  wc(C2  +  D2)]/[(1-cC)2  +  c2D2]  .  (22) 

If  Ei"  =  0  then  D  =  0,  the  above  equation  becomes  that  for  case  1.  Therefore,  the  scattering  contribution  to 
the  imaginary  part  of  <  e  >  has  a  (3  dependence.  Because  P  is  of  the  order  of  (k2a)3,  which  is  much  less  than 
one  for  the  long  wavelength  limit,  we  know  the  scattering  contribution  is  fairly  small  when  compared  with 
the  absorption  contribution  which  is  frequency  independent  if  Ej  and  e2  are  also  frequency  independent.  The 
real  part  of  the  effective  dielectric  constant  is  less  affected  by  the  absorption  due  to  the  fact  that  ej"  is,  in 
general,  fairly  small. 
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ABSTRACT 

Propgation  of  electromagnetic  waves  in  a  parallel-plate  waveguide  wholly  filled  wvh  a  chiral  medium  is 
examined.  The  dispersion  equation  derived  leads  to  two  sets  of  modes.  Propagation  constants  for  the  two  sets 
have  been  numerically  obtained. 

1.  INTRODUCTION 

Although  the  phenomenon  of  chirality  is  known  chiefly  at  the  molecular  level,  it  ha s  been  suggested  [1] 
that  particles  endowed  with  chirality  can  exist  at  even  lower  frequencies,  say,  in  the  GHz  range.  This  is 
because  chirality,  or  handedness,  is  a  geometric  property:  for  example,  the  electromagnetic  response  of  a 
right-handed  helix  is  different  from  that  of  a  left-handed  one  [2].  Furthermore,  by  embedding  such  chiral 
particles  in  a  low-loss  dielectric  medium,  that  medium  too  will  possess  handedness.  With  advances  in 
polymer  science,  it  is  becoming  increasingly  possible  that  such  artificial  materials  can  be  manufactured  with 
ease,  and  their  properties  tailored  by  altering  the  sizes  and  concentration  of  the  embedded  chiral  particles. 

Significant  advances  have  taken  place  recently  in  the  formulation  of  a  frequency-domain  electromagnetic 
theory  for  chiral  media,  and  these  have  been  summarised  by  us  elsewhere  [1].  With  these  developments,  it  is 
time  that  aspects  relating  to  the  application  of  chiral  media  for  practical  problems  be  explored.  To  that  end, 
we  have  already  obtained  the  eigenmodes  of  a  perfectly  conducting  sphere  filled  with  a  homogeneous, 
isotropic,  chiral  medium  [3].  In  continuation  of  our  aim,  we  study  here  the  modes  of  a  parallel-plate 
waveguide  filled  with  a  chiral  material.  It  is  our  conjecture  that  this  geometry  will  be  of  use  in  the 
development  of  integrated  circuitry  with  chiral  substrates. 

Consider  a  source-free  region  occupied  by  an  isotropic  chiral  medium  in  which  the  usual 
constitutive  relations  D  =  eE  and  B  =  pH  do  not  hold  due  to  their  incompatibility  with  the  handedness  of 
the  medium.  Instead,  the  relations 

D=eE  +  (teVxE  ;  B  =  nH  +  (5^VxH  (1) 

hold,  and  satisfy  the  requirements  of  time-reversal  symmetry  and  reciprocity.  Following  [4],  the 
electromagnetic  field  is  transformed  to 

^  =  ^l+aR^2  ’  ^=^2  +  aL^l  ’ 

where  the  left-  (LCP)  and  the  right-  (RCP)  circularly  polarized  fields,  Qj  and  Q2,  respectively,  must  satisfy 
the  Helmholtz  equations 

(V2  +  yj)  Q,  =0  ;  {V2  +  i?)Q2  =  0  .  (3) 

along  with  the  rotational  conditions 

vxQ,=y,Q,  ;  vxq2=-y2q2  .  (4) 

Needless  to  say,  these  fields  are  also  divergence-free,  vide 

V’Q,  =  0  ;  V.Q2  =  0  .  (5) 
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In  these  equations,  the  two  wavenumbers  are  given  by 

Y,  =k/ {1  —  kp}  ;  ^  =  k/{l  +  kP)  ,  (6) 

and 

aL  =  -i(e/p)1'2  ;  aR  =  -i(n/E),/2  .  (7) 

An  exp[-ioat]  time  dependence  has  been  assumed  throughout  this  work,  while  k  =  ©(pe}*/2  is  simply  a 
shorthand  notation. 

2.  MODAL  ANALYSIS 

Consider  the  region  bounded  by  the  perfectly  conducting  plates  z  =  ±d,  which  is  wholly  filled  with  a 
chiral  medium.  It  is  well-known  that  although  LCP  and  RCP  waves  can  propagate  independently  in  an 
unbounded  chiral  region,  at  a  bimaterial  interface  mode-conversion  occurs  [4],  Therefore,  pure  LCP  or  RCP 
modes  cannot  exist  within  a  bounded  chiral  volume.  Using  a  representation  given  by  us  earlier  [S],  the 
electromagnetic  field  inside  the  parallel-plate  waveguide  can  be  adequately  expressed  by 

Qi  *  IA+  Yl"1  {-“i  ex  +  K  ez  _  ty  ey )  «p[i<XjZ] 

+  Ai_  Yi 1  {<*i  ex  +  <  ez  ~  iYt  ey }  exp[-ictj z]]  exp  [iKx]  ,  (8a) 

%  =  (-°2  ex  +  K  ez  +  ^2  ey>  CXptiCCjZ] 


+  A2_y'1  {a2  ex  +  k  e2  +  iy2  e  }  exp[-ia2z]]  exp  [iicx] 


(8b) 


which  satisfy  the  phase-matching  conditions  via  the  horizontal  wavenumber  k.  In  these  equations, 
dj  =  +V(Yi2  -  k2)  and  ot2  =  +'!('$'  -  K2),  while  Aj-j.  and  A2±  are  the  unknown  field  coefficients. 

The  boundary  conditions  require  that  the  tangential  components  of  the  electric  field  be  identically  zero  on 
the  surfaces  z  =  ±d.  The  use  of  (2)  and  (8)  along  with  the  boundary  conditions  leads  to  the  dispersion  equation 

0=  [{“jYj1  +  a2^1}  {1  -exPti2  (Yj  +Y2)  d]}  +  {ajj1  -  c^y"1}  {exp[i2yj  d]}  -exp[i2y2 d]}]  x 


[{aiV  +(hh}  U-exp[i2  (Yj  +  Y2)d]}  -  {aj'1  -a^1}  {exp[i2y,  d]}  -exp[i2y2d]}]. 

(9) 

Since  (9)  contains  two  factors  on  its  right  hand  side,  it  is  clear  that  two  kinds  of  modes  can  exist.  The 
dispersion  equation  for  the  first  set  of  modes  is  given  by 

0=  [{aj'1  +  a2y21}  {1  -expfBty  +y2)d]}  +  {aj"1  -a^j1}  {exp[i2yj  d])  -exp[i2y2d]}]  , 

(10a) 

and  for  these  modes  it  is  easily  seen  that 

A1+/A1_  =  A2+/A2_  =  -1  ,  (10b) 


^A2+/A1+=sin(a1d)/sin(a2d)  .  (10c) 

As  a  consequence  the  field  components  Ex  and  Hx  are  even  with  respect  to  z,  whereas  the  field  components 
Ey,  Hy,  E 2  and  are  odd  with  respect  to  z. 

Cm  the  other  hand,  the  dispersion  equation  for  the  second  set  of  modes  is  given  by 
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0=[{a17;1  +  ct2Y21}  (l-exp(i2(Y1+T2)d]}-{a1Y[l-a2Y21)  {exp[i2Yj  d]} -cxp[i2Y2  d]{]  , 

(11a) 

and  for  these  modes  it  can  be  shown  that 

A1+/Ai_  =  A2+/A2_=1  ,  (lib) 

aRA2+/A1+  =  cos(a1d)/cos(a2d)  .  (11c) 

As  a  consequence  the  field  components  Ex  and  Hx  are  odd  with  respect  to  z,  whereas  the  field  components 
Ey,  Hy,  Ez  and  Hz  are  even  with  respect  to  z. 


0.0  2.0  4.0  6.0  8.0  10.0  12.0 

Normalized  Frequency  (k*d) 

Figure  1.  Solutions  Kd  of  the  dispersion  equations  (12c)  and  (13c)  are  virtually  identical  for  p/d  £  10-4. 

They  also  correspond  (almost  exactly)  to  the  TE-  and  TM-  polarised  fields  when  p/d  =  0. 

3.  NUMERICAL  RESULTS 

With  the  developments  of  the  previous  section,  the  modal  fields  of  the  first  set  can  be  compactly  set 
down  as 

Q,  =exp[itcx]Y^  {-<*5  cos(ajZ)ex  +  iKsin(a,z)  ez  +YIsin(a,z)  ey)  ,  (12a) 

cxp[iKx]  y2'  {-a2cos(02z)  ex  +  iKsin(a2z)  ez -  Y2sin(a2z)  ey }  sin(ajd)  / sin(a2d)  ; 

(12b) 

the  horizontal  wavenumber  k  can  be  determined  for  this  set  by  solving  the  equation 
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aly2/a2y\+m(ald)/t2n(a1d)  =  0  .  (12c) 

In  a  similar  vein,  the  modal  fields  of  the  second  set  can  be  written  as 

Q,  =exp[iicx] y”1  {-ctjSin^z) ex  +  Kcos(a,z) ez -iy1cos(a,z)  ey}  ,  (13a) 

\  Q2  =  cxp[iKx]  y"1  {-ictjS^c^z)  ex  +  kcos^z)  ez  +  iy2cos(a2z)  ey )  cos(ajd)  / cos(a2d)  ; 

(13b) 

the  horizontal  wavenumber  k  can  be  determined  far  this  set  by  solving  the  equation 


“lV0^!  +cot(aid)/cot(a2d>  =  0  .  (13c) 

The  solutions  Kd  of  the  dispersion  equations  (12c)  and  (13c)  were  obtained  numerically  on  a  Macintosh 
II  minicomputer  as  functions  of  the  normalized  frequency  kd  <  10.0  for  various  values  of  p/d;  while  k  < 
min(yj,Y2].  was  keP£  Jess  than  0.5  for  Figs.  1-3,  it  is  less  than  0.99  for  Figs.  4  and  5.  Fig.l  shows  the 
calculations  for  3/d  =  10'^.  For  P/d  S  10"^,  the  roots  for  the  two  sets  of  modes  did  not  appreciably  differ 
from  each  other;  furthermore  o^d  and  a^i  were  approximately  equal  to  integral  multiples  of  n/2.  The  effect 
of  chirality  became  numerically  appreciable,  however,  at  P/d  =  10‘3,  although  the  differences  between  the 
roots  for  the  two  sets  are  still  small  enough  to  be  appreciated  on  a  graph. 

Shown  in  Figs.  2  and  3  are  the  roots  Kd  of  the  dispersion  equations  (12c)  and  (13c),  respectively,  for  p/d 
=  10*2.  The  root  structures  are  now  different  for  the  two  sets,  and  particularly  so  when  kd  is  high.  When  p/d 
increases  even  further,  the  differences  are  even  more  telling,  as  illustrated  in  Figs.  4  and  5.  Thus  it  is  only  at 
the  higher  frequencies,  and  for  a  higher  degree  of  chirality  as  characterised  by  larger  values  of  |P|,  that  the 
effect  of  the  chirality  of  the  medium  becomes  significant. 


0.0  2.0  4.0  6.0  8.0  10.0  12.0 


Normalized  Frequency  (k*d) 

Figure  2.  Solutions  Kd  of  the  dispersion  equation  (12c)  of  the  first  set  at  p/d  =  10*2;  kp  S  0.5. 
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Normalized  Frequency  (k*d) 

ire  3.  Solutions  icd  of  the  dispersion  equation  (13c)  of  the  second  set  at  (3/d  =  10*2;  k(J  <,  0.5. 


Normalized  Frequency  (k*d) 

Figure  4.  Solutions  icd  of  the  dispersion  equation  (12c)  of  the  rust  set  at  (i/d  *  10* k(3  £  0.99. 
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Figure  5.  Solutions  ted  of  the  dispersion  equation  (13c)  of  the  second  set  at  p/d  =  10'1;  kp  <  0.99. 
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Errata: 
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Equations  (9),  (10a),  (11a)  and  (13a)  of  the  subject  paper  [1]  contained  typographical  errors.  The  correct 
versions  are  as  follows: 


0=  [{a  Yj1  +  a^1}  {l-exp^a^  a2)d]}  +  {a^1  -a^1}  {exp[i2a  d]}  - exp[i2a2d] }]x 
[{a^1  +  a^1}  {1  - exp[i2(ai  +  a^d]}  -  {a^1  -  a^1}  {expfttc^d]}  - expfi2a2d] }] . 

(9) 

0  =  [{“jYj1  +  a^1)  (1  -exp[i2  (<Xj+  a2)d]}  +  {a^1  -  a^1 }  {exp^ajd]}  —  exp[i2ct2d] } ] , 

(10a) 

0  =  [{c^Yi1  +  1 1  - exp[i2  (ccj+a^d]}  -  {a^1  -  a^1}  {exp[i2a,d]}  -exp[i2a2d]}], 

(11a) 


Qj  =  exp[iKx]  y/  {-iajSi^OjZ)  ex  +  tccos(ajZ)  ez  -  iy1  costc^z) ey) ,  (13a) 
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Abstracl—Tht  coupled  dipole  approximation  method  has  been  ex¬ 
tended  in  order  to  compute  the  scattering  characteristics  of  three- 
dimensional,  homogeneous,  lossless,  anisotropic  objects.  Numerical 
results  are  given  for  uniaxial  spheres  made  of  titanium  dioxide. 

Introduction  and  Theory 

The  Mie  formulation  [1]  of  scattering  by  a  homogeneous,  isotropic 
sphere  has  enjoyed  an  enormous  amount  of  popularity  ever  since  its 
inception  in  1909,  but  no  formulation  of  comparable  magnitude  exists 
for  anisotropic  spheres.  Recently,  Graglia  and  Uslenghi  [2]  have 
formulated  integral  equations  for  scattering  by  anisotropic  objects 
using  microscopic  arguments;  however,  they  have  been  successful  in 
the  numerical  implementation  of  their  formulation  only  for  two- 
dimensional  problems  [3].  A  method  for  the  solution  of  electromag¬ 
netic  scattering  by  transversely  isotropic  cylinders  due  to  Monzon 
and  Damaskos  [4]  should  also  be  mentioned,  as  should  Moon’s 
formalism  [S]  for  the  scattering  of  elastic  waves  by  uniaxial, 
piezoelectric  cylinders  made  of  barium  titanate. 

The  oscillations  of  a  small  gyrotropic  sphere  have  been  handled  by 
Walker  [6]  by  considering  quasistatic  electric  and  magnetic  modes. 
For  the  electric  modes  the  internal  E-field  was  expressed  as  V^,,  and 
for  the  magnetic  modes  the  internal  H-field  was  modeled  as  Vi^m. 

Manuscript  received  March  4,  1988;  revised  August  12,  1988. 

The  authors  are  with  the  Department  of  Engineering  Science  and 
Mechanics  and  the  Center  for  the  Engineering  of  Electronic  and  Acoustic 
Materials.  Pennsylvania  State  University,  University  Park,  PA  16802. 

IEEE  Log  Number  8927232. 


Walker’s  method  was  iteratively  refined  by  Yan'-Shen  [7],  [8],  but 
remained  valid  for  very  small  spheres.  This  work  was  extended  by 
Pistol’kors  [9]  who  considered  the  self-oscillation  modes  with  no 
azimuthal  dependence  of  a  small  ferrite  sphere.  Further  research 
using  similar  approximations  has  also  been  reported  in  two  papers  by 
Wolff  [10],  [1 1].  A  more  comprehensive  approach  [12],  published  in 
Russian,  on  scattering  by  ferrite  spheres  appears  to  be  flawed  in  its 
selection  of  the  basis  functions  for  representing  fields  in  the  spherical 
coordinate  system. 

A  semimicroscopic  approach  [15],  in  which  the  scatterer  is 
modeled  by  mutually  interacting  point-polarizable  groups,  has  gained 
popularity  in  the  last  few  years  (e.g.,  [16],  [17]).  This  technique, 
dubbed  the  coupled  dipole  approximation  procedure,  has  been 
successfully  applied  to  compute  the  scattering  characteristics  of 
nonspherical,  lossless  dielectric  objects.  This  procedure  has  also  been 
successful  in  predicting  the  Perrin  matrix  for  scattering  by  bacteria 
[18].  Here  we  report,  for  the  first  time,  the  application  of  this  method 
for  the  computation  of  the  scattering  by  homogeneous,  anisotropic 
scatterers  with  a  relative  permittivity  dyadic  «.  Numerical  results  will 
be  presented  for  spheres  made  of  titanium  dioxide,  which  is  an 
uniaxial  medium. 

Detailed  discussions  of  the  presented  method  as  applicable  to 
isotropic  scatterers  are  available  elsewhere  [15]— [17];  therefore,  here 
we  shall  mention  only  its  features  arising  due  to  the  anisotropic 
properties  of  the  scatterer.  The  scatterer  is  modeled  by  a  three- 
dimensional  array  of  point-polarizable  groups  arranged  on  a  cartesian 
lattice  embedded  in  free  space.  Each  one  of  these  groups  is 
characterised  by  a  polarizability  tensor  p,  such  that  when  an  electric 
field  E  hits  it,  the  scattered  field  can  be  adequately  expressed  as  being 
due  to  an  electric  dipole  moment  p,  with 

P  =  P  ’  E.  (1) 

It  has  been  shown  by  Bedeaux  and  Mazur  [19]  that  the  polarizability 
of  a  small  volume  u  and  dielectric  dyadic  (e0e)  immersed  in  free  space 
is  given  by 

f-3ueo[«-S]  •  [.  +  23]-'.  (2) 

For  our  purposes  we  set  u  =  V/N,  in  which  V  is  the  volume  of  the 
scatterer  and  N  is  the  total  number  of  point-polarizable  groups  it  has 
been  modeled  by.  It  is  to  be  emphasized  that  from  a  geometric 
viewpoint,  the  polarizability  p  could  be  that  of  an  equivoluminal 
sphere.  Still,  for  our  purposes,  each  sphere  is  to  be  visualized  solely 
as  a  point-polarizable  group  located  at  the  center  of  the  sphere;  hence, 
any  discussion  of  intersecting  spheres  is  moot.  It  is  hoped  that  with  a 
sufficiently  large  N,  a  convergent  result  would  be  obtained  (as  is  the 
case  with  the  method  of  moments),  and  which  would  be  checked  with 
experimental  results  for  adequacy  and  accuracy. 

It  also  needs  be  mentioned  that  the  term  [<  -  £J]  •  [«  +  2£]~ 1  has 
not  been  obtained  through  a  simple-minded  extension  of  the 
Clausius-Mossotti  relation  for  small,  isotropic,  dielectric  spheres 
[20].  Instead,  Bedeaux  and  Mazur  [19]  invoked  a  statistical  theoiy  for 
the  wave  vector  and  the  frequency-dependent  dielectric  tensor  for 
their  treatment  of  a  nonpolar  fluid.  Hence,  their  theory  is  not 
restricted  to  small  values  of  the  polarizability  per  unit  volume.  The 
only  restriction  here  is  that  the  medium  not  be  lossy,  for  it  is  not  dear 
if  a  lossy  spherical  unit  can  be  adequately  represented  by  a  complex 
electric  dipole. 

The  electric  field,  E^,^,  exciting  the  mth  spherical  subunit  is 
composed  of  the  incident  field  E;«  (I’m)  us  well  as  the  scattered  fields 
due  to  all  the  other  subunits  evaluated  at  r„,  the  location  of  the  mth 
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subunit.  Thus 

Einc(f/>i) =  2  ®Li«»  '  =  1,  2,  ■  ■  • ,  M,  (3a) 

»-IA-  -M 

where 

Q.«<i  =  -  (1  -  V4«0/?a,) 

•  exp  [ikRm„]{gm„ P 

*/n/r  ®jiw»  ^/wr  p  }■  (3b) 

In  these  equations,  Sm„  is  the  Kronecker  delta,  k  =  tovfio^ol  is  the 
free-space  wavenumber,  Rm„  =  rm  -  r„,  em„  =  Rmn/Rmn,  gmn  = 
3(Ar/?m„]~J  -  li[kRmn\-'  -  1  and  3 hmn  =  g„„  -  2.  Equation  (4b) 
can  be  soi  for  the  exciting  fields,  whence  the  far-zone  scattered 
field  due  to  the  anisotropic  scatterer  can  be  easily  derived  as 

Fsc(/')=  lim  {kr  exp  [-/Arr]E*(r)} 

kr~<*> 

=  (k3/4rt0)  2  exp  l-ikrm  ■  r/r]f 

m*  1,2,-  •  •  M 


Fig.  1.  Comparison  of  the  scattering  patterns  of  an  isotropic  sphere,  of 

radius  a  and  relative  permittivity  7.2,  using  Mie  tneory  ( - )  and  the 

present  method  ( - ).  The  plane  wave  Einc  =  e, exp  [ikz]  is  incident  with 

ko  =  0.5. 


•  [/-/"2rrl  •  .  (4) 

Equations  (3)  give  the  procedure  its  name.  The  exciting  field  is 
determined  at  the  location  of  each  point-polarizable  group  via  the 
second  term  in  (3b),  which  permits  interactions  between  the  groups. 
These  interactions  are  exact,  given  that  the  groups  are  dipolar. 

Numerical  Results  and  Discussion 

Equations  (3b)  and  (4)  were  programmed  and  solved,  for  a  uniaxial 
sphere,  on  a  DEC  VAX  1 1/730  minicomputer  using  single  precision 
complex  arithmetic.  Uniaxial  dielectric  media  are  characterized  by  an 
anisotropic  permittivity,  with  the  relative  permittivity  stated  in  dyadic 
notation  as 


t  =  ei^|  +  («|,-«i)cc,  (5) 

where  is  the  unit  dyadic  and  c  is  the  unit  vector  along  the  optic  axis; 
without  loss  of  generality  it  can  be  assumed  that  c  =  e*.  We  chose  the 
scatterer  to  be  a  sphere  of  radius  a  (i.e.,  V  =  4ra}/3)  made  up  of 
titanium  dioxide  (e±  =  0.913,  =  7.197)  on  which  plane  waves 

are  incident  at  a  normalized  frequency  ka  =  0.5.  We  compared  the 
scattering  pattern  F*  for  a  dielectric  sphere  of  relative  permittivity 
7.2  with  rigorous  calculations  made  using  the  Mie  theory,  and 
observed  agreement  within  5  percent  when  179  spherical  subunits 
were  utilized;  see  Fig.  1.  For  the  results  presented  in  the  sequel, 
therefore,  179  spherical  subunits  were  used  to  model  the  uniaxial 
spheres. 

Shown  in  Fig.  2  are  computations  made  for  the  titanium  dioxide 
sphere  when  a  plane  wave,  Einc  =  ex  exp  [  ikz  ]  is  incident  on  it.  The 
scattering  pattern  in  the  x-z  plane  turns  out  be  predominantly  9- 
polarized,  whereas  in  the  y-z  plane  it  is  ¥>-polarized.  For  comparison, 
the  data  from  Mie  theory  for  isotropic  spheres  of  relative  permittivi¬ 
ties  5.913  and  7. 197  are  also  plotted.  The  depolarized  scattered  fields 
(i.e.,  the  <?- polarized  one  in  the  x-z  plane  and  the  0-polarized  one  in 
the  y-Z  plane),  which  are  due  to  anisotropy  of  the  scatterer,  have  not 
been  plotted  since  they  are  several  orders  of  magnitude  smaller  for 
the  calculations  pertaining  Fig.  2. 

In  Fig.  3,  the  scattering  pattern  is  plotted  when  the  incident  plane 
wave  is  given  by  E,„.  =  es  exp  [ikx\\  the  depolarized  component  of 
the  scattered  field  is  ^-polarized  in  the  x-y  plane,  and  0-polarized  in 
the  y-z  plane.  The  scattering  pattern  in  the  x-z  plane  is  still  0- 
polarized;  but,  in  the  y-z  plane,  the  predominant  component  of  the 
scattering  pattern  is  also  the  0-polarized  one,  not  the  ^-polarized  one. 
This  suggests  the  tremendous  effect  of  the  anisotropy  when  the 


Fig.  2.  Scattering  pattern  of  a  TiOj  sphere  ( - )  of  radius  a  on  which  a 

planewave,  E«*  =  e,  exp  [ikz]  is  incident  with  ka  =  0.5.  Top:  e.-F*  in 
the  y-z  plane.  Bottom:  e*’F„  in  the  x-z  plane.  For  comparison.  Mie  theory 
calculations  for  isotropic  spheres  of  relative  permittivities  5.913  (  -  •  -  • 
-  •)  and  7.197  ( - )  are  also  plotted. 

incident  electric  vector  is  parallel  to  the  optic  ax's.  This  understand¬ 
ing  is  further  underscored  in  Fig.  4  where  the  scattering  pattern  is 
plotted  for  the  incident  plane  wave  Einc  =  ey  exp  [ikx];  the 
depolarized  component  of  the  scattered  field  is  0-polarized  in  the  x-z 
plane,  and  ys-polarized  in  the  y-z  plane.  It  should  be  noticed  that  the 
depolarized  component  manifests  itself  significantly  only  in  the  y-z 
plane,  but  its  magnitude  is  small. 

This  coupled  dipole  approximation  method  is  simple,  and  does  not 
require  the  solution  of  any  integrodifferential  equation.  As  is 
apparent  from  (3),  the  presented  method  will  work  for  any 
constitutive  dyadic  c,  provided  it  is  lossless;  furthermore,  scatterers 
of  any  shape  can  be  accommodated.  Thus,  this  procedure  will  work 
well  for  any  three-dimensional,  lossless,  anisotropic  scatterer.  With 
increasing  frequency,  however,  the  number  of  spherical  subunits 


9  in  che  x-z  or  the  y-z  plane 

Fig.  3.  Scattering  pattern  of  Ti02  sphere  of  radius  a  on  which  a  plane  wave, 
Einc  =  exp  [ikx]  is  incident  with  ka  =  0.5.  Code:  e^-F^  in  the  x-z  plane 
( - ),  e^F*  in  the  y-z  plane  ( - ),  and  in  the  y-z  plane  ( - ). 
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Fig.  4.  Scattering  pattern  of  a  Ti02  sphere  of  radius  a  on  which  a  plane 

wave,  Ei„e  =  e,  exp  [ifcjr]  is  incident  with  ka  =  0.5.  Code:  e,-FK  ( - ), 

e,-F*  in  the  y-z  plane  ( - ),  and  e,  F^  in  the  y-z  plane  ( - ). 


required  for  a  convergent  solution  increases,  thereby  straining  the 
memory  resources  of  any  computer.  On  the  other  hand,  for  asphericai 
objects,  fewer  spherical  subunits  will  be  needed,  making  this 
procedure  very  attractive.  The  danger  of  having  to  invert  extremely 
large  matrices  can  be  surmounted  by  resorting  to  an  order  of 
scattering  procedure,  in  which  the  interactions  are  built  in  through  a 
Feynmann  series,  as  has  been  done  for  isotropic  particles  by  Singham 
and  Bohren  [21].  In  addition,  scattering  by  inhomogeneous  objects, 
like  coated  spheres  and  spheroids,  can  also  be  handled  as  has  been 
done  for  the  isotropic  case  [16].  It  should  also  be  noted  that  this 
procedure  works  best  when  the  internal  Field  does  not  vary  rapidly  in 
space. 
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ABSTRACT 

Electromagnetic  waves  can  discriminate  between  objects  of  different  handedness  due  to  their 
transverse  nature,  which  implies  that  the  origin  of  chirality  need  not  necessarily  be  molecular  as 
in  the  case  of  optically  active  media.  Effectively  chiral  composites  may,  therefore,  be  constructed 
by  embedding  chiral  microstructures  in  non-chiral  host  media. 


INTRODUCTION 

Many  organic  molecules  occur  as  stereoisomers  in  enantiomorphic  pairs,  Le„  one  isomer  is 
the  mirror  image  of  the  other  one,  but  the  two  of  them  are  not  congruent  with  each  other.  The 
baas  far  the  difference  in  the  physical  properties  of  the  mirror-conjugates  lies  in  the  handedness 
or  the  chirality  possessed  by  their  molecular  configurations.  When  an  electromagnetic 
disturbance  travels  through  a  medium  consisting  of  chiral  molecules,  it  is  forced  to  adapt  to  the 
handedness  of  the  molecules.  Electromagnetic  (EM)  waves  can  recognize  the  handedness  of  a 
chiral  object  primarily  due  to  their  transverse  nature,  i.e.,  the  applicable  vector,  infinite-medium 
Green's  function  contains  the  vector  distance  between  the  source  and  the  field  points.  From  a 
microscopic,  not  molecular,  point  of  view,  therefore,  it  is  possible  that  an  effectively  chiral 
composite  medium  may  be  constructed  by  embedding  macroscopic  chiral  objects,  all  of  the  same 
handedness,  in  a  dielectric  host  The  microstiucture  size  should  be  large  enough  that  the  EM 
wave  in  the  matrix  can  appreciate  the  handedness  of  the  micros tructure;  at  the  same  time,  the 
microstructure  size  should  be  small  enough  so  that  the  composite  medium,  consisting  of  achiral 
phases,  is  effectively  chiral.  By  varying  the  concentrations  and  the  sizes  of  the  chiral  inclusions, 
the  properties  c£  the  composite  medium  may  be  altered  to  suit  desired  polarization  characteristics. 
It  has  been  observed  by  us  that  the  introduction  of  chirality  In  an  otherwise  achiral,  dielectric 
scattering  volume  radically  alters  its  scattering  and  absorption  characteristics  [1].  Particularly 
with  the  proliferation  of  research  on  novel  polymers,  it  is  possible  that  such  mati-riult  can  turn 
out  to  be  of  considerable  importance  for  many  electromagnetic  applications  It  has  been 
demonstrated  [23]  that  it  is  possible  to  produce  polymers  which  possess  no  true  axisynunetric 
centers,  and  helical  strands  —  either  left-  or  right-handed  -  can  be  produced. 


WHAT  IS  CHIRALITY? 

In  order  to  understand  the  characteristics  of  chiral  media,  one  has  to  begin  with  the  notions 
of  the  intrinsic  polarization  and  the  intrinsic  magnetization  vectors.  Free  space  is  vacuum.  i.c.,  it 
does  not  contain  any  matter.  The  constitutive  relations  for  the  electromagnetic  fields  in  free  space 
can  be  adequately  expressed  in  the  form 


D-EqE;  B-p0H.  (I) 

Fields  in  material  media  induce  bound  charges  too.  Collectively,  the  constitutive  equations  in 
material  media  have  the  form 
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D  -  CqE  +  P;  B  *  m,H  +  M, 


(2) 


in  which  P  and  M,  respectively,  are  the  polarization  and  the  magnetization.  The  polarization  P  is 
nothing  but  the  avenge  electric  dipole  moment  per  unit  volume,  while  M  is  the  avenge  magnetic 
dipole  moment  per  unit  volume. 

In  non-magnetic  dielectric  media,  a  prominent  contribudon  to  P  comes  from  the  charge 
separation  in  the  atom,  and  is  designated  as  electronic  polarization.  This  is  because  an  applied 
electric  field  causes  displacement  of  the  electrons  with  respect  to  the  atomic  nucleus.  Atomic  or 
ionic  polarization  may  be  caused  by  the  displacement  of  atoms  or  ions  in  a  molecule  by  electric 
fields.  Some  media  may  also  have  permanent  electric  dipole  moments  because  the  centers  of  the 
negative  and  the  positive  charges  axe  not  co-incident  Without  an  applied  electric  field,  these 
permanent  dipoles  are  randomly  aligned,  and  their  effect  is  nulled  on  the  average;  however,  when 
an  electric  field  is  applied,  then  these  dipole  moments  line  up  to  give  rise  to  orientational 
polarization.  In  linear  media,  these  various  contributions  can  be  summed  up  as 


P-XeCoE: 


D»£oIl+XeJE*fiE  (3) 


X*  being  the  electric  susceptibility,  and  e  the  (real)  dielectric  constant 

Each  of  the  three  types  of  polarizations  are  frequency-dependent;  however,  as  the  frequency 
Increases  the  polarizadon  may  not  be  able  to  follow  the  oscillations  of  the  applied  field. 
Orientational  polarizadon  concerns  itself  with  the  most  massive  entities;  hence,  it  is  the  first  to  be 
unable  to  reach  its  equilibrium  value.  The  loss  of  orientational  polarization  occurs  at  lower 
frequencies  for  macnxnoleculcs  and  large-chain  polymers,  and  it  occurs  at  higher  frequencies  for 
small  molecules.  As  the  frequency  increases  even  further,  fust  the  atomic,  and  then  the 
electronic,  polarizations  also  fall  off.  Typically,  the  times  for  orientational,  atomic,  and  electronic 
polarizations  to  reach  their  equilibrium  values  are  10*9  s,  10*t3  s  and  10*15  s,  respectively. 
Highly  elongated  and  aligned  molecules  tend  to  give  rise  to  anisotropic  properties. 

Electronic  polarizability  is  due  to  the  bound  electrons.  Therefore,  it  increases  when  the 
atomic  radius  also  increases.  This  implies  that  polymers  with  delocalized  electrons  will  have 
larger  electronic  polarizations;  for  example,  conducting  polymers  with  electron-withdrawing 
groups.  On  the  other  hand,  polymers  possessing  highly  asymmetric  positive  and  negative  charge 
densities  will  have  large  orientational  polarizabilities.  Too  some  extent,  these  considerations 

permit  one  to  tailor  the  frequency-dependent  dielectric  constants  of  polymers. 

Under  the  influence  of  the  externally  applied  field,  the  polar  molecules  rotate  towards  an 
equilibrium  distribution.  If  the  polar  molecules  ate  massive,  or  if  the  frequency  is  very  high,  the 
rotatory  motion  lags  behind  the  applied  field  and  equilibrium  is  never  attained.  Thus,  the 
polarization  is  no  longer  in  sync  with  the  applied  field,  giving  rise  to  a  conduction  current 
density.  This  results  in  thermal  dissipation  of  energy.  Thus,  ohmic  losses  in  dielectric  nHii  give 
rise  to  a  conductivity  o,  which  is  incorporated  in  the  time-harmonic  constitutive  equations  by 
making  e  complex.  Le. 


D-eE. 


e*col1+Xe]+j«/«. 


W 


00  “  2xf  being  the  circular  frequency.  It  would  also  follow  that  conduction  losses  increase  with 
temperature,  since  the  conduction  charges  have  enhance  mobilities  at  higher  temperatures. 
Furthermore,  mere  massive  molecules  are  less  able  to  follow  the  applied  field,  which  means  that 
the  ratio  tanfi  -  Im[e}/Re{e)  is  generally  higher  for  such  substances.  Unsaturated  polymers  will 
also  have  higher  tano,  because  or  the  increased  availability  of  electrons  for  conduction.  It  should 

also  be  noted  that  while  charge  transport  in  metals  takes  place  via  electrons  only,  charge  transport 

in  polymen  may  even  involve  ions. 

Thereare  microscopic  currents  due  to  electron  spin  and  the  motion  of  electrons  around  the 
nucleus.  These  microscopic  currents  act  as  sources  of  macroscopic  magnetic  fields,  thereby 
endowing  the  electrons  and  the  atoms  with  magnetic  dipole  moments;  but  these  currentsdonot 
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produce  macroscopic  charge  transport.  In  magnetic  media,  the  effect  of  the  magnetic  dipole 
moments  cannot  be  ignored,  giving  rise  to  the  magnetization  M.  Again,  for  linear  media,  it  can 
be  stated  that 


M  =  WoH; 


B-Roll+Xm]H  =  nH 


(5) 


Xm  being  the  magnetic  susceptibility,  and  p  the  (real)  permeability.  Magnetic  ohmic  losses  can  be 
incorporated  by  making  p.  complex  in  the  same  fashion  as  the  electric  ohmic  losses. 

Chirality  was  first  observed  as  optical  activity,  which  is  the  rotation  of  the  plane  of 
polarization  in  certain  linear  isotropic  media.  Pasteur  interpreted  the  observations  of  Biot  and 
Arago  by  imagining  that  the  arrangement  of  atoms  within  an  optically  active  material  is 
mirror-asymmetric  [4].  Phenomenological  studies  by  Drude  [51  indicated  that  the  rotation  of  the 
plane  of  polarization  is  predicted  by  Maxwell's  equations  provide  P  has  an  additional  term 
proportional  to  VxE.  Further  studies  [6,7]  have  shown  that  even  M  must  have  a  term 
proportional  to  VxH.  Thus,  the  constitutive  equations  for  homogeneous  chiral  media  have  the 
form 


D  »  e[E  +  pVxE];  B  =  p[H  +  aVxH],  (6) 


in  which  a  and  P  are  the  chirality  parameters.  Since  the  curl  is  a  pseudo- vector  under  a  reflection 
of  the  coordinate  system,  the  mirror-asymmetry  of  chiral  media  is  immediately  apparent  in  (6). 
The  nonlocal  character  of  (6)  needs  to  be  noticed,  because  the  polarization  P  (resp.  magnetization 
M)  depends  not  only  on  E  (resp.  H)  but  also  on  the  circulation  of  E  (resp.  H).  In  a 
not-too-rigorous  manner,  one  may  even  observe  that  P  (resp.  M)  has  a  component  due  to  the 
time-rate  of  change  of  H  (resp.  E),  vide  Faraday  and  Ampere-Maxwell  laws.  Reciprocity 
demands  that  a  =  p,  a  constraint  which  is  generally  obeyed  by  optically  active  media  [8]. 

Thus,  optically  active  media  (mostly  organic  materials)  are  naturally  chiral,  at  optical 
frequencies;  interesting  examples  include  the  famous  Watson-Crick  double-helix  representation 
of  the  DNA  molecule.  But  the  underlying  principle  behind  chirality  is  the  mirror-asymmetry  of 
the  constituent  microstructure.  Therefore,  artificially  chiral  media  can  be  constructed  by 
embedding  chiral  microstnictures  in  non-chiral  host  media,  such  composites  being  effectively 
chiral  even  at  the  sub-optical,  and  even  in  the  high  mm-wave,  frequencies.  The  microstructure 
size  should  be  large  enough  compared  to  the  wavelength  in  the  matrix  medium  so  that  the  spatial 
variation  of  the  EM  field  can  sense  its  handedness;  at  the  same  time,  it  should  be  small  enough  so 
that,  at  least  in  some  frequency  range,  the  composite  should  appear  to  be  effectively  chiral. 

In  a  non-chiral  medium,  the  simplest  EM  waves  that  can  propagate  are  linearly-polarized 
plane  waves.  But  in  chiral  media  the  simplest  waves  are  either  left-circularly  polarized,  or  they 
are  right-circularly  polarized,  plane  waves;  it  is  important  to  note  that  the  left-  and  the 
right-handed  waves  travel  with  different  phase  velocities  [7].  Furthermore,  when  a  linearly 
polarized  wave  hits  a  chiral  scatterer  embedded  in  a  non-chiral  medium,  the  scattered  field  has 
both  LCP  and  RCP  components.  Now,  let  us  envision  the  propagation  of  a  plane  wave  in  an 
artificial  chiral  medium  which  contains  a  certain  density  of  chiral  inclusions  in  a  host  phase. 
Multiple  scattering  [9]  between  the  chiral  inclusions  will  create  LCP  and  RCP  fields.  In  the  right 
frequency  range  only  LCP  or  RCP  fields  will  be  allowed  to  propagate  in  the  forward  direction, 
thereby  rendering  the  composite  effectively  chiral. 


ANTI-REFLECTION  COATINGS 

There  is  considerable  interest  in  developing  low-weight  coatings  which  will  reduce 
reflections  from  metallic  surfaces.  In  view  of  our  previous  studies  [1],  it  appeared  that  artificial 
chiral  media  may  be  just  the  thing  to  fill  the  bill.  To  further  investigate  this  premise,  we 
considered  the  planewave  reflection  characteristics  of  a  chiral  coating,  of  thickness  d,  on  top  of  a 
perfectly  conducting  plane.  The  chiral  layer  has  a  permittivity  e,  a  permeability  ii  =  |i0,  as  well 


the  chirality  parameter  8.  vide  the  constitutive  relations  [7] 


D  -  e(E  +  JJVxE];  B  =  p0[H  +  pVxH].  (7) 

The  chiral  dielectric  medium  is  low- loss  so  that  its  relative  permittivity  e/e0  is  complex  with 
Imje/eJ  «  Rc(e/e„).  The  thickness  d  of  the  coating  was  assumed  to  be  2  mm  for  all  of  our 
numerical  studies,  and  the  reflection  coefficients,  RTE  and  RTM  [10],  were  computed, 
respectively,  for  TE- and  TM- polarized  plane  wave  incidence  cases. 

A  design  of  a  wideband  anti-reflection  chiral  coating  cannot  afford  to  have  both  e/e0  and  P 
independent  of  frequency.  Once  it  became  clear  from  numerical  experimentation  that  (i)  chirality 
is  ineffective  in  the  reduction  of  the  reflection  efficiency  if  is  purely  real,  and  (ii)  that  either  p 
or  £/£0  or  both  must  be  frequency  dependent  in  order  to  obtain  a  wideband  anti-reflection 
coating,  the  design  of  such  a  coating  became  more  of  an  optimization  problem.  Shown  in  Figs. 
1  -  3  is  a  sample  design  with  £/£g  assumed  to  be  constant,  whereas  p  =  p(f)  assumed  is  shown  in 
Fig.  1.  In  Figs.  2  and  3,  the  enhancement  of  the  absorption  efficiency  of  over  a  30-300  GHz 
frequency  range  by  incorporating  this  frequency  dependent  0  =  (1(f)  is  illustrated,  for  the  TE-  and 
TM-  polarization  incidence  cases,  respectively;  90  is  the  angle  of  incidence  with  respect  to  the 
normal.  The  relative  permittivity  e/e^  is  set  at  5.0  i0.05  for  all  frequencies  considered.  It 

should  be  noted  from  these  two  figures  that  whereas  the  design  objective  of  achieving  Rtc  and 


1,  the  reflection  efficiencies  hover  around  92%  if  P  were  to  be  set  equal  to  zero. 

From  our  numerical  studies,  several  conclusions  can  be  drawn.  Firstly,  in  reducing 
reflection  if  a  larger  constant  p  is  used,  then  the  values  of  eJtg  tend  to  decrease  over  the  entire 
frequency  band  of  interest,  a  goal  which  appears  to  be  desirable  for  a  material  scientist. 
However,  this  also  tends  to  reduce  the  bandwidth  over  which  the  desired  absorption  efficiencies 
can  be  achieved.  Secondly,  chirality  in  the  absence  of  a  lossy  e  is  of  no  use  whatsoever  in 
reducing  reflected  power  density:  conducting  polymers  constitute  feasible  matrix  media. 
Therefore,  chirality  serves  only  as  an  enhancement  factor  for  absorption,  but  of  itself  it  is  not  an 
absorbing  mechanism.  Thirdly,  and  very  importantly,  both  e/e0  and  P  should  be  frequency 
dependent  This  last  conclusion,  however,  was  not  verified  here  because  of  the  complexities  of 
multivariate  optimization  problems. 


EQUIVALENT  DIPOLE  MOMENTS  OF  HELICAL  ENSEMBLES 

Thus,  the  significance  of  chiral  media  cannot  be  denied.  As  has  been  mentioned  earlier, 
effectively  chiral  media  can  be  constructed  by  embedding  chiral  microsmictures  in  non-chiral 
host  media.  Such  microstnictures  may  be  macromolecular  polymers  with  helical  conformations 
[2,3,11].  What  is  of  interest  for  EM  use  is  the  optimization  of  the  polymer’s  chirality:  the 
parameters  governing  the  helix  geometry  must  be  examined  for  their  effect  on  an  incident  EM 
wave.  We  decided  to  model  helical  polymers  by  an  arrangement  of  spherical  beads  (in 
themselves,  large  molecules  for  our  purposes)  suspended  on  a  helical  strand  which  is 
indistinguishable  from  the  surrounding  space.  Whereas  the  spheres  are  also  sufficiently  small  to 
be  modeled  as  point  electric  dipoles,  the  overall  size  of  the  finite  helical  arrangement  can  be  large 
enough  to  be  in  the  high-frequency  regime.  The  helix  on  which  the  tiny  spheres  are  located  is 
given,  in  a  cartesian  co-ordinate  system,  by  the  radius  vector 


r<5)  -  *l«x  cos4  +  uy  h  sin^|  +  uz  %  «  {— , «.),  (g) 


where  a  is  the  radius  and  P  is  the  pitch  of  the  helix;  the  handedness  parameter  h  »  +1  if  the  helix 
curls  up  in  the  +z  direction  according  to  the  right-handed  rule,  and  h  ■  -1,  if  otherwise;  and  u., 
'V  »nd  u,  are  the  unit  orthogonal  vectors.  Let  the  helix  be  finite  in  extent,  having  2N+1 
complete  rotations,  N  being  a  positive  integer  or  zero.  On  each  of  the  2N+1  rings  of  this  finite 
helix,  there  are  2M+1  spheres  arranged  over  equal-A^  segments.  M  being  a  positive  integer 
greater  than  zero.  Each  of  the  spheres  in  this  arrangement  has  a  radius  b  which  is  small  enough 


Figure  1  The  frequency  profile 
P  *  P(f)  used  for  the  results 
shown  in  Figs.  2  and  3.  The 
parameter  p  carries  the  unit  of 
meter. 


Figure  2  Reflection  efficiency 
Rte  plotted  as  a  function  of  0O 
and  frequency  f  for  a  2  mm  thick 
anti-reflection  coating.  The 
complex  permitdvity  e/Eq  =  5.0  + 
i0.05.  In  the  upper  figure  the  3 
*  p(0  of  Fig.  1  is  used  while  in 
the  lower  figure  P  is  set  equal  to 
0.0  m. 
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Figure  3  Same  as  Fig.  2,  but 
Rtm  »*  plotted. 
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that  no  two  of  them  ever  touch.  The  spheres  possess  a  dielectric  constant  e,  and  if  a  Held 
illuminates  the  m*  sphere,  an  electric  dipole  moment  [12] 

pm  *  a  £m  -  4«e0b3  (£  -  £<,)(£  +  2^)- 1  Em  (9) 


is  induced  on  it,  with  a  being  the  isotropic  electric  polarizability  of  the  spheres. 

The  EM  field  incident  on  the  helical  ensemble,  E^,  can  be  any  arbitrary  field  so  long  as  its 
source  is  not  located  anywhere  inside  or  on  the  minimum  sphere  circumscribing  the  helix.  But 
the  field  E„  actually  incident  on  the  m*  sphere  is  not  Ej^r^)  alone;  it  also  consists  of  the  fields 
re-radiated  by  all  of  the  other  spheres  as  well  [9,13].  With  this  reasoning,  the  system  of  3Q 
simultaneous  equations. 


=  Ejnc(rm)  +  ak^rtEo)'1  ^Tuwm  [R  mn"  *  cxp[jkRmn)  * 

*fgmnnmn(nmn*£n)  *  (W)fBmn  *  2)  (10) 


must  be  solved,  in  order  to  obtain  the  various  exciting  fields  £„.  Here,  R^,  ■  lrm  -  rnl,  * 
(rm  -  rnVRmn.  g mn  *  3(kRmn)'2  -  3j(kRmn)*'  - 1,  while  the  wavenumber  k  -  Once  the 

soluuon  of  (10)  has  been  obtained,  the  total  scattered  field  for  lcr  -*  ••  can  be  computed  from 


4neoEsc(r)  -ohV  1  expQkr}  Lm{expl-jkrmT/r]  [^n-rff^nVr2]}.  (11) 


Multiple  scauering  between  the  spheres  would  ensure  that  Em  *  E^r^.  It  is  to  be 
expected  that  the  excidng  fields  should  reflect  the  chirality  of  the  helical  ensemble,  because  the 
ensemble  of  unit  vectors  (n^j  used  in  (10)  is  handed,  even  though  the  scalar  distances  {RmnJ 
do  not  depend  on  h.  This,  indeed  turned  out  to  be  the  case.  From  our  computations  a  very 
interesting  conclusion  can  be  drawn;  Provided  k^  II  us  or  II  ux,  then  (i)  only  uy*Em 
changes  sign  with  h  when  u/E^  »  0,  and  fii)  only  u^xEm  changes  sign  with  h  when  uyx  Eine 
-  0.  Symmetries  of  similar  natures  are  also  present  when  k^  II  Uy.  but  they  are  not  so  easily 
describable. 

Since  the  source  of  EK  is  not  concentrated  at  any  given  location,  it  can  be  conveniently 
expressed  as  being  due  to  a  set  of  multipoles  concentrated  at  the  origin.  The  two  lowest  order 
multipoles  would  be  an  electric  dipole,  p^,  and  a  magnetic  dipole,  m,-,,.  Since  the  helical 
ensembles  are  much  smaller  than  the  ambient  wavelength,  it  appears  that  optimally  chiral  helical 
strands  must  have  maximal  values  of  m^^.  The  magnetic  dipole  moment  can  be  calculated  by 
taking  the  projection  of  EM  on  E*.  in  the  limit  kr  as 


of*  d0  sin0  oP*  dtp  Em(r)  •  Em*(r)  =  of*  d9  sin0  of2*  dtp  E^r)  •  Em*(r),  (12) 


in  which  is  the  electric  field  radiated  by  the  equivalent  magnetic  dipole. 


4j“oEm(r)  ■  *  k2 1**>11'2  r2  cxpljkr]  {rxnteqvt).  (13) 


Several  studies  were  carried  out  for  different  planewave  incidence  conditions  and  the  results  of 
these  investigations  are  summarized  now.  From  studying  the  magnitude  of  m_,  as  a  function  of 
the  ratio  t/P  for  several  different  incident  plane  waves,  we  observe  that  when  art*  -  0.24,  optimal 
chirality  is  achieved  because  m^  is  maximal.  The  computed  optimal  value  of  the  ratio  a/P  is, 
interestingly  enough,  not  far  from  values  reported  for  various  helical  polymers  f  14, 13).  It  is 
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important  to  realize  that  chirality  will  not  be  observed  unless  a/P  *  0  or  P/a  *  0,  because  in  these 
two  respective  limits,  the  ensemble  of  scatterers  is  either  linear  or  circular,  and,  thus,  devoid  of 
handedness.  On  the  other  hand,  the  ratio  a/P  cannot  be  very  large  either,  because  polymers  in  the 
helical  configuration  deviate  only  slightly  from  their  linear  counterpans  [16].  Thus,  that  the 
optimal  value  of  a/P  turns  out  to  be  0.24  is  quite  satisfactory. 

Furthermore,  it  is  obvious  that  by  increasing  a,  m^*,  can  also  be  enhanced:  but  the  spheres 
are  and  a  is  only  weakly  dependent  on  c/Cq.  Even  so,  it  is  recommended  that  a  be  as  large 
as  possible.  Next,  it  was  also  found  that  increasing  the  number  of  spheres  per  ring  (i.e., 
increasing  M)  is  also  helpful  in  increasing  chirality  by  enhancing  the  number  of  dipole-dipole 
interaction  paths.  Short  helices  appeared  to  be,  however,  preferable  to  large  ones;  in  other  words 
N  should  be  as  small  as  possible. 


EXPERIMENTAL  RESULTS 

As  pan  of  our  experimental  effort  to  test  the  viability  of  chiral  composites,  we  devised  a 
free-space  measurement  method,  using  a  dielectric- loaded  horn  antenna  and  the  HP8510A 
network  analyzer  to  determine  the  reflection  from  an  artificially  chiral  slab  which  is  backed  by  a 
metal  plate.  A  41.5  mm  thick  chiral  slab  was  made  of  Eccogel  90,  into  which  were  put  in  a 
certain  number  of  miniature  right-handed  metallic  springs;  the  metal  content  of  the  exposed  area 
due  to  the  introduction  of  the  springs  was  set  at  3.2%  v/v.  The  horn  antenna  was  arranged  so  as 
to  simulate  a  linearly  polarized  planewave  normally  incident  on  the  slab.  The  full  details  of  our 
experimental  efforts  will  appear  elsewhere  [17],  but,  for  the  purpose  of  illustration,  shown  in 
Fig.  4  is  the  reflected  power  measured  over  the  14.5-17.5  GHz  frequency  range.  As  control,  yet 
another  sample  with  3.2%  metallic  content,  but  with  each  spring  replaced  by  an  equi-voluminal 
steel  sphere  was  used;  the  measurements  with  the  control  sample  are  also  shown  in  Fig.  4. 
Almost  over  the  whole  frequency  range,  it  appears  that  the  chiral  sample  reduces  the  reflected 
power  more  than  the  control  sample,  thus  illustrating  the  effect  of  chirality.  It  should  be  noted 
that  the  chiral  slab  is  solid,  in  contrast  to  the  extant  absorbing  coatings  which  are  largely  porous 
and  brittle  (e.g.,  the  lining  of  anechoic  chambers). 


Figure  4  Measured  reflected  power  when  a  linearly  polarized  planewave  is  normally  incident  on 
a  chiral  slab  backed  by  a  metal  plate.  The  identifier  'rh'  denotes  3.2%  v/v  metallic  content  in 
Eocogel  90  in  the  form  of  miniature  right-handed  metallic  springs.  The  control  sample  has  each 

spring  replaced  by  a  steel  sphere  OSS'). 
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